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In  this  dissertation  we  investigate  the  multiple  resource  constrained  single  project 
time/cost  trade-off  problem  for  the  nonpreemptive  case.  Research  in  this  area  has 
concentrated  on  problems  where  the  activity  duration  is  fixed  for  each  mode,  both  in  the 
single  mode  case  and  in  the  multiple  mode  case.  Our  problem,  which  considers  the 
multiple  mode  case,  bridges  the  gap  between  the  traditional  time/cost  trade-off  problem 
and  the  resource  constrained  project  scheduling  problem  by  allowing  for  duration 
reduction  through  higher  direct  cost.  Although  the  backtracking  procedure  proposed  by 
Patterson,  Slowinski,  Talbot  and  Weglarz  and  improved  by  Sprecher  can  solve  the  newly 
defined  problem,  its  unwieldy  computational  requirements,  both  as  an  exact  and  as  a 
heuristic  procedure,  render  it  impractical.  Thus,  we  present  an  efficient  optimal  solution 
procedure  and  a heuristic  procedure.  Our  optimal  solution  procedure  is  a branch  and 
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bound  algorithm.  We  report  computational  results  for  the  optimal  solution  procedure  and 
the  heuristic  procedure  and  compare  these  with  the  backtracking  algorithm.  Our  results 
show  that  our  optimal  solution  procedure  and  heuristic  procedure  perform  better  than  the 
backtracking  algorithm. 


CHAPTER  1 
INTRODUCTION 


We  considered  the  multiple  resource  constrained  single  project  scheduling 
problem.  In  this  problem,  each  activity  is  subject  to  technological  precedence  relations, 
the  execution  of  each  activity  requires  some  set  of  resources,  and  the  amount  of 
resources  available  is  constant.  The  project  has  a due  date.  If  the  project  is  completed 
by  the  due  date,  the  total  project  cost  is  the  sum  of  all  activity  execution  costs,  otherwise 
the  total  project  cost  is  increased  by  the  project  penalty  cost.  The  project  penalty  cost 
is  the  product  of  the  penalty  cost  per  period  by  the  project  delays  (the  total  elapsed 
periods  from  the  project  due  date  to  the  project  completion  time). 

The  Time/Cost  Trade-off  Problem  (TCTP)  and  the  Resource-Constrained 
Time/Cost  Tradeoff  Problem  (RCTCTP)  are  well-known  problems  in  project  scheduling 
problems  that  are  closely  related  to  the  problem  that  this  dissertation  investigates.  The 
differences  among  the  three  problems  should  be  investigated  from  the  underlying 
assumptions. 

The  typical  TCTP  assumes  that  an  activity  duration  can  be  chosen  at  any  duration 
between  its  normal  and  crash  (the  shortest  possible)  durations  and  that  the  activity  cost 
function  is  linear  between  the  normal  and  crash  durations  assuming  a positive  crash  cost. 
The  objective  of  TCTP  is  to  find  the  project  cost  curve  that  specifies  the  minimum 
project  cost  for  all  the  possible  project  durations.  In  TCTP  resource  requirements  for 
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the  execution  of  an  activity  are  not  considered.  TCTP  for  a particular  project  completion 
time  can  be  formulated  as  a linear  program  as  shown  by  Fulkerson  [1961]  and  Kelly 
[1961].  The  complete  project  cost  curve  for  the  entire  project  completion  times  can  be 
found  by  network  flow  algorithms  presented  by  Fulkerson  [1961]  and  Kelly  [1961]. 

RCTCTP  was  introduced  by  Talbot  [1982],  In  his  work,  activity  execution 
cannot  be  interrupted  by  another  activity  (nonpreemptive  case),  and  each  activity  has  a 
set  of  possible  durations.  Different  sets  of  resources  are  required  for  different  durations. 
A duration-resource  combination  is  called  a mode.  Further,  he  distinguished  resources 
into  three  categories  by  using  the  terminology  introduced  by  Weglarz  [1979]  and 
Slowinski  [1980], [1981].  If  the  amount  of  resource  available  during  each  period  is 
limited  and  independent  of  previous  consumption  of  the  resource,  then  the  resource  is 
renewable.  An  example  of  a renewable  resource  is  a machine.  If  the  total  amount  of 
resource  available  over  the  life  or  part  of  the  life  of  the  project  is  limited,  then  the 
resource  is  nonrenewable.  An  example  of  a nonrenewable  resource  is  money.  If  a 
resource  is  constrained  by  period  and  by  the  total  availability  during  the  project  life,  then 
the  resource  is  doubly  constrained.  An  example  is  a monetary  budget.  A budget  for  a 
project  is  often  limited  and  the  budget  per  period  can  be  restricted.  Since  each  mode  has 
a unique  set  of  resource  requirements,  it  incurs  a unique  cost.  Talbot  formulated  this 
problem  as  an  integer  program  and  presented  a two-stage  solution  procedure.  The 
objective  of  this  problem  is  minimization  of  the  project  cost  or  minimization  of  the 
project  completion  time. 
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Our  problem  is  different  from  the  two  problems  mentioned  above,  and  part  of  the 
difference  comes  from  the  treatment  of  activities.  We  assume  that  activity  duration  is 
determined  by  the  resource  set  to  be  used  and  the  degree  of  crashing.  As  an  example, 
let  us  consider  a ground  evening  activity.  Suppose  that  evening  ground  requires 
bulldozer(s)  and  that  a big  and  two  small  bulldozers  are  available.  Assume  that  the  job 
can  be  done  in  ten  days  by  using  two  small  bulldozers  or  in  seven  days  by  using  a big 
bulldozer,  assuming  eight  hours  of  operation  per  day  in  both  cases.  Thus,  activity 
duration  is  a function  of  the  resource  set  to  be  used  (two  small  bulldozers  or  one  big 
bulldozer).  The  difference  in  the  durations  comes  from  the  variation  in  the  productivity 
of  the  resource  set.  As  one  may  imagine,  there  can  be  more  than  two  modes  for  the 
ground  evening  activity.  For  example,  in  the  case  of  one  big  bulldozer,  we  can  assign 
five,  six,  seven  or  eight  workers  to  support  the  activity.  A second  example  is  allocation 
of  three  or  four  dollies  (for  carrying  soil).  If  we  specify  all  possible  resource 
combinations  for  executing  an  activity,  there  will  be  many  modes,  but  it  may  not  be 
worthwhile  to  investigate  all  possibilities.  Therefore,  in  specifying  modes,  we  will 
consider  only  the  most  important  resources  such  as  bulldozer,  remicon  truck,  specialized 
labor  and  so  on.  These  resources  are  reusable  and  relatively  scarce,  and  it  is  not  easy 
to  acquire  additional  units.  Resources  such  as  dollies  or  unskilled  part-time  labor  are  not 
considered  in  the  construction  of  modes  because  the  additional  units  of  resources  can  be 
acquired  by  higher  direct  costs.  The  other  factor  that  determines  an  activity  duration  is 
the  degree  of  crashing  to  be  used.  The  first  type  of  crashing  is  assigning  more 
supporting  resources  to  the  activity.  In  the  case  of  the  ground  evening  activity,  the 
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activity  duration  with  one  big  bulldozer  can  be  shortened  by  assigning  more  supporting 
labor  or  equipment  such  as  dollies.  Assigning  these  resources  increases  activity  cost. 
The  second  type  of  crashing  is  using  overtime.  We  assumed  that  the  job  takes  ten  days 
with  two  small  bulldozers  based  on  eight  hours  of  operation  per  day.  This  implies  that 
the  activity  requires  160  small  bulldozer  hours  of  operation  (2  small  bulldozer  * 8 hours 
* 10  days).  If  we  add  a second  shift  (a  total  of  sixteen  hours  per  day),  then  the  same  job 
can  be  done  in  only  five  days.  Or  if  we  use  just  two  hours  of  overtime  per  day,  it  takes 
eight  days.  Using  overtime  or  the  second  or  third  shifts  may  increase  the  activity  cost 
since  overtime  labor  and/or  additional  shifts  usually  cost  more  than  normal  shift.  As  we 
can  imagine  many  combinations  of  the  two  types  of  crashing,  there  will  be  many  possible 
ways  of  crashing  an  activity  in  a mode.  Let  us  refer  to  the  lowest  cost  of  completing  one 
activity  with  a particular  mode  as  the  mode  normal  cost  and  the  corresponding  duration 
as  the  mode  normal  duration.  In  a similar  way,  let  us  refer  to  the  shortest  duration  of 
completing  an  activity  with  a given  mode  as  the  mode  crash  duration  and  the 
corresponding  activity  cost  as  the  mode  crash  cost.  It  is  obvious  that  mode  crash  cost 
is  equal  to  or  higher  than  the  mode  normal  cost  since  we  define  the  mode  normal  cost 
as  the  least  expensive  way  of  completing  an  activity  with  a particular  mode.  While  any 
duration  between  the  normal  and  crash  durations  can  be  realized  by  varying  the  degree 
of  overtime  or  allocating  additional  resources,  the  corresponding  cost  function  is  not 
necessarily  linear.  However,  for  simplicity  we  will  assume  that  the  cost  function 
between  the  normal  and  crash  durations  is  linear.  Since  multiple  modes  are  assumed  for 
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each  activity,  there  exist  multiple  normal  duration/costs  and  crash  duration/costs  for  each 
activity. 

The  differences  among  the  two  problems  and  ours  are  summarized  as: 

TCTP:  Activity  duration  and  cost  is  a function  of  monetary  spending  while  all  the 

resource  constraints  are  ignored. 

RCTCTP:  Activity  duration  and  cost  is  a function  of  the  mode  to  be  used  while 

resource  constraints  are  considered. 

OURS:  Activity  duration  and  cost  is  a function  of  the  mode  to  be  used  and  the 

degree  of  crashing  implemented  while  resource  constraints  are  considered. 
We  believe  that  our  treatment  of  activities  connects  the  Time/Cost  Tradeoff 
Problem  and  the  Resource  Constrained  Project  Scheduling  Problem  in  a realistic  way. 
Of  course,  if  we  find  all  the  possible  durations  between  each  mode’s  normal  and  crash 
durations  and  convert  these  durations  into  individual  modes,  then  our  problem  becomes 
RCTCTP.  In  this  case,  the  number  of  modes  per  activity  can  increase  tremendously. 
One  of  the  advantages  of  our  approach  is  that  the  number  of  modes  per  activity  is  smaller 
than  that  in  RCTCTP  after  the  conversion  mentioned  above.  A second  advantage  is  a 
more  natural  interpretation  of  the  activity  duration/cost  function  since  the  activity  cost 
and  the  duration  are  determined  by  the  mode  and  the  degree  of  crashing  used. 

To  our  knowledge,  there  is  no  work  for  the  nonpreemptive  case  in  the  literature 
that  considered  activity  duration  and  cost  as  a function  of  the  mode  and  the  degree  of 
crashing  to  be  used.  The  objective  of  our  problem  is  to  find  a feasible  schedule  that 
minimizes  the  total  project  cost.  A schedule  consists  of  specifying  each  activity’s  mode, 
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duration  and  start  time.  Feasibility  is  achieved  by  satisfying  the  precedence  relations  and 
resource  constraints  in  all  periods  during  the  project  life. 

In  the  next  chapter,  we  present  a survey  of  the  project  scheduling  problems.  This 
survey  focuses  on  the  resource  constrained  project  scheduling  problems,  time/cost 
tradeoff  problem  and  resource  constrained  time/cost  tradeoff  problems.  In  Chapter  3, 
we  present  the  formulation  of  our  problem.  Chapter  4 is  devoted  to  an  exact  solution 
procedure  including  the  explanation  of  some  newly-defined  concepts.  Our  exact  solution 
procedure  is  a branch  and  bound  algorithm,  where  a linear  programming  relaxation  is 
used  for  bounding.  In  Chapter  5,  an  efficient  heuristic  procedure  is  introduced.  Chapter 
6 includes  the  computational  results  for  both  the  exact  and  heuristic  procedures. 
Conclusions  and  a discussion  of  future  research  are  presented  in  Chapter  7. 


CHAPTER  2 
LITERATURE  REVIEW 

One  of  the  earliest  complete  surveys  of  project  scheduling  problems  was  done  by 
Davis  [1973],  and  the  most  recent  survey  was  presented  by  Icmeli  [1992],  Icmeli 
focused  on  literature  since  1973  and  summarized  how  the  basic  project  scheduling 
problems  (such  as  RCPSP,  TCTP  and  the  payment  scheduling  problem)  are  combined 
to  generate  a new  class  of  project  scheduling  problems.  The  author  concluded  that  "the 
future  of  the  project  scheduling  literature  appears  to  be  developing  in  the  direction  of 
combining  the  fundamental  problems  and  developing  efficient  exact  and  heuristic  methods 
for  the  resulting  problems."  [1992,  p.  7], 

Since  a thorough  and  recent  literature  review  for  the  project  scheduling  problems 
is  available,  we  will  concentrate  on  the  development  in  the  three  areas  that  are  closely 
related  to  the  problem  studied  in  this  dissertation.  The  three  areas  are  the  time/cost 
trade-off  problem  (TCTP),  the  resource  constrained  project  scheduling  problem  (RCPSP), 
and  the  resource  constrained  time/cost  trade-off  project  scheduling  problem  (RCTCTP). 

2.1  Time/Cost  Trade-Off  Problem 

The  typical  TCTP  assumes  that  an  activity  duration  can  be  shortened  from  a 
normal  duration  up  to  a crash  duration  through  higher  direct  cost  assuming  a linear  cost 
function.  While  precedence  relations  are  considered,  the  resource  requirements  for 
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executing  each  activity  are  not  explicitly  considered.  The  major  question  that  can  be 
answered  by  TCTP  is  "what  activities  (if  any)  should  be  expedited  to  achieve  the  desired 
project  duration  at  a minimal  additional  cost?"  as  described  by  Tufekci  [1982,  p.  109]. 

Fulkerson  [1961]  and  Kelly  [1961]  formulated  the  time/cost  trade-off  problem  as 
a linear  program.  It  should  be  noted  that  if  we  are  concerned  only  with  a specific  project 
completion  time,  then  the  question  can  be  tackled  by  linear  programming.  The  primary 
concern  of  TCTP  is  not  a single  project  duration,  but  it  is  all  the  feasible  project 
durations.  Both  Fulkerson  and  Kelly  tried  to  identify  the  activities  to  be  shortened  and 
the  minimum  project  cost  for  each  feasible  project  duration  by  parameterizing  the  project 
duration.  The  linear  program  was  dualized  and  converted  to  a flow  problem.  The 
algorithm  starts  with  the  maximum  project  duration  (which  is  the  length  of  the  longest 
path  from  the  start  node  to  the  end  node  while  each  activity  duration  is  assumed  to  be 
normal),  and  successive  reductions  on  the  project  duration  are  attempted  by  an  iterative 
flow  algorithm.  Kelly  showed  in  the  same  paper  [1961]  how  his  algorithm  could  handle 
activity  utility  functions  that  are  piecewise  linear,  nondecreasing  and  concave;  the 
objective  was  the  maximization  of  utility  among  all  the  feasible  project  schedules  of  the 
same  project  duration.  By  utilizing  the  characteristics  of  the  piecewise  linear, 
nondecreasing  and  concave  function,  Kelly  showed  how  to  reduce  the  nonlinear  problem 
into  a linear  problem  by  breaking  activities  into  subactivities.  Clark  [1961]  used  an 
incremental  allocation  scheme  to  find  the  minimum  cost  of  successive  reductions  in  the 
project  duration  where  the  cost  function  was  assumed  to  be  positive  and  decreasing. 
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The  research  on  TCTP  has  continued  in  two  ways:  improving  the  efficiency  of 
the  algorithm  and  generalizing  the  activity  cost  function. 

Three  excellent  papers  that  describe  efficient  algorithms  are  worthwhile  to 
introduce.  Phillips  and  Dessouky  [1977],  [1979]  introduced  a minimal  cut  concept.  The 
minimal  cut  was  located  directly  by  a cut-seeking  procedure.  They  considered  only  the 
critical  activities  in  the  flow  network.  Their  algorithm  requires  the  existence  of  a 
feasible  flow  in  the  network  before  the  algorithm  can  be  applied.  Through  their 
computational  results,  they  demonstrated  that  their  algorithm  outperformed  the  traditional 
out-of-kilter  algorithm.  A more  efficient  algorithm  was  introduced  by  Tufekci  [1982]. 
He  introduced  an  algorithm  that  finds  the  Phillips  and  Dessouky  network  cut  by  using 
a labeling  algorithm.  In  his  algorithm,  the  maximal  flow  found  at  the  k-th  iteration  is 
a feasible  flow  for  the  flow  network  at  the  (k-(-l)-th  iteration.  By  preserving  the  flow 
and  labels  on  the  flow  network,  the  required  number  of  augmentations  for  obtaining  a 
maximal  flow  and  a minimal  cut  at  the  succeeding  step  is  reduced,  and  the  computational 
effort  is  improved  over  the  algorithm  suggested  by  Phillips  and  Dessouky. 

Work  on  the  generalization  of  the  activity  cost  function  has  continued. 
Lamberson  and  Hocking  [1970]  and  Berman  [1964]  provided  efficient  algorithms  to  deal 
with  (nonlinear)  convex  activity  cost  functions.  Since  the  objective  is  the  minimization 
of  the  sum  of  activity  costs  and  activity  costs  are  convex,  any  locally  optimal  solution 
is  a global  optimal  solution.  Falk  and  Horowitz  [1972]  extended  a global  optimal 
solution  procedure  for  concave  or  piecewise  linear  and  convex  cost-time  functions. 
When  the  activity  cost  function  is  concave,  a linear  cost  curve  is  used  for 
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underestimating,  and  the  resulting  problem  is  a linear  programming  problem.  Whenever 
the  solution  of  the  linear  programming  problem  causes  an  underestimation  gap  between 
the  linear  cost  function  value  and  the  concave  cost  function  value  of  some  activity,  the 
feasible  region  of  the  activity  duration  is  divided  into  two  subsets,  and  a tighter 
underestimating  linear  cost  curve  for  each  subset  is  derived.  The  division  of  the  feasible 
region  and  the  derivation  of  the  tighter  underestimating  cost  curve  are  incorporated  in  a 
branch  and  bound  algorithm.  Robinson  [1975]  introduced  a dynamic  programming 
algorithm  that  can  handle  TCTP  with  arbitrary  but  nonincreasing  cost  functions.  The 
objective  is  the  minimization  of  the  project  duration.  Elmaghraby  [1993]  also  used  a 
dynamic  programming  technique  for  the  project  completion  time  minimization  problems 
that  have  a monotonically  nonincreasing  cost  function.  In  his  paper,  a node  reduction 
technique  was  used,  and  as  a result,  the  complexity  of  the  dynamic  programming 
approach  was  drastically  reduced  from  0(Uh)  to  0(UC),  where  h is  the  number  of 
common  activities  and  c is  the  minimum  number  of  reduced  nodes.  An  optimal 
procedure  for  the  discrete  time/cost  trade-off  problem  was  introduced  by 
Demeulemeester,  Herroelen  and  Elmaghraby  [1993].  This  procedure  uses  dynamic 
programming  and  branch  and  bound  techniques.  Nair,  Prasad  and  Aneja  [1993] 
presented  an  optimization  procedure  based  on  a labelling  scheme  for  the  problem  in  that 
the  activity  cost  function  is  nonincreasing,  piecewise  linear  and  convex. 

Siemens  [1971]  presented  a heuristic  procedure,  Siemens  Approximation  Method 
(SAM),  for  the  project  scheduling  problem.  The  activity  cost  function  was  convex.  The 
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author  demonstrated  the  algorithm’s  efficiency  by  comparing  the  solutions  from  SAM 
with  solutions  found  by  a linear  programming  method. 

2,2  Resource  Constrained  Project  Scheduling  Problem 

In  a survey  paper,  Davis  [1973]  classified  resource  allocation  problems  into  three 
categories:  Time/Cost  Trade-off  Problems,  Leveling  Resource  Demand  and  Project 
Scheduling  under  Fixed  Resource  Limits.  He  author  characterized  his  categorization  in 
terms  of  resource  allocation.  In  TCTP,  the  duration  of  some  or  all  of  the  activities  in 
a project  can  be  shortened  by  the  allocation  of  more  resources.  In  this  problem, 
allocation  of  more  resources  is  possible  at  the  expense  of  higher  activity  direct  cost,  and 
the  amount  of  available  resources  is  assumed  unlimited.  The  Leveling  Resource  Demand 
problem  occurs  when  there  are  enough  resources  to  concurrently  schedule  all  jobs  sharing 
the  same  resource  types,  but  it  is  desired  to  keep  the  resource  utilization  ratio  at  a certain 
level.  The  objective  of  this  problem  is  to  smooth,  as  much  as  possible,  the  resource 
profile(s)  over  time  during  the  project  life.  The  Project  Scheduling  Problems  under 
Fixed  Resource  Limits  occurs  when  there  are  only  fixed  amounts  of  resources  available 
in  each  period  in  the  project  life.  Davis  remarked  that  this  problem  is  representative  of 
a class  of  combinatorial  problems,  which  are  characterized  by  a factorial  growth  in 
computational  effort  as  problem  size  increases.  Davis  also  discussed  the  relationship 
between  this  problem  and  Assembly  Line  Balancing.  He  noted  that  ALB  "is  concerned 
with  repetitive  operations  and  large  numbers  of  identical  products,  while  the  project 
scheduling  problem  is  usually  a one-time,  one-of-a-kind  operation."  [1973,  p.  299]. 
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Davis  also  compared  the  project  scheduling  problem  to  the  job  shop  problem  and  pointed 
out  that  the  job  shop  problem  has  "the  continuous  nature  of  work  input  and  flow  in  the 
shop.  . . . Thus,  from  a managerial  point  of  view,  there  is  no  end  to  the  scheduling 
effort  and  this  had  led  many  job  shop  researchers  to  concentrate  on  the  development  of 
sequencing  rules  which  are  statistically  superior  over  the  long  run,  i.e.,  in  terms  of 
average  resource  utilization,  average  delay,  etc."  [1973,  p.  300].  He  concluded  that  in 
spite  of  these  differences,  there  exist  conceptual  similarities  between  the  project 
scheduling  problem  and  the  assembly  line  balancing  problem  or  between  the  project 
scheduling  problem  and  the  job  shop  problem  to  stimulate  interest  in  the  reciprocal 
applicability  of  solution  procedures.  He  classified  project  scheduling  problems  in  terms 
of  the  number  of  resource  types:  (1)  only  one  resource  type  which  is  common  to  all 
jobs,  (2)  more  than  one  resource  type  per  period,  but  only  one  type  per  job  with 
availabilities  and  requirements  limited  to  one  resource  unit  (traditional  job  shop  problem 
assumptions)  and  (3)  more  than  one  resource  type  per  job  (multiple  resource  case). 
Davis  concluded  that  there  was  conflicting  evidence  as  to  which  sequencing  heuristic  was 
dominating  and  that  there  was  no  optimal  procedure  that  had  been  demonstrated  as 
computationally  feasible  for  the  sorts  of  large,  complex  projects  that  occurred  in  practice. 

Since  RCPSP  usually  is  of  the  multiple  resource  variety,  we  will  concentrate  on 
the  multiple  resource  constrained  project  scheduling  problem.  Multiple  resource  case  can 
be  characterized  by  constant  resource  levels,  variable  resource  levels,  single  mode  case 
and  multiple  mode  case. 
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2.2.1  RCPSP  Formulation 

Manne  [1960]  formulated  the  job  shop  scheduling  problem  in  that  an  integer 
variable  indicated  a job  start  time.  Wagner  [1959]  formulated  the  machine  scheduling 
problem  as  an  integer  linear  program.  0-1  variables  were  used  to  indicate  whether  a job 
was  assigned  to  a specific  order  position  on  a specific  machine.  Neither  formulation 
included  multiple  resource  constraints.  Bowman  [1959]  used  0-1  variables  to  indicate 
whether  a job  was  being  processed  in  each  period  during  the  scheduling  horizon.  One 
of  the  earliest  successful  formulations  was  made  by  Pritsker,  Walters  and  Wolfe  [1969]. 
In  their  0-1  integer  programming  formulation  for  the  multiproject  scheduling  and  the  job 
shop  scheduling  problems  with  limited  resources,  0-1  variables  indicated  whether  a job 
was  completed  in  a period.  They  suggested  that  "for  a given  type  of  scheduling 
environment,  alternative  formulations  often  may  be  devised.  An  efficient  formulation, 
however,  will  depend  upon  a judicious  choice  of  definition  for  the  variables."  [1969,  p. 
94],  A drawback  of  the  0-1  formulation  of  Pritsker  et  al.  is  that  the  number  of  variables 
increases  very  rapidly  as  the  problem  size  increases.  Patterson  and  Huber  [1974] 
formulated  RCPSP  as  a 0-1  integer  linear  program.  Their  formulation  is  similar  to  that 
of  Pritsker  et  al.  [1969]  in  that  the  0-1  variables  are  used  to  indicate  whether  a job  is 
completed  in  a given  period.  Patterson  and  Roth  [1976]  briefly  reviewed  the  literature 
with  respect  to  the  definition  of  problem  variables  and  concluded  that  "a  judicious  choice 
of  variable  definition  will  lead  to  a formulation  of  the  sequencing  problem  which  is 
readily  solvable,  whereas  an  efficient  formulation  may  result  in  the  fewest  number  of 
zero-one  variables  or  the  fewest  number  of  constraints  being  required  to  represent  a given 
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sequencing  problem.  . . . Thus,  variable  definition  is  as  important  as  it  relates  to  the 
resultant  problem  structure  as  it  is  to  the  size  of  the  zero-one  problem  (measured  by  the 
number  of  zero-one  variables  and/or  constraints  required)."  [1976,  p.  450].  In  their 
formulation,  a 0-1  variable  indicates  whether  a job  is  started  in  a specific  period  or  not. 
Their  compact  formulation  allowed  them  to  use  an  implicit  enumeration  algorithm  in  a 
special  way.  Talbot  and  Patterson  [1978]  presented  a conceptual  formulation  for  RCPSP. 
Resource  constraints  were  not  explicitly  included,  so  it  was  not  possible  to  solve  the 
problem  with  a standard  integer  programming  technique. 

2.2,2  RCPSP  with  Constant  Resource  Level  and  Single  Mode 

In  this  section,  we  will  consider  the  multiple  resource  constrained  project 
scheduling  problem.  The  objective  of  this  problem  is  the  minimization  of  the  project 
completion  time,  and  the  resource  availability  per  period  is  assumed  constant.  General 
assumptions  for  this  problem  include  a fixed  activity  duration  (a  single  mode), 
nonpreemptive  activities,  resource  constraints  and  precedence  constraints.  If  any  of  the 
assumptions  are  altered  by  any  work  listed  in  this  section,  the  modification  will  be  noted. 

Since  the  1960s,  many  techniques  such  as  integer  linear  programming,  dynamic 
programming  and  implicit  enumeration,  have  been  used  to  solve  RCPSP  with  constant 
resource  levels.  As  Davis  discussed  in  his  survey  [1973],  solving  this  problem  directly 
by  using  integer  linear  programming  is  computationally  impractical.  Davis  and  Heidorn 
[1971]  proposed  an  optimal  solution  procedure  using  bounded  enumeration  that  was 
originally  devised  by  Gutjahr  and  Nemhauser  [1964]  for  the  assembly  line  balancing 
problem.  The  advantages  of  the  procedure  are  that  the  resource  requirements  can  vary 
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over  the  job  duration  and  that  various  assumptions  with  respect  to  the  job  continuity  are 
allowed,  with  no  additional  computational  effort.  The  algorithm  generates  feasible 
(activity)  subsets,  sets  of  tasks  that  can  be  processed  simultaneously.  The  major 
drawback  of  their  algorithm  is  that  the  number  of  feasible  subsets  tends  to  increase 
sharply  as  the  problem  size  increases.  Johnson  [1967]  and  Schrage  [1970]  used  branch 
and  bound  algorithms  to  solve  the  problem.  The  problem  considered  by  both  was  a 
single  resource  constrained  problem.  The  first  computationally  competitive  solution 
procedure  was  introduced  by  Stinson,  Davis  and  Khumawala  [1978].  They  formulated 
RCPSP  as  a conceptual  linear  program  and  presented  a branch  and  bound  based 
procedure,  which  resembles  that  of  Johnson  [1967].  But  Stinson  et  al.  derived  a 
stronger  lower  bound,  which  considers  precedence  and  resource  constraints,  and 
employed  a node  selection  heuristic.  The  stronger  bound  improved  computational 
efficiency,  and  the  node  selection  heuristic  strengthened  both  the  dominance  and  the 
lower  bound  pruning  of  the  branch  and  bound  tree.  The  dominance  rule  is  partially 
inherited  from  the  works  of  Johnson  [1967]  and  Schrage  [1970].  The  lower  bound 
pruning  is  based  on  the  critical  path  for  the  unscheduled  activities  while  ignoring 
resource  constraints.  For  the  efficient  use  of  the  dominance  rules,  which  rely  upon 
comparisons  of  nodes,  they  used  skiptracking.  In  contrast  to  backtracking,  skiptracking 
forces  the  branch  and  bound  tree  to  expand  in  breadth  and  to  move  downward  in  a 
somewhat  uniform  fashion.  When  a problem  is  very  large,  a node  selection  rule  based 
on  a single  attribute  has  limited  effectiveness.  In  1974,  Ashour,  More  and  Chiu  [1974] 
pointed  out  this  difficulties  and  proposed  a node  selection  rule  based  on  a series  of  partial 
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schedule  attributes,  which  was  termed  a decision  vector  by  Ashour  et  al.  Stinson  et  al. 
adopted  the  idea  of  Ashour  et  al.  and  used  a decision  vector  consisting  of  eight  partial 
schedule  attributes.  Computational  results  with  up  to  sixty  activities  and  six  resource 
classes  were  presented. 

Christofides  et  al.  [1987]  proposed  a depth-first  search  branch  and  bound 
technique.  Four  lower  bounds  were  examined.  The  first  bound  was  based  on  the  longest 
path.  The  second  bound  was  based  on  linear  programming  relaxation  with  multiple  types 
of  cutting  planes.  The  third  bound  was  based  on  a Lagrangean  relaxation  of  the 
formulation.  The  idea  of  using  Lagrangean  relaxation  in  order  to  obtain  bounds  had 
already  been  tried  by  Fisher  [1973].  Relaxation  of  the  resource  constraints  results  in  a 
relaxed  problem,  which  is  easy  to  solve.  The  last  bound  was  based  on  disjunctive  arcs. 
The  idea  of  using  disjunctive  arcs  was  proposed  by  Balas  for  the  job  shop  problem 
[1969]  and  the  project  scheduling  problem  [ 1970]  and  was  used  by  Gorenstein  [1972]  for 
the  project  (job)  sequencing  problem  with  resource  constraints.  Christofides  et  al.  [ 1987] 
reported  computational  results  on  the  performance  of  each  bound  on  problems  that  have 
up  to  twenty-five  activities  and  three  types  of  resources.  The  first  bound  was  useless 
when  resource  constraints  were  tight.  The  bounds  generated  by  using  linear  relaxation 
with  valid  cutting  planes  were  generally  better  than  those  generated  by  using  Lagrangean 
relaxation. 

Demeulemeester  and  Herroelen  [1992]  proposed  a branch  and  bound  procedure 
that  is  computationally  the  fastest  for  RCPSP.  Their  procedure  adopted  a depth-first 
search.  In  the  branch  and  bound  tree,  nodes  represent  resource  and  precedence  feasible 
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partial  schedules.  Branches  from  a parent  node  correspond  to  exhaustive  and  minimal 
combinations  of  activities,  where  the  delay  of  these  activities  resolves  resource  conflicts 
(resource  constraint  violation  at  some  period)  at  each  parent  node.  Precedence  and 
resource-based  bounds  are  used.  When  the  depth-first  procedure  does  not  allow  the  tree 
pruning  by  schedule  dominance  that  was  used  by  Stinson  et  al.  [1978],  two  dominance 
rules  are  used  to  prune  the  search  tree:  a left-shift  dominance  rule  and  a dominance  rule 
based  on  the  concept  of  a cutset.  The  left-shift  dominance  rule  shares  a conceptual 
similarity  with  the  one  used  by  Schrage  [1970],  Herroelen  [1972]  and  Stinson  et  al. 
[1978].  At  every  time  instance,  a cutset  is  defined  as  the  set  of  all  unscheduled  activities 
for  which  all  predecessor  activities  have  already  been  scheduled.  By  comparing  cutsets, 
any  partial  schedule  that  is  verified  not  to  be  optimal  is  fathomed.  The  computational 
results  with  Patterson’s  1 10  test  problems  showed  that  their  procedure  outperformed  the 
procedure  of  Stinson  et  al. 

2.2.3  RCPSP  with  Variable  Resource  Level  and  Single  Mode 

In  this  section,  we  will  review  the  development  of  RCPSP  for  the  case  in  that  the 
amount  of  available  resources  per  period  is  not  constant  and  each  activity  has  a single 
mode  with  fixed  duration. 

Patterson  and  Huber  [1974]  formulated  the  problem  as  a 0-1  integer  linear 
program  and  used  a minimum  and  a maximum  bounding  techniques,  which  are  called  the 
horizon-varying  approach.  The  bounding  technique  was  originally  developed  by  Gutjahr 
and  Nemhauser  [1964],  The  minimum  bounding  technique  starts  with  the  smallest 
possible  project  duration  and  checks  whether  there  exists  a feasible  solution  that  satisfies 
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this  duration.  If  a feasible  solution  exists,  then  the  procedure  terminates.  If  not,  the 
project  duration  is  increased  by  one  period,  and  the  process  is  repeated  until  the  first 
feasible  solution  (which  is  also  optimal)  is  found.  The  maximum  bounding  technique  is 
analogous  to  the  minimum  bounding  technique,  but  it  starts  with  a sufficiently  long 
project  duration  and  terminates  when  there  is  no  feasible  solution  to  satisfy  this  project 
duration.  An  attempt  to  combine  two  bounding  algorithms  and  perform  a binary  search 
was  tried  but  proved  fruitless.  Through  computational  experiments,  the  authors 
concluded  that  their  algorithm  is  more  effective  than  0-1  programming  approach  without 
a bounding  technique. 

Patterson  and  Roth  [1976]  formulated  the  problem  differently.  In  their 
formulation,  a 0-1  variable  indicates  whether  a job  is  started  in  a certain  period.  The 
problem  was  solved  by  an  implicit  enumeration  procedure  (a  modified  Balasian 
algorithm)  by  exploiting  the  special  structure  of  their  formulation.  The  presented 
computational  results  clearly  showed  that  the  implicit  enumeration  procedure  performed 
better  than  the  horizon-varying  approach. 

Talbot  and  Patterson  [1978]  presented  another  solution  procedure  based  on  the 
implicit  enumeration  algorithm  introduced  by  Balas  [1965].  Their  proposed  algorithm 
differs  from  the  general  Balasian  implicit  enumeration  in  the  sense  that  nonnegative 
integer,  instead  of  0-1  binary,  variables  are  used  and  that  variables  are  augmented 
sequentially  by  a job  number  rather  than  feasibility  tests.  A network  cut  concept  is 
employed  to  accelerate  the  enumeration.  A network  cut  is  an  integer  time  period  that  is 
used  to  evaluate  partial  feasible  solutions  as  good  or  bad.  The  evaluation  of  a partial 
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feasible  solution  enabled  the  algorithm  to  fathom  inferior  partial  solutions  before 
proceeding  further.  The  computational  experiments  revealed  that  their  algorithm  with 
the  network  cut  outperformed  the  algorithm  proposed  by  Patterson  and  Roth  [ 1 97 6 J . 

2.2.4  RCPSP  with  Multiple  Modes 

The  RCPSP  with  multiple  modes  was  introduced  by  Talbot  [1982]  and  followed 
by  Patterson,  Slowinski,  Talbot  and  Weglarz  [1989];  the  computational  results  of  the 
latter  were  reported  by  Patterson,  Slowinski,  Talbot  and  Weglarz  [1990],  Both 
algorithms  can  handle  scheduling  problems  with  a variety  of  objective  functions  and  with 
constant  or  variable  resource  levels.  The  algorithm  presented  by  Patterson  et  al.  [1989] 
differs  from  Talbot’s  [1982]  in  the  usage  of  the  precedence  tree  that  guides  the  search. 
The  tree  of  rank -ordered  activity  numbers  identifies  all  precedence-feasible  sequences  of 
activity  numbers  in  the  solution  space.  Recently,  Sprecher  [1994]  applied  five  bounding 
rules  to  the  algorithm  introduced  by  Patterson  et  al.  [1989]  and  verified  the  efficiency 
of  those  bounding  rules  by  presenting  favorable  computational  result. 

2.2.5  Heuristics  for  RCPSP 

The  first  systematic  effort  to  evaluate  the  efficiency  of  existing  heuristics  for 
RCPSP  with  a single  mode  was  done  by  Davis  and  Patterson  [1975],  Based  on  the 
previous  research,  they  selected  the  most  effective  eight  heuristics:  (1)  minimum  job 
slack,  (2)  resource  scheduling  method,  (3)  minimum  late  finish  time,  (4)  greatest 
resource  demand,  (5)  greatest  resource  utilization,  (6)  shortest  imminent  operation,  (7) 
most  jobs  possible,  and  (8)  select  jobs  randomly.  None  of  the  heuristic  rules 
outperformed  the  others  consistently,  but  the  minimum  job  slack  rule  most  often 
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produced  the  smallest  deviation  from  the  optimal  solution.  Boctor  [1990]  compared  13 
heuristic  procedures  for  a single  mode  case  and  concluded  that  the  minimum  slack  rule, 
the  minimum  latest  finish  time  rule  and  the  resource  scheduling  method  were  the  most 
efficient  heuristics.  He  also  compared  the  efficiency  of  the  multiheuristic  procedures. 
The  term  multiheuristic  procedure  refers  to  the  usage  of  more  than  one  heuristic  and 
taking  the  best  solution.  Boctor  [1993]  compared  21  heuristics  for  the  multiple  mode 
case  with  240  test  problems  on  various  conditions  and  concluded  that  more  effort  is 
needed  to  design  more  efficient  heuristics.  Drexl  and  Gruenewald  [1993]  introduced  a 
stochastic  heuristic  method  that  is  based  on  a weighted  random  selection  technique,  which 
was  successfully  used  by  Drexl  [1991].  They  also  presented  computational  results  that 
favor  the  stochastic  method  over  the  deterministic  method. 

2.2.6  Variation  of  RCPSP 

One  of  the  most  widely  studied  variations  of  RCPSP  is  the  net  present  value 
maximization  problem.  The  net  present  value  maximization  problem  without 
consideration  of  resource  constraints  was  introduced  by  Russell  [1970].  Grinold  [1972] 
successfully  transformed  the  problem  into  a linear  programming  problem.  RCPSP  with 
net  present  value  maximization  problem  was  introduced  by  Doersch  and  Patterson  [1977], 
and  subsequent  work  was  done  by  Patterson,  Slowinski,  Talbot  and  Weglarz  [1989], 
Patterson,  Slowinski,  Talbot  and  Weglarz  [1990]  and  Yang,  Talbot  and  Patterson  [1992]. 
While  earlier  work  in  RCPSP  with  net  present  value  maximization  problems  used  a 
partial  feasible  schedule  concept,  Icmeli  [1992]  used  a branch  and  bound  algorithm, 
where  the  bounding  problem  was  obtained  by  dropping  all  the  resource  constraints. 
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Mason  and  Moodie  [1971]  used  a branch  and  bound  algorithm  to  solve  a project 
scheduling  problem  that  combines  the  leveling  and  single  resource  constrained  problems. 
The  objective  was  the  minimization  of  the  total  project  cost,  which  was  the  sum  of  the 
resource  level  change  cost  and  the  project  delay  cost.  Smith-Daniels  et  al.  [1987]  dealt 
with  project  scheduling  with  materials  ordering  problems  and  used  a decomposition 
procedure  to  solve  it.  Weiss  [1988]  considered  the  scheduling  problem  for  the 
simultaneous  management  of  multiple  resource  constrained  project  networks  in  parallel. 
The  problem  is  formulated  as  an  integer  programming  problem,  which  is  similar  to  the 
just-in-time  materials  requirement  planning  schedule  model.  First,  the  problem  is  relaxed 
to  a linear  programming  problem,  and  then  the  Dantzig-Wolfe  Decomposition  is  used  to 
obtain  the  optimal  solution  to  the  linear  programming  problem.  The  LP  relaxation 
generates  a set  of  feasible  schedules  to  the  original  problem.  Finally,  feasible  schedules 
are  combined  to  form  a solution  to  the  original  problem. 

2.3  Resource  Constrained  Time/Cost  Trade-Off  Problem 

The  marriage  of  TCTP  and  RCPSP  forms  a new  project  planning  problem, 
RCTCTP.  For  the  preemptive  case,  Slowinski  [1980],  [1981]  and  Weglarz  [1981] 
formulated  and  solved  the  RCTCTP.  For  the  nonpreemptive  case,  Talbot  [1982] 
considered  the  resource  constrained  project  scheduling  problem.  Job  performance  in  this 
problem  is  a function  of  discrete  resource  consumption.  His  research  considered  all  three 
resource  types— renewable,  nonrenewable  and  doubly  constrained.  Since  each  activity 
duration  is  a function  of  resources  used,  and  since  the  usage  of  nonrenewable  resources 
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can  be  converted  to  a monetary  value,  several  different  objectives,  such  as  minimization 
of  project  completion  time  and  project  cost  minimization,  can  be  studied.  The  problem 
is  formulated  as  a 0-1  integer  linear  program.  The  presented  solution  procedures  are 
almost  identical  for  the  different  objectives  that  can  be  used.  A two  stage  algorithm  is 
proposed.  In  stage  one,  a network  is  relabeled  using  one  of  the  eight  scheduling  rules. 
The  label  on  each  activity  specifies  the  ordering  for  the  schedule  in  stage  two.  In  stage 
two,  each  activity  is  assigned  to  the  earliest  period  without  violating  resource  and 
precedence  constraints  or  without  exceeding  bounds  such  as  latest  finish  time.  If  the 
assignment  is  successful,  the  procedure  moves  to  the  next  activity.  If  not,  backtracking 
occurs,  and  the  most  recently  assigned  activity  is  rescheduled  to  a later  time  in  order  to 
make  enough  resources  available  for  the  activity  in  consideration.  This  process  is  applied 
to  all  the  activities  and  all  the  modes  until  an  improved  solution  is  found  or  until 
optimality  is  verified.  Whenever  an  improved  solution  is  found,  bounds  are  tightened, 
and  the  whole  process  is  repeated.  Since  the  suggested  methodology  finds  improved 
feasible  solutions  successively,  a good  feasible  solution  can  be  obtained  before  the 
completion  of  total  enumeration.  The  solution  procedure  is  used  as  a good  heuristic 
procedure  when  excessive  computational  effort  is  required  to  solve  the  problem 
optimally.  Talbot’s  report  on  computational  results  showed  that  the  algorithm  provides 
optimal  solutions  to  small  size  problems  and  good  heuristic  solutions  to  large  problems. 
Patterson,  Slowinski,  Talbot  and  Weglarz  [1989J  provided  new  life  for  this  problem.  As 
mentioned  in  Section  2.2.4,  the  algorithm  proposed  by  Patterson  et  al.  can  handle  general 
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objectives.  As  a result,  the  algorithm  can  solve  RCTCTP.  By  using  the  precedence  tree, 
the  algorithm  presented  by  Patterson  et  al.  outperforms  Talbot’s. 


CHAPTER  3 
THE  PROBLEM 


It  is  assumed  that  the  project  consists  of  I>0  nondummy  activities  and  has  a 
known  and  positive  due  date,  DUE_DATE>0.  The  notation  follows  the  Activity-on- 
Node  format,  where  nodes  represent  activities.  Labeling  of  nodes  (activities)  also 
follows  the  convention  in  that  a succeeding  node  gets  a higher  number  than  all  of  its 
predecessor  nodes.  Two  dummy  nodes  are  introduced:  node  0 to  denote  the  start  of  the 
project  and  node  (1  + 1)  to  denote  the  completion  of  the  project.  Both  dummy  activities 
have  zero  fixed  duration  and  incur  no  cost.  All  of  the  nondummy  activities  require  some 
resources  and  some  positive  cost  and  time  for  processing. 

3.1  Assumptions  and  Objective  of  the  Problem 

The  followings  are  the  assumptions  of  the  problem. 

(i)  Precedence  Constraints:  Each  activity  can  start  only  when  all  of  its 
technologically  preceding  activities  are  completed. 

(ii)  Resource  Requirement:  Each  activity  requires  some  set  of  resources  while  it  is 
processed.  The  types  and  amounts  of  resources  for  each  mode  are  known  and 
fixed  for  the  execution  of  each  activity.  Further,  the  types  and  numbers  of 
resources  in  a mode  do  not  change  by  the  degree  of  crashing  used.  The  total 
available  resource  levels  are  known  and  constant  during  the  project  life.  Once 
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an  activity  starts  with  a mode,  then  the  activity  takes  a fixed  amounts  of  resources 
during  the  whole  activity  duration,  where  the  amounts  and  the  types  of  the 
resources  are  predetermined  by  the  mode  chosen.  All  resources  are  reusable. 

(iii)  Mode  Restriction:  For  each  activity,  only  one  mode  is  chosen  from  the  available 
mode  set.  In  other  words,  if  activity  i is  started  using  mode  m,  then  activity  i 
must  be  processed  using  mode  m for  its  duration. 

(iv)  Variable  Duration:  Duration  of  each  activity  is  a function  of  mode  to  be  used  and 
the  degree  of  crashing  implemented. 

(v)  No  Preemption:  Once  an  activity  starts,  then  the  activity  cannot  be  interrupted. 
It  must  continue  to  be  processed  (during  the  operation  time)  until  its  completion. 

(vi)  Due  Date:  The  project  has  a known  due  date.  If  the  project  is  not  completed  by 
its  due  date,  then  a known  and  constant  positive  penalty  cost  is  incurred  in  each 
period  after  the  due  date  until  the  project  is  completed. 

(vii)  Integrality  of  Scheduling:  Each  activity  must  start  at  the  beginning  of  some 
period,  not  in  the  middle  of  the  period,  and  it  must  also  finish  its  operation  at  the 
end  of  some  period. 

(viii)  Integrality  of  Normal  and  Crash  Durations:  The  normal  and  the  crash  durations 
of  each  activity  are  assumed  to  be  positive  integers. 

The  objective  of  the  problem  is  to  schedule  all  of  the  activities  in  the  project  in 
a way  that  minimizes  project  cost  while  satisfying  the  precedence  and  resource 
constraints.  If  the  project  is  completed  by  the  due  date,  the  project  cost  is  the  sum  of 
all  activity  execution  costs,  otherwise  the  project  cost  is  increased  by  the  project  penalty 
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cost.  Scheduling  an  activity  means  specifying  the  activity  start  time,  duration  and  mode 
to  be  used.  Satisfying  the  precedence  constraints  means  that  the  start  time  of  an  activity 
should  be  no  earlier  than  the  completion  time  of  all  the  preceding  activities.  Satisfying 
the  resource  constraints  means  that  the  sum  of  required  resources  at  any  time  during  the 
execution  of  the  project  does  not  exceed  the  level  of  available  resources. 


3.2  Mathematical  Formulation  of  the  Problem 


The  mathematical  formulation  of  the  problem  is  presented. 

i 


[P]  min  £ £ xi(m)  * { NCi(m)  + 

i=l  meM(i) 

+ xI+1  * Penalty  * ( FI+1  ■ 

st  Fj  + Zj  < Fj, 

X!  /C  Xi(m)ri(m)k  * ’ 

ieS,  meM(i) 

F0  = 0 

XT  ^(m)  = 1 ’ 

meM(i) 

Xi(m)  ^ {0,1}, 

^i(m)y  i(m)c  — ^i(m)^i  — i(m)n’ 

Fj  integer. 


^(m)  yi(m)n  ^i(m)  } 

DUE_DATE  ) 

(1) 

(i,j)  e H 

(2) 

for  all  k and  t 

(3) 

(4) 

for  all  i,  m 

(5) 

for  all  i,  m 

(6) 

for  all  i,  m 

(7) 

i = 0,  ...,  1 + 1 

(8) 

Z; 


integer. 


i = 1,  ...,  I 
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where, 


vi+i 


0 if  FI+1  <;  DUE_DATE 

1 otherwise 


(9) 


Zq  = 0 


(10) 


w+i 


= 0 


i: 


H: 

M(i): 

i(m): 

^i(m) ' 

z{: 

y i(m)n‘ 
Y i(m)c" 
Ci(m)* 


NCj, 


i(m)- 


Fs: 


an  activity,  i = 0,  1,  ...,  1 + 1.  Activity  0 is  a dummy 
activity  denoting  the  start  of  the  project,  and  activity  (1+ 1) 
is  a dummy  activity  denoting  the  completion  of  the  project, 
set  of  ordered  pairs  of  activities  indicating  precedence 
relationship. 

the  set  of  modes  available  to  activity  i 
mode  of  activity  i,  m £ M(i) 

0-1  mode  variable;  1 if  activity  i assumes  mode  m,  0 
otherwise. 

duration  of  activity  i 

normal  mode  duration  of  activity  i using  mode  m 
crash  mode  duration  of  activity  i using  mode  m 
(nonnegative)  marginal  mode  crashing  cost  of  activity  i 
using  mode  m defined  between  normal  and  crash  mode 
durations  of  the  activity  i with  using  mode  m 
normal  mode  cost  of  activity  i using  mode  m 
start  time  of  activity  i,  i = 0,  1,  ...,  1 + 1. 
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ri(m)k:  the  units  of  resource  type  k required  to  process  activity  i 

using  mode  m 

K:  the  number  of  resource  types  required  for  the  project 

bk:  the  amount  of  resource  k available  per  period 

St:  set  of  activities  in  progress  in  time  interval  ]t-l,t]  such  that 

St  = {i  | Fi  < t-1  < Fi  + zj. 

DUEDATE:  project  due  date 

Penalty:  penalty  cost  per  period  incurred  by  delaying  the  project 

beyond  the  DUE  DATE. 

xI+1:  a 0-1  variable  indicating  whether  the  due  date  is  violated  or 

not. 

The  objective  function  (1)  is  the  minimization  of  the  total  project  cost,  which  is 
the  summation  of  activity  costs  and  the  total  penalty  cost  if  the  project  is  not  done  by  its 
due  date.  Individual  activity  cost  is  a function  of  the  mode  to  be  used  and  the  duration 
of  the  activity.  Let  us  consider  the  cost  of  activity  i.  The  set  of  available  modes  for 
activity  i is  M(i).  Among  M(i),  only  one  mode  is  selected.  Once  a mode  is  selected, 
an  appropriate  duration  is  also  selected.  Selecting  only  one  mode  among  M(i)  is 
guaranteed  by  constraint  (5).  Suppose  that  the  select  mode  for  activity  i is  m*.  Then, 
xi(mt)  takes  value  1,  and  all  the  other  mode  variables  for  activity  i will  be  set  to  0.  The 
activity  duration  variable,  Zj,  is  bounded  between  the  normal  mode  duration  (the  longest) 
and  the  crash  mode  duration  (the  shortest)  of  activity  i using  the  chosen  mode,  as  shown 
in  constraint  (7).  Since  all  of  the  other  xi(m)  are  set  to  0,  constraint  (7)  is  valid  for  any 
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Zj  bounded  by  yi(m,)n  and  yi(m»)c.  If  z x is  equal  to  the  normal  duration  of  mode  m*,  the 
activity  cost  is  the  normal  cost  of  mode  m*,  NCi(mt).  If  the  duration  is  shorter  than  the 
normal  duration,  this  implies  that  crashing  is  implemented,  and  the  activity  cost  increases 
by  the  marginal  mode  crashing  cost,  ci(mt),  for  each  period  the  duration  is  shortened.  The 
completion  time  of  activity  i is  the  start  time  of  activity  i,  Fi?  plus  its  duration,  zx. 
Constraint  (2)  ensures  that  the  technologically  succeeding  activity  starts  after  the 
completion  of  the  immediate  preceding  activity.  The  constraint  (3)  is  a conceptual 
formulation  of  resource  constraints.  For  each  activity  in  S„  only  one  mode  variable  is 
set  to  1 , and  the  amount  of  required  resources  is  specified  by  the  mode  chosen  for  each 
activity  in  St.  The  durations  of  the  two  dummy  activities  are  set  to  zero  in  constraint 
(10).  In  constraint  (9),  a 0-1  variable  is  set  to  zero  if  the  project  due  date  is  met  and  is 


set  to  1 otherwise. 


CHAPTER  4 

AN  EXACT  SOLUTION  PROCEDURE 


The  nonlinearity  of  the  problem  defined  in  chapter  3 comes  from  the  introduction 
of  the  mode  variables  and  also  from  the  consideration  of  resource  constraints.  If  each 
mode  has  a fixed  activity  duration,  then  explicit  representation  of  resource  constraints  is 
easier,  although  still  not  practical.  Since  an  activity  duration  can  be  any  number  of 
periods  between  each  mode’s  normal  and  crash  durations,  explicit  representation  of 
resource  constraints  becomes  more  difficult.  Our  proposed  solution  procedure  is  a 
branch  and  bound  technique,  which  starts  with  the  conceptual  formulation  and  solves 
linear  programs  as  bounding  problems  at  the  nodes  of  the  branch  and  bound  tree. 

Functions  that  are  necessary  for  the  construction  of  the  relaxed  problem  are 
defined  in  section  4.1,  and  the  branch  and  bound  procedure  is  presented  in  section  4.2. 
In  section  4.3,  some  refinements  to  our  solution  procedure  are  presented. 

4.1  Activity  Cost  Functions 

Three  types  of  activity  cost  functions  are  needed  for  our  exact  solution  procedure. 
These  are  mode  cost  function,  activity  cost  function  with  multiple  modes  and  convex 
activity  cost  function. 
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Definition  4. 1 Mode  Cost  Function 

Suppose  that  mode  m is  available  for  activity  i.  The  mode  cost  function  of 
activity  i using  mode  m,  Ci(m)(Zi),  is  defined  as  follows: 


NCj(m)  + Ci(m)  ^ y i(m)n 

+ CO 


Zi>  fOT  y«m)c  * Zi  * y i(m)n 
otherwise. 


where  0 < yi(m)c  < yi(m)n  and  ci(m)>  0.  (See  Figure  4.1.) 


Figure  4.1  Mode  Cost  Function 

The  mode  cost  corresponding  to  any  duration  shorter  than  the  mode  crash  duration 
or  longer  than  the  mode  normal  duration  is  set  to  infinity  since  activity  i cannot  be 
processed  in  such  a duration  using  mode  m.  Since  we  define  the  mode  normal  cost  as 
the  lowest  activity  cost  using  mode  m,  C^Z;)  > NCi(m)  holds  for  all  yi(m)c  < z{  < yi(m)n. 
We  assume  that  a constant  (and  positive)  crashing  cost,  ci(m)  > 0,  exists  for  any  duration 


reduction  from  the  mode  normal  duration. 
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Definition  4.2  Activity  Cost  Function  with  Multiple  Modes 

Suppose  activity  i has  a set  of  modes,  M(i)  = {1,2,..,  | M(i)  | },  available.  The 
Activity  Cost  Function  of  activity  i with  mode  set,  M(i),  Cj(Zj),  is  defined  as  follows: 

Ci(Zi>  Xi(l)>  — • » Xi(|M(i)|)  ) = X!  Xi(m)  * ^i(m)(Zi) 

meM(i) 

where  £ = 1 , xi(m)  6 {0,1}  for  all  m 6 M(i)  and  Ci(m)(Zi)  is  defined  in 

ieM(i) 

Definition  4.1.  (See  Figure  4.2.) 


Ci 


Figure  4.2  Activity  Cost  Function  with  Multiple  Modes 

Since  an  activity  is  processed  using  only  one  mode,  the  activity  cost  is  identical 
to  the  chosen  mode’s  cost.  The  variable,  xi(m),  is  used  to  indicate  which  mode  is 
selected.  If  mode  m is  not  selected,  then  xi(m)  is  set  to  zero,  and  its  mode  cost  is  not 
included  in  the  activity  cost.  If  mode  m*  is  chosen  for  the  execution  of  the  activity,  then 
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xi(m*)  is  set  to  1,  and  the  activity  cost  is  identical  to  m*’s  mode  cost.  In  this  case,  activity 
duration  must  be  bounded  by  mode  m*’s  normal  and  crash  durations  in  order  for  activity 
cost  to  be  finite. 

The  activity  cost  function  is  hard  to  incorporate  into  a linear  programming 
problem.  Thus,  we  will  define  a convex  activity  cost  function  that  underestimates  the 
activity  cost  function  for  all  activity  durations  and  modes  and  that  can  be  easily 
incorporated  into  a linear  program.  The  Convex  Activity  Cost  Function  is  defined  as 
follows: 


Definition  4.3  Convex  Activity  Cost  Function 

Suppose  activity  i has  a nonempty  set  of  modes  available,  M(i)  = 
{1,..,  | M(i)  | }.  Let  Sj  be  a nonempty  set  in  R2  such  that  Sj  = {(d^C1),  (d2,C2),  ..., 
(d2 1 M(i)  I ,C2 1 M(i)  I )},  where  (dm,Cm)  = (yi(m)n,Ci(m)(yi(m)n))  for  1 < m < | M(i)  | and 
(dm,Cm)  = (yi(m)0,Ci(m)(yi(m)c))  for  | M(i)  | + 1 < m < 2 | M(i)  | . Then,  the  Convex 
Activity  Cost  Function  of  activity  i,  CQ,  is  defined  as  follows: 


cq(Zi)  = 


min  | a*  C | a*  d = Zj  } 

for  min_dm  <,  z{  <,  max_dm 
+ °°  otherwise. 


where  a = [a„  a2,  ...  , a2 , M(i)  ( ]‘,  0 < am  < 1 for  m = 1,  ...,  2 | M(i)  | , a,  + a2 
+ •••  + a2|M(i)|  = 1,  d = [d1,  ...,  d2 1 M(i>  I ]l,  C = [C\  ...,  C2 ' M<i> ' ]‘,  rnin_dm  = min 
{dm  | 1 < m < 2 | M(i)  | } and  max_dm  = max  {dm  | 1 < m < 2 | M(i)  | }. 
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Sj  is  a collection  of  each  mode’s  (normal  duration,  the  corresponding  cost)  and 
(crash  duration,  the  corresponding  cost).  We  use  the  convex  hull  of  S*  for  the 
construction  of  CQ.  Notice  that  each  element  in  CONV  {Sj}  is  represented  as  (ati,  c/C) 
for  some  a,  where  a,  d and  C are  defined  in  Definition  4.3.  CQ  is  defined  as  finite 
between  min  {dm  | 1 < m < 2 | M(i)  | } and  max  {dm  | 1 < m < 2 | M(i)  | } and 
is  set  to  infinity  elsewhere.  Between  rnin_dm  and  rnax_dm,  CQ(Zi)  is  the  minimum  value 
of  a'C,  where  a‘d  = Z;.  The  convex  hull  of  Sj  and  the  Convex  Activity  Cost  Function, 
CQ,  are  illustrated  in  Figure  4.3  and  4.4. 


Figure  4.3  CONV  {SJ 

The  following  theorem  shows  that  CQ  is  a function  underestimating  Q for  all 


durations  and  mode  selections. 
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Convex  Activty  Cost 


Figure  4.4  Convex  Activity  Cost  Function 

Theorem  4.4  CQ  is  a Function  Underestimating  Q 

Suppose  that  activity  i has  a nonempty  set  of  modes,  M(i)  = {1,  2,  ..., 
| M(i)  | }.  Then,  CCj(Zj)  < C^,  xi(1),  ...,  xi(|  M(i)  ( })  for  all  z{  > 0 and  for  all  (xi(1),  ..., 
xi(  |mo)  | >),  where  Q and  CQ  are  defined  in  Definitions  4.2  and  4.3,  respectively. 

Proof:  Assume  that  all  the  assumptions  of  Theorem  4.4  hold.  Also  assume  that  Sj  is 
defined  in  Definition  4.3.  Suppose  that  (z*,xi(1)*,...,xi( , M(i) | >*)  is  an  arbitrary  duration 
and  mode  selection  such  that  xi(1)*  + ...  + xi(|M(i)()*  = 1 and  each  xi(m)*  E {0,1}  for 
all  m E M(i).  It  is  obvious  that  there  exists  x^*  such  that  xi(m#)*  = 1.  From  the 
definition  of  Q,  Cj(Zj*,xi(1)*,  ...,  xi( , M(i)  | >*)  is  equal  to  C^fa*).  z * must  be  one  of  the 
following  three  cases:  (i)  yi(m,)c  < z*  < yi(m.)n,  (ii)  z*  < yi(m*)c  and  (iii)  z*  > yl(m*)n. 

Suppose  that  (i)  holds.  Notice  that  Cj(mH.)  is  a linear  function  between  yi(m*)c  and 
yi(m*)n.  Since  a linear  function  is  also  a convex  function,  (Zj^C^CZi*))  can  be  expressed 
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as  a convex  combination  of  (yi(m*)n, Ci(m.)(yi(m*)ri))  and  (yi(m*)C,Ci(mt)(yi(m,)c)).  Since 
(y i(m*)n’ ^-'i(m*)(y i(m*)n))  and  (yi(m.)c,Ci(m*)(yi(m*>c))  G S„  there  exists  a*  such  that  (z^C^./z  *)) 
= (a'kFaT),  where  a,  d and  C are  defined  in  Definition  4.3.  Since 
C^*,  xi(1)*,  ...,  xj( , M(j) | )*) 

= Ci(m.)(z*) 

= aHC 

> min  {c/C  | c/d} 

= CQCzJ  (a) 

Suppose  that  (ii)  or  (iii)  holds.  Then,  from  the  definition  of  Ci(m),  C^^Zj*)  is 
infinite.  Therefore,  if  (ii)  or  (iii)  holds,  then  CCjfo*)  < Ci(m,)(zi:|t)  = C^Zj*,  xi(1)*,  ..., 

xi(|M(i)|)*) (b) 

Since  (Zj*,xi(1)*,.--,xi(  | M(i)  | }*)  is  arbitrarily  chosen,  (a)  and  (b)  prove  Theorem  4.4. 

Q.E.D. 

We  will  present  a procedure  for  constructing  the  CQ  in  Procedure  4.5  and 
provide  Theorem  4.9  for  the  validity  of  Procedure  4.5. 

Procedure  4.5  Procedure  for  Constructing  the  Convex  Activity  Cost  Function 

We  assume  that  we  know  the  mode  normal  duration/cost  and  mode  crash 
duration/cost  of  activity  i for  all  m G M(i). 

Step  1.  For  each  mode  in  M(i),  find  two  points:  (the  longest  duration  with  finite 

mode  cost,  corresponding  cost)  and  (the  shortest  duration  with  finite  cost. 
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corresponding  cost).  Denote  the  former  point  as  (dml,  Cml)  and  the  latter 
point  as  (d”2,  C”2). 

Step  2.  Collect  all  those  points  and  sort  them  in  decreasing  order  of  duration.  If 

ties  happen,  select  the  one  with  the  lowest  cost  and  discard  all  the  pairs 
that  have  the  same  duration. 

Step  3.  Relabel  the  sorted  points  as  (d\C'),  (d2,C2),  ...,  (dJ,CJ)  where  dJ  is  the  j-th 
longest  duration,  and  Cj  is  its  corresponding  cost. 

Step  4.  Let  k = 1,  where  k is  the  breakpoint  index  in  CC;. 

Let  (yjSCQ1)  = (d1^1). 

Step  5.  If  yk  = dJ,  we  have  all  the  breakpoints  in  CCj. 

Let  K(i)  = k,  where  K(i)  is  the  number  of  breakpoints  in  CCj. 
Go  to  Step  6. 

Otherwise,  find  j*  such  that 


£j*  _ 

= maximum 

dJ*  - yk  di  < yk 


Cj  - Ck 

j j k 

dJ  - ys 


If  j*  is  not  unique,  then  choose  the  one  with  the  largest  index 
(thus,  the  smallest  duration). 

Let  k = k + I and  (yik,CCik)  = (dj*,Cj*). 

Repeat  Step  5. 
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Step  6.  Define  CQ  as  follows: 


for  Zj  = yf 

for  y*  ± z{<  yf'1 2 3 4,  k=2,..,K(i) 


otherwise 


An  example  for  Procedure  4.5  is  provided  in  Example  4.6. 

Example  4.6  An  Example  for  Procedure  4.5 

Assume  that  activity  i has  two  modes.  The  normal  duration/cost  and  crash 
duration/cost  of  each  mode  are  summarized  below: 

Mode  (Normal  duration,  cost)  (Crash  duration, cost) 


Step  2.  (10,30),  (7,32),  (6,38),  (3,41) 

Step  3.  j (dj,Cj) 


1 (10,30) 

2 (7,32) 

3 (6,38) 

4 (3,41) 

Step  4.  k = 1.  (y,1,  CQ1)  = (10,30) 


1 (10,30) 

2 (7,32) 


(6,38) 

(3,41) 


Step  1.  For  mode  1,  (10,30),  (6,38)  and  for  mode  2,  (7,32),  (3,41). 
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Step  5. 


Step  5. 


Step  5. 


Step  6. 


Since  10  3,  find  j*. 

j 

(dj,Cj)  (0  - 

• Ci1)/(dj  - 

y.1) 

2 

(7,32) 

-0.6667 

< = 

max 

3 

(6,38) 

-2.0000 

4 

(3,41) 

-1.5714 

Thus,  j 

i*  = 2.  k = 

: 2,  (yi2,Ci2) 

= (7,32). 

Since  ' 

1 3,  find 

j*- 

j 

(dj,C)  (0  - 

■C2)/(dJ  - 

yi2) 

3 

(6,38) 

-6.0000 

4 

(3,41) 

-2.2500 

< = 

max 

Thus,  j 

j*  = 4.  k = 

: 3,  (yf.Cfl 

= (3,41). 

Since  2 

S = 3,  K(i) 

= k = 3. 

Go  to  Step 

6. 

Summary  for  the  breakpoints  i 

in  CQ 

k (yi\Cik) 


1 (10,30) 

2 (7,32) 

3 (3,41) 

CCj  is  set  as  follows: 


30.00 

for  zi  = 10 

46.00 

- 0.6667  * Zj 

for  7 £ Zj  < 10 

47.75 

- 2.2500  * z. 

for  3 <;  Z;  < 7 

+ oo 

otherwise 

order  to  prove  the  validity  of  Procedure  4.5,  the  following  propositions  are 


introduced. 
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Proposition  4.7  CC^y/)  = CQ1 

Suppose  that  activity  i has  a nonempty  set  of  modes,  M(i)  = {1,..,  | M(i)  | }. 
Then,  CQfyj1)  = CQ1,  where  CQ  is  defined  in  Definition  4.3  and  CQ1  and  y3  are  found 
by  using  Procedure  4.5. 

Proof:  Assume  that  activity  i has  a nonempty  set  of  modes,  M(i)  = {1,..,  | M(i)  | } and 
CQ  is  defined  in  Definition  4.3.  Further  assume  that  CQ1  and  y/  are  found  by  using 
Procedure  4.5.  Assume  that  dm,  1 < m < 2 | M(i)  | , d,  C and  a are  defined  in 
Definition  4.3.  Notice  that  Steps  1 to  3 in  Procedure  4.5  imply  that  y3  = max  {dm  | 

1 < m < 2 | M(i)  | }.  Suppose  that  c/d  is  equal  to  y/.  Then,  all  ara,  1 < m < 

2 | M(i)  | , such  that  dm  < y3  must  be  0 since  there  is  no  dm  greater  than  y/.  Observe 
the  followings: 

CQfo1) 

= min  {a'C  | a'd  = y/} 

from  the  definition  of  CQ 
= min  {Cm  | dm  = y/,  1 < m < 2 | M(i)  | } 

since  all  am,  1 < m < 2 | M(i)  | , such  that  dm  < y3  must  be  0. 

Notice  that  Steps  1 to  3 in  Procedure  4.5  imply  that  CQ1  is  min  {Cm  | dm  = y3, 
1 < m < 2 | M(i)  | }.  This  proves  Proposition  4.7.  Q.E.D. 

If  dm  = y/  for  all  m,  then  it  is  easy  to  see  that  Proposition  4.7  proves  the  validity 
of  Procedure  4.5.  Suppose  that  there  exists  at  least  one  dm  such  that  dm  < y/,  for  some 
m E { 1 , . . ,2  | M(i)  | }.  Then,  Procedure  4.5  finds  multiple  breakpoints.  In  other 
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words,  K(i)  > 2.  In  order  to  prove  the  validity  of  Procedure  4.5,  we  need  the  following 
proposition. 

Proposition  4.8  + bk(dm)  < Cm  for  all  m and  k. 

Assume  that  all  the  assumptions  of  Proposition  4.7  hold.  Also  assume  that  K(i) 
> 2,  where  K(i),  a*  and  bk  are  defined  in  Procedure  4.5  and  (dm,Cm)  are  defined  in 
Definition  4.3.  Then,  ak  + bk(dm)  < Cm  for  all  m E {1,  ..,  2 | M(i)  | },  k E {2,  ..., 
K(i)}. 

Proof:  Assume  that  all  the  assumptions  of  Proposition  4.8  hold. 

Claim  1:  a?  + b;2(dm)  < Cm  for  all  m E { 1 , . . ,2  | M(i)  | } 

Steps  1 to  3 of  Procedure  4.5  imply  that  dm  < y/  for  all  m E { 1 . ,2  | M(i)  | }. 
From  the  definitions  of  a/"  and  bk,  a^  + b2  * y/  = Q1  holds.  Since  Q1  is  defined  as  the 
minimum  {Cm  | dm  = y/,  1 < m < 2 | M(i)  | },  the  following  holds: 

aj2  + b2(dm)  < Cm  for  all  m such  that  dm  = y/ (1-1) 

Further,  aj2  + bj2(dm)  < Cm  for  all  m such  that  dm  < y/.  By  contradiction, 
assume  that  there  exists  (dm*,Cm*)  E Sj  such  that  aj2  + b2(dm*)  > Cm*.  Then,  the 
following  must  hold: 


Since  a7  is  defined  as  shown  in  Step  5 of  Procedure  4.5,  (1-2)  contradicts  the 


definition  of  a?.  Thus,  the  following  holds: 

aj2  + bj2(dm)  < Cm  for  all  m such  that  dm  < y*. 


(1-3) 
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Notice  that  dm  < y/  for  all  m.  Thus,  (1-1)  and  (1-3)  prove  Claim  1. 

Claim  2:  a*  + b,k(dm)  < Cm  for  all  m E {1,..,2  | M(i)  | },  3 < k < K(i) 

Suppose  that  k = 3.  Steps  5 and  6 of  Procedure  4.5  imply  that  as3  + bi3(yi2)  = 
&i2  + bj2(y2).  From  Step  5 of  Procedure  4.5,  it  is  easy  to  see  that  there  is  no  m E (1, 
...,  2 | M(i)  | },  which  satisfies  a/*  + h?  * dm  = Cm  and  dm  < y2  at  the  same  time. 
Therefore,  Claim  1 implies  that 

a.2  + b2  * d™  < cm  for  all  (dm5cm)  E S,  such  that  dm  < y2 (2-1) 

From  (2-1),  it  is  easy  to  see  that  a/  < a^  since  a(3  + bi3(yi2)  = af  + bj2(y2) 
holds.  Thus,  the  following  inequalities  hold: 

aj3  + b;3  * dm  < a3  + b 2 * dm  < Cm  for  dm  > y2  (2-2) 

The  first  inequality  holds  since  a(3  < a(2  and  a^  + bj3(y2)  = aj2  + b2(y2),  and 
the  second  inequality  holds  from  Claim  1.  Let  us  consider  the  following  inequality: 

aj3  + b 3 * dm  < Cm  for  dm  < y 2 (2-3) 

The  proof  for  (2-3)  is  similar  to  the  proof  for  (1-2).  Since  each  dm  must  be  equal 
to  or  greater  than  y^'1  (dm  > yL1)  or  less  than  y^'1  (dm  < y^'1),  (2-2)  and  (2-3)  prove 
Claim  2 when  k = 3. 

The  proof  for  Claim  2 when  k > 3 can  be  derived  inductively  from  the  proofs 
of  (2-2)  and  (2-3). 

Since  k = 2,  ...,  K(i),  Claim  1 and  Claim  2 prove  Proposition  4.8.  Q.E.D. 

Theorem  4.9  Validity  of  Procedure  4.5 

Procedure  4.5  produces  a valid  CQ,  where  CQ  is  defined  in  Definition  4.3. 
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Proof:  Assume  that  activity  i has  a nonempty  set  of  modes,  M(i)  = {1,  | M(i)  | }. 

Also  assume  that  a*  and  bk,  k = 2,  K(i),  yk,  k = 1,  K(i),  are  defined  in 

Procedure  4.5  and  S*,  d,  C and  a are  defined  in  Definition  4.3.  It  is  easy  to  see  that  y/ 
= max  {dm  | 1 < m < 2 | M(i)  | } and  y*®  = min  {dm  | 1 < m < 2 | M(i)  | }. 
Let  (B'd^C)  be  an  arbitrary  point,  where  B = [Bj,  ...,  R2\  mo>  | ]*,  0 < Bm  < 1 for  all  m 
G {1,  ..,  2 | M(i)  | } and  B,  + B2  + ...  + B2|M(i)|  = 1. 

Notice  that  Proposition  4.7  proves  the  validity  of  Procedure  4.5  at  zt  = ytl. 
Suppose  that  B‘d  < y/.  Then,  there  exists  some  k*  G {2,  ...,  K(i)}  such  that  yk*  < Bld 
< ykM  since  y*®  < yfm  < ...  < y/  and  y*®  < a'd  < y*1  for  all  a.  Observe  the 
followings: 

af  + b,k*  ( B‘d) 


2|M(i)| 

= 4'  * •>!"  * E » d” 

m=l 

2 |M(i)  | 

E Pm*(>f  *bT*d-) 

m=l 

2|M(i)| 

s E Pm  * (cm ) 

m=l 

= plc  


from  the  definitions  of  8 and  d 

2 |M(i)  | 

SillCe  E Pm  = 1 

m=l 

since  a^  + bk*  * dra  < Cm  for  all  m 

from  Proposition  4.8. 

(a) 


(a)  proves  that  Procedure  4.5  constructs  a function  which  underestimates  CQ 
defined  in  Definition  4.3  for  y(K®  < zt  < y3.  Suppose  that  (zi;  ak  + bk  z)  is  an 
arbitrarily  chosen  point,  where  k G {2,  ...,  K(i)}  and  yk  < z,  < yk_1.  Notice  that 
(yk,CCk)  G Sj  for  k = 1,  ...,  K(i),  where  CCk  is  defined  in  Procedure  4.5.  Thus,  there 
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exists  some  a such  that  a‘d  = zt  and  c/C  = a;k  + bk  zr  Therefore,  Procedure  4.5 
constructs  CQ  defined  in  Definition  4.3  for  y*®  < zx  < yQ  Notice  that  Proposition  4.7 
proves  CCifyj1)  = CC,1  and  CQ(Zj)  is  defined  as  infinity  for  zx  < rnin_dm  (=  y,K(1))  and 
Zj  > rnax_dm  (=  y/).  Thus,  Procedure  4.5  produces  a valid  CQ  defined  in  Definition 
4.3.  Q.E.D. 

Lemma  4.10  CQ  is  a Convex  Linear  Piecewise  Function 

CQ  is  a convex  linear  piecewise  function  whose  domain  is  {zx  E R'  | y*®  < 
Zj  < y/},  where  CQ  is  defined  in  Definition  4.3  and  y>  and  y*®  are  defined  in 
Procedure  4.5. 


Proof:  Theorem  4.9  verifies  that  Procedure  4.5  produces  a valid  CQ.  Assume  that  all 
the  notations  that  are  not  defined  in  this  proof  follow  the  definitions  in  Definition  4.3  and 
Procedure  4.5.  Notice  that  Procedure  4.5  constructs  CQ  such  that 

CC,1  for  z.  = y1 


CQ(Zi)  = ' 


k , k 
ai  + bi  Zi 


+ 00 


for  Yik  ^ Zj  < yf'1,  k = 2,..,K(i)  * - W 
otherwise 


It  is  easy  to  see  that  ak+1  < ak  for  all  k E {2,  ..,  K(i)-1 } and  ajk+1  + bk+1  * yk  = a^ 
+ bk  * yk  for  all  k E {2,  ..,  K(i)-1}.  And  + b3  * y>  = CQ1.  Thus,  CQ  is  a 
convex  linear  piecewise  function  whose  domain  is  {z{  E R1  | ysK®  < zs  < y/}. 

Q.E.D. 


So  far  we  have  considered  activity  cost  functions  having  no  bound  on  an  activity 
duration  and  no  mode  restriction.  Suppose  that  an  activity  duration  must  be  between 
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some  periods.  Or  suppose  that  only  certain  mode  can  be  selected  for  the  processing  of 
an  activity.  We  will  denote  LB(Zj)p  and  UB(Zi)p  as  the  lower  and  the  upper  bounds  of  the 
duration  of  activity  i,  respectively.  We  will  also  denote  M(i)p  as  a set  of  modes  available 
to  activity  i for  the  processing  assuming  that  availability  of  each  mode  is  determined  by 
the  mode  restriction.  M(i)p  is  a subset  of  M(i).  The  mode  cost  with  bounds  and  mode 
restriction  is  defined  as  follows: 

Definition  4.11  Mode  Cost  Function  with  Bounds  and  Mode  Restriction 

Suppose  that  m G M(i).  Assume  that  LB(Zi)p  and  UB(Zj)p  are  the  lower  and  the 
upper  bounds  of  the  duration  of  activity  i,  respectively,  where  0 < LB(Zj)p  < UB(Zj)p. 
Suppose  that  there  exists  mode  restriction  such  that  only  mode  G M(i)p  can  be  chosen 
for  the  processing  of  activity  i.  The  mode  cost  with  activity  duration  bounds  and  mode 
restriction,  C^,  is  defined  as  follows: 


Figure  4.5  shows  Ci(m)(Zj),  where  m G M(i)p,  yi(m)c  < UB(zi)p  < yi(m)nand  LB(Zj)p 

= 0. 

In  a similar  way,  we  can  define  Activity  Cost  Function  with  Duration  Bounds  and 


y i(m)c  ^ ^ y i(m)n 

LB(z/  Zj  s UB(Zj)p 

m e M(i)p 


+ oo 


otherwise. 


Mode  Restriction. 
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Figure  4.5  An  Example  of  Ci(m)(p) 

Definition  4.12  Activity  Cost  Function  with  Bounds  and  Mode  Restriction 

Suppose  that  activity  i has  a set  of  available  modes,  M(i)  = {1,..,  | M(i)  | } and 
a set  of  updated  available  modes,  M(i)p,  where  M(i)p  is  a subset  of  M(i).  Assume  that 
LB(Zi)p  and  UB(Zj)p  are  the  lower  and  the  upper  bounds  of  the  duration  of  activity  i, 
respectively,  where  0 < LB(Zi)p  < UB(Zi)p.  Then,  the  Activity  Cost  Function  of  activity 
i with  duration  bounds  and  mode  restriction  is  defined  as  follows: 

Cj(zi7  xi(l)>  xi(  | M(i)  | ))<P>  — 52  Xi(m)  * ( Q(m)(Zi)  ) 

meM(i) 

where  £ xi(m)  = 1 , xi(m)  € (0,1}  for  all  m E M(i)p,  xi(m)  = 0 for  all  m <£ 

meM(i)p 

M(i)p  and  Ci(m)(p)(zi)  is  defined  in  Definition  4.11. 


One  example  with  LB(Zi)p  = 0 and  M(i)p  = {2}  is  presented  in  Figure  4.6. 
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C® 

1 


Figure  4.6  An  Example  of  C^,  where  M(i)p  = {2} 

The  Convex  Activity  Cost  Function  with  Bounds  is  defined  as  follows: 

Definition  4.13  Convex  Activity  Cost  Function  with  Bounds  and  Mode  Restriction 

Suppose  that  activity  i has  a nonempty  set  of  modes  available,  M(i)  = 
{1,..,  | M(i)  | } and  a set  of  updated  available  modes,  M(i)p,  where  M(i)p  is  a subset  of 
M(i).  Assume  that  LB(Zj)p  and  UB(Zi)p  are  the  lower  and  the  upper  bounds  of  the 
duration  of  activity  i,  respectively,  where  0 < LB(Zj)p  < UB(Zi)p.  Let  Sj’  be  a set  in  R2 
such  that  S(’  = {(min  {yi(m)n,  UB(Zi)p},  Ci(m)(rnin  (yKm)n,  UB(Zi)p}))  | m E M(i)p  and 
Ci(m)(rnin  {yi(m)n,  UB(Zi)p})  is  finite  } U {(max  {yi(m)c,  LB(Zi)p},  Ci(m)(max  {yi(m)c,  LB(Zi)p})) 

| m E M(i)p  and  Cj(m)(rnax  {yi(m)c,  LB(Zi)p})  is  finite  }.  Assume  that  each  element  in 
S/  is  relabeled  such  that  S*’  = {(d^C1),  (d2,C2),  (dJ,Cj),  (dJ,CJ)},  where  J = 

| S;’  | . Then,  the  Convex  Activity  Cost  Function  of  activity  i with  duration  bounds  and 
mode  restriction,  CC^,  is  defined  as  follows: 
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CCj(z.)  = 


min  | a*  C | a1  d = Zj  } 

for  min_dJ  <.  <.  max_d) 

+ OO  otherwise. 


where  max_dj  = max  {dj  | 1 < j < J},  min_dj  = min  {dj  | 1 < j < J},  a = [a„ 
oc2,  ...  , ttj]1,  0 < ckj  < 1 for  j = 1,  J,  a,  + a2  + ...  + a,  = 1,  d = [d1,  dJ]‘ 
and  C = [Cj,  CJ]‘. 


We  defined  the  lowest  cost  of  completing  one  activity  with  a particular  mode  as 
the  mode  normal  cost  and  the  corresponding  duration  as  the  mode  normal  duration.  And 
we  defined  the  shortest  duration  of  completing  an  activity  with  a given  mode  as  the  mode 
crash  duration  and  the  corresponding  activity  cost  as  the  mode  crash  duration.  When  we 
defined  CQ  in  Definition  4.3,  we  used  Si5  which  is  a collection  of  each  mode’s  (normal 
duration,  the  corresponding  mode  cost)  and  (crash  duration,  the  corresponding  mode 
cost).  When  there  exist  activity  duration  bounds,  LB(z,)p  and  UB(Zi)p,  and/or  mode 
restriction,  the  lowest  cost  of  completing  activity  i with  mode  m is  Cj(m)(rnin 
{UB(Zi)p,yi(m)n})  if  UB(Zj)p  > yi(m)c,  LB(Zj)p  < yi(m)n  and  m G M(i)p.  In  a similar  way, 
the  shortest  duration  of  completing  activity  i with  mode  m is  max  {LB(zi)p,  yj(m)c}  if 
UB(Zi)p  > yi(m)c,  LB(zi)p  < yi(m)n  and  m G M(i)p.  Notice  that  if  UB(Zi)p  < yi(m)c  or 
LB(Zi)p  > yi(m)n,  or  if  m is  not  in  M(i)p,  then  there  exists  no  zi  whose  mode  cost  is  finite. 
It  is  obvious  that  if  UB(Zi)p  > yi(m)c,  LB(Zj)p  < yj(m)n  and  m G M(i)p,  then  min 
{UB(zi)p,yi(m)n}  and  max  {LB(Zj)p,  yi(m)c}  are  the  longest  and  shortest  durations  with  finite 
mode  cost,  respectively.  In  defining  CC^,  we  use  S ’,  which  is  a collection  of  each 
mode’s  (the  longest  duration  with  finite  mode  cost,  the  corresponding  mode  cost)  and 
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(the  shortest  duration  with  finite  mode  cost,  the  corresponding  mode  cost).  Suppose  that 
a mode,  m*,  is  in  M(i),  but  not  in  M(i)p.  Then,  mode  m*  does  not  provide  any  point 
to  Sj’  since  Ci(mt)(p)(z^)  is  always  infinite.  If  a mode  6 M(i)p  has  no  duration  in  which 
the  mode  cost  is  finite,  then  the  mode  does  not  provide  any  point  to  S'. 

Figure  4.7  illustrates  CC^,  where  M(i)  = M(i)p  = {1,2},  LB(Zi)p  = 0 and 

UB(Zi)p  = Zj*. 


Figure  4.7  An  Example  of  CQ®,  where  M(i)p  = {1,2} 

The  procedure  for  constructing  CCi(p>  is  identical  to  Procedure  4.5  except  for  Step 
1.  The  modified  Step  1 of  Procedure  4.5  is  as  follows: 

Modification  4.14  The  Modified  Step  1 of  Procedure  4.5 
Step  1.  For  each  mode  in  M(i)p, 

if  UB(Zj)p  > yi(m)c  and  LB(Zj)p  < yi(m)n,  then 

let  (dml,  Cml)  = (min  (UB(Zi)p,  yi(m)n),  Ci(m)(dm1))  and 


(dm2,  Cm2)  = (max  (LB(Zj)p,  yi(m)c),  C,(m)(dm2)). 
if  UB(Zi)p  < yi(m)c  or  LB(z,)p  > yi(m)n,  then  do  nothing. 
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The  proof  for  the  validity  of  Procedure  4.5  with  Modification  4.14  is  similar  to 
the  Proof  of  Theorem  4.9. 

4.2  The  Branch  and  Bound  Procedure 

We  propose  a branch  and  bound  solution  procedure  for  our  problem.  In  our 
proposed  branch  and  bound  procedure,  node  0 refers  to  (P).  Let  us  denote  the  node 
problem  of  node  p,  feasible  region  of  the  node  p and  the  optimal  solution  value  of  node 
problem  of  node  p as  (Pp),  F(PP)  and  V(PP),  respectively.  (P°)  is  (P).  Consider  some 
arbitrary  nodes  A and  B on  the  branch  and  bound  tree.  If  there  is  a branch  from  node 
A to  node  B,  then  we  will  say  that  node  A is  the  parent  node  of  B and  that  node  B is  a 
child  node  of  A.  A node  may  have  multiple  children  nodes  (or  a null  child  node),  but 
each  node  has  just  one  parent  node,  except  node  0.  P-Node(p)  will  denote  the  parent 
node  of  node  p,  while  C-Node(p)  will  denote  the  set  of  children  nodes  of  node  p. 

Suppose  that  we  are  at  some  node  p.  Since  (Pp)  is  not  easy  to  solve,  we  will 
solve  a relaxation  of  (Pp),  which  we  will  denote  as  (RPP).  Suppose  that  we  derive  a 
solution  of  (Pp)  from  the  optimal  solution  of  (RPP).  If  the  solution  of  (Pp)  satisfies  all  the 
constraints  of  (Pp),  then  we  have  solved  (Pp).  If  not,  then  we  will  branch  (or  create 
children  nodes)  from  node  p such  that  F(PP)  = U ( F0”)  } . As  long  as  each 

qeC-Node(p) 

child  node  q problem,  q E C-Node(p),  has  the  same  objective  function  of  its  parent  node 
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p problem,  and  if  F(PP)  = (J  { F(Pq)  } holds,  then  V(PP)  = 

qeC-Node(p) 

minimum { V(Pq)  } holds.  In  other  words,  (Pq)  is  defined  as  a node  problem  that  has 

qeC-Node(p) 

the  same  objective  function  as  (Pp),  and  the  feasible  space,  F(Pq),  is  defined  as  shown 
above.  Now,  if  we  fail  to  derive  V(PP)  by  solving  its  relaxed  problem,  then  solving  (Pp) 
turns  into  solving  all  of  the  (Pq),  where  q E C-Node(p).  This  process  is  repeated  until 
we  solve  (P°),  which  is  (P). 

In  section  4.2.1,  three  node  problems  are  defined,  and  the  conditions  for  deriving 
V(PP)  are  explained.  The  branching  rules  that  are  used  when  we  fail  to  find  V(PP)  are 
explained  in  section  4.2.2,  and  our  complete  branch  and  bound  algorithm  is  presented 
in  section  4.2.3. 

4.2.1  Node  P Problems  and  Optimality  Conditions 

At  each  node,  a node  problem  is  defined,  and  its  relaxed  problem  is  solved.  We 
will  use  the  relaxed  problem  solution  to  test  whether  we  can  derive  the  optimal  solution 
of  the  node  problem.  First,  the  problem  at  node  p,  (Pp),  is  defined  as  follows: 

Definition  4.15  Problem  at  Node  p,  (Pp) 

The  Problem  at  node  p,  (Pp),  is  defined  as  follows: 

i 

(Pp)  min  ^2  xi(m)  * | NCi(m)  + ci(m)  yi(m)n  - ci(m)zi  J 

i=l  meM(i) 


+ xI+1  * Penalty  * ( FI+1  - DUE_DATE  ) 


(1) 
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st 


where, 


Fj  + Zj  < Fj, 

(ij)  e Hp 

(2) 

XI  Xi(m)ri(m)k  ^ \ » 

ieSt 

for  all  k and  t 

(3) 

•n 

o 

II 

O 

(4) 

XI  Xi(m)  * ’ 

meM(i)p 

i = 1,  ..., 

I 

(5) 

xKm)  G {°4}, 

i = 1,  ..., 

I,  m E M(i) 

(6) 

^i(m>y  i(m)c  — ^i(m)^i  — ^i(m)y  i(m)n> 

i = 1,  ..., 

I,  m E M(i) 

(7) 

Fj  integer. 

i = 0,  ..., 

1+1 

(8) 

Zj  integer. 

i = 1,  ..., 

I 

LB(Zi)p  < Zj  < UB(Zi)p, 

i = 1,  ..., 

I 

(9) 

LB(Fi)p  < Fj  < UB(Fi)p, 

i = 1,  ..., 

1+1 

(10) 

f 0 if  FI+1  <;  DUE_DATE 
1+1  [ 1 otherwise 

(ID 

o 

II 

(12) 

Zi+I  — 0 

j* 

f 

II 

o 

for  all  m 

£ M(i)p  for  all  i 

(13) 

Hp:  set  of  ordered  pairs  of  activities  indicating  the  precedence 

relationship  at  node  p.  HPNode(p>  is  a subset  of  Hp  for  all  p, 

p ^ 0. 
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M(i)p:  the  set  of  modes  available  to  activity  i at  node  p.  M(i)p  is 

a subset  of  M(i)P  Node(p)  for  all  p,  p ^ 0. 

LB(Zi)p:  the  lower  bound  of  at  node  p.  LB(Zi)p  > LB(zj)PNode(p) 

for  all  p,  p ^ 0. 

UB(Zi)p:  the  upper  bound  of  Zj  at  node  p.  UB(zi)p  < UB(zi)P  Node(p) 

for  all  p,  p ^ 0. 

LB(Fi)p:  the  lower  bound  of  F;  at  node  p.  LB(Fi)p  > LB(Fi)p"Node(p) 

for  all  p,  p ^ 0. 

UB(Fi)p:  the  upper  bound  of  Fj  at  node  p.  UB(Fj)p  < UB(Fi)PNode(p) 

for  all  p,  p / 0. 

Other  notations  follow  the  definitions  given  in  section  3.2. 

Node  0 is  called  the  origin  node,  and  (P°)  is  identical  to  the  Problem  defined  in 
section  3.2  such  that  H°  = H,  M(i)°  = M(i),  LB(Zi)°  = 0,  UB(zi)°  = + oo,  for  all  i £ 
{1,  ...,  I},  LB(Fi)°  = 0,  and  UB(Fj)°  = + oo  for  all  i.  The  problem  at  node  p,  p ^ 
0,  and  the  problem  at  node  0 differ  in  precedence  relations,  bounds  on  each  activity 
duration  and  a set  of  mode  available  for  each  activity  at  the  node.  Let  node  q be  a child 
node  of  node  p.  Then,  Hq  is  a union  of  Hp  and  a set  of  additional  precedence  relations, 
if  any.  LB(Zj)q  > LB(Zj)p,  and  UB(Zj)q  < UB(Zi)p.  M(i)q  is  a subset  of  M(i)p.  Since 
M(i)p  consists  of  modes  available  for  activity  i at  node  p,  only  modes  in  M(i)p  can  be 
chosen  for  the  processing  of  activity  i at  node  p.  Therefore,  all  xi(m)  such  that  m 


£ M(i)p  will  be  set  to  0.  This  restriction  is  presented  in  constraint  (13).  The  following 
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Proposition  shows  that  the  optimal  solution  value  of  a child  node  problem  is  greater  than 
or  equal  to  the  optimal  solution  value  of  the  parent  node  problem  if  there  exists  a feasible 
solution  of  the  child  node  problem. 

Proposition  4.16  V(Pq)  > V(PPNode(q)) 

Suppose  that  node  q is  a child  node  of  node  p.  Assume  that  the  node  q problem 
has  a feasible  solution.  Then,  the  optimal  solution  value  of  the  node  q problem  is  greater 
than  or  equal  to  the  optimal  solution  value  of  the  node  p problem,  where  node  problem 
is  defined  in  Definition  4.15. 

Proof:  Assume  that  all  the  assumptions  of  Proposition  4.16  hold.  It  is  easy  to  see  that 
the  objective  function  of  the  node  q problem  is  identical  to  that  of  the  node  p problem. 
From  the  definitions  of  Hq,  M(i)q,  LB(Zj)q,  UB(Zi)q,  LB(Fj)q  and  UB(Fj)q,  it  is  easy  to  see 
that  any  feasible  solution  of  the  node  q problem  is  a feasible  solution  of  the  node  p 
problem.  This  proves  Proposition.  Q.E.D. 

Since  the  problem  at  node  p is  hard  to  solve,  we  construct  the  relaxed  problem 
at  node  p that  can  be  solved  with  moderate  effort  and  that  provides  a lower  bound  for 
the  V(Pp). 

Definition  4.17  The  Relaxed  Problem  at  Node  p,  (RPP) 

The  relaxed  problem  at  node  p,  (RPP),  is  defined  from  (Pp)  by  dropping  the 
resource  and  integrality  constraints  and  by  replacing  the  activity  cost  function  with 
bounds  and  mode  restriction,  Q®,  with  the  convex  activity  cost  function  with  bounds  and 
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mode  restriction,  CC^.  Also  all  the  constraints  that  contain  mode  variable,  xi(m),  are 
deleted. 


i 

(RPP)  min  CCf^Zj) 

i=l 

+ xI+1  * Penalty  * ( F]+1  - DUE_DATE  ) 
st  Fj  + Zj  < Fj,  (i,j)  E Hp 

F0  = 0 

LB(z/  < z,  < UB(Zi)p,  i = 1,  I 

LBCFJp  < Fi  < UB(Fj)p,  i = 1,  ...,  1 + 1 

( 0 if  Fn+1  <;  DUE_DATE 

Xn+1  \ 1 otherwise 

Zo  = 0 

Zi+i  = 0 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 
(7) 


Since  CQ®  < Q for  all  durations  and  for  all  mode  selections  (xi(1),  ...,  xi(  ( M(i)  | }) 
for  each  activity  and  the  constraints  of  (RPP)  is  a subset  of  constraints  of  (Pp),  V(RPP)  < 
V(PP)  follows.  Thus,  V(RPP)  provides  a lower  bound  for  V(PP).  James  Kelly  [1961] 
showed  how  to  handle  activity  utility  functions  that  are  piecewise  linear,  nondecreasing 
and  concave,  where  all  the  constraints  are  linear  and  the  objective  is  the  maximization 
of  the  utility  of  a project.  He  reformulated  the  activity  utility  functions  so  that  the 
problem  is  reduced  to  a linear  program.  Notice  that  the  objective  of  the  relaxed  problem 
is  the  minimization  of  the  sum  of  the  convex  activity  cost  and  penalty  cost  and  that  each 
convex  activity  cost  function  is  piecewise  linear.  Also  notice  that  all  the  constraints  are 
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linear  except  constraint  (6)  that  is  relevant  to  the  penalty  cost  function.  Note  that  the 
penalty  cost  function  can  be  expressed  as  maximum  {0,  Penalty  * FI+1  - Penalty  * 
DUEDATE}  for  all  FI+1  > 0.  Since  Penalty  is  assumed  to  be  positive,  the  penalty  cost 
function  is  a convex  linear  piecewise  function.  One  example  is  presented  in  Figure  4.8. 


Penalty  Cost 


Figure  4.8  Penalty  Cost  Function 

Thus,  we  can  reduce  the  relaxed  problem  into  a linear  programming  problem. 
The  relaxed  problem  at  node  p in  a linear  programming  formulation  is  called  a linear 
programming  problem  at  node  p,  (LPP),  and  the  formulation  is  presented  in  Definition 
4.18. 

Definition  4.18  Linear  Programming  Problem  at  Node  p,  (LPP) 

The  linear  programming  problem  at  node  p,  (LPP),  is  constructed  by  transforming 
a relaxed  problem  at  node  p into  a linear  programming  problem. 
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(LPP)  min 


i 

E 

i=l 


cc 


K(i)  - 1 

E 

k=l 


k k 
Ci  Z: 


+ Penalty  * (BDUR  - DUEDATE  - zpena)ty) 

(1) 

Fj  + Zj  < Fj,  (ij)  e Hp 

(2) 

F0  = 0 

(3) 

K(i) — 1 

z,  = y l-  E 1 

k=l 

(4) 

0 < zk  < yk  - yk+1  k = 1,  K(i)  - 1,  i = 1,  . 

I 

(5) 

LB(Zi)p  < z,  < UB(Zi)p  i = 1,  I 

(6) 

LB(Fi)p  < F,  < UB(Fi)p  i = 1,  1 + 1 

(7) 

o 

II 

c? 

(8) 

Zi+i  — 0 

Fi+i  ^ B DUR  - zpenally 

(9) 

0 < zpenalty  < B DUR  - DUE  DATE 

(10) 

where  CC/: 


convex  activity  function  value  at  the  first  breakpoint  (the  longest 
duration)  on  CC,® 

duration  value  at  the  k-th  breakpoint  (k-th  longest  duration)  on 


CC,® 

K(i):  the  number  of  breakpoints  in  CC,® 

ck:  the  marginal  crashing  cost  of  the  duration  between  the  k-th  and  the 

(k+l)-th  breakpoint  in  CC,®  such  that 
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k 


Ci 


pk+1  /-i  k 

hi L_ri_  fork  = 1 K(i)-1. 

k+1  k 


Z|k:  crashing  variable  of  activity  i over  the  interval  between  yk+1  and 

yik. 

B DUR:  an  arbitrary  period  long  enough  that  no  project  completion  time 

will  be  later  than  the  period.  An  example  for  B DUR  is  the  sum 
of  the  longest  durations  (with  finite  cost)  of  each  activity  in  the 
project. 

zpenod:  crashing  variable  for  penalty  cost  over  the  interval  between 

DUE  DATE  and  B DUR. 


The  LP  formulation  of  the  convex  activity  cost  function  can  be  easily  understood 


by  a small  example.  Let  us  consider  the  same  activity  shown  in  Example  4.6.  The 
summary  for  the  CCp*  is  as  follows: 


k (yk,  Ck) 


1 (10,30) 

2 (7,32) 

3 (3,41) 


CCjCzj)® 


30.00  for  Zj  = 10 

46.00  - 0.6667  * z{  for  7 £ < 10 

47.75  - 2.2500  * for  3 £ z.  < 7 


+ oo 


otherwise 
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Since  activity  i has  three  breakpoints,  CC&  contains  two  intervals.  Thus,  we 
need  two  nonnegative  crashing  variables,  z/  and  z2.  The  length  of  each  interval  is  the 
difference  between  breakpoints,  z'  and  z2  are  bounded  as  follows: 


0 < Z;1  < ; 

y > . y2  = 10  - 7 = 3 

0 < Z:2  < 

y2  . y3  = 7 - 3 = 4. 

Each  crashing  variable  has  an  associated  crash  cost.  For  example,  z^’s  crash  cost 
is  0.6667,  and  z^’s  crash  cost  is  2.25.  Let  us  express  the  activity  i duration,  zh  as  y/  - 
z2  - z2.  If  Zj  = y/  (=  10  = the  longest  duration),  then  z'  and  z7  must  be  zero  since 
both  are  nonnegative  (no  crash  occurs),  and  the  cost  is  set  to  CC,1  (the  lowest  one). 
Suppose  that  for  some  reason,  z-,  must  be  crashed  by  one  period  from  y*.  Since  z{  = y/ 
- z^  - z2,  this  implies  that  the  value  z/  or  z2  must  be  increased  from  0.  Since  our 
problem  is  a minimization  problem,  LP  always  picks  the  crashing  variable  whose  crash 
cost  is  the  lowest  before  choosing  others.  This  is  repeated  until  z,  reaches  7,  in  other 
words,  z/  becomes  3.  If  further  crashing  is  needed,  then  z2  is  selected  since  z'  cannot 
be  increased  any  more.  This  is  repeated  until  zt  = y2  (=  4 = the  shortest  duration). 
In  summary,  CC^  of  Example  4.6  is  expressed  in  LP  formulation  as  follows: 

CCi(zi)(p)  = 30  + 0.6667  * z>  + 2.2500  * Zj2 
where  z{  = 10  - z/  - z2,  0 < z/  < 3 and  0 < z2  < 4. 

Figure  4.9  is  provided  for  the  above  example. 

The  Penalty  Cost  Function  in  an  LP  formulation  can  be  understood  in  a similar 
way.  But,  unlike  the  convex  activity  cost  function,  the  first  breakpoint,  (the  longest 
project  duration  with  finite  penalty  cost,  its  corresponding  cost),  does  not  exist.  Thus, 
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(p) 

i 


Figure  4.9  CQ®  in  LP  Formulation 

we  consider  a sufficiently  long  project  duration,  BDUR,  and  calculate  the  corresponding 
cost,  which  is  Penalty  * (B  DUR  - DUEDATE).  Now,  we  have  three  breakpoints  in 
the  penalty  cost  function— (B  DUR,  Penalty*(B_DUR-DUE_DATE)),  (DUE  DATE,  0) 
and  (0,0).  Thus,  two  crashing  variables,  z^^1  for  the  interval  between  B DUR  and 
DUE  DATE  and  zpenalty2  for  the  interval  between  0 and  DUE  DATE,  are  introduced,  and 
Fn+1  is  set  to  B_DUR  - zpenaity'  -zpenalty2.  Notice  that  the  marginal  crashing  cost 
corresponding  to  z^,^1  is  -Penalty  and  that  the  marginal  crashing  cost  corresponding  to 
penalty2  is  0-  Thus,  the  Penalty  Cost  function  can  be  expressed  as  follows: 

Penalty  Cost  Function  (FN+1) 

= Penalty  * (B  DUR  - DUE  DATE)  + (-  Penalty)  * z^1  + 0 * zpenalty2 
= Penalty  * (B  DUR  - DUE  DATE  - z^1) 
where  FN+1  = B DUR  - z^,^1  - zpenaHy2,  0 < zpcniiity’  < (B  DUR  - DUE  DATE)  and  0 
< zpenalty2  < DUE_DATE.  More  compact  formulation  for  the  penalty  cost  function  can 
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be  made  after  dropping  zpenalty2  and  by  replacing  zpenalty‘  with  Zpena)ty.  Since  zpenaity2  > 0, 
Fn+i  = BDUR  -Zpe^1  - zpenalty2  implies  that  FN+1  > BDUR  -Zp^1  = BDUR  -zpenalty. 
Thus,  the  new  formulation  for  the  penalty  cost  function  is  as  follows: 

Penalty  Cost  Function  (FN+1) 

= Penalty  * (B  DUR  - DUEDATE  - zpenalty) 
where  FN+1  > B DUR  - zpenalty  and  0 < zpenalty  < (B  DUR  - DUE  DATE).  As  an 
example,  Figure  4. 10  is  provided. 

Notice  that  (RPP)  is  a relaxed  problem  of  (Pp)  and  that  (LPP)  is  an  LP  formulation 
of  (RPP).  Thus,  it  is  easy  to  see  that  V(LPP)  < V(PP).  Since  (LPP)  is  a linear 
programming  problem,  it  can  be  solved  with  moderate  effort.  Suppose  that  V(LP)P  is 
infinite.  Then,  we  conclude  that  (Pp)  has  no  feasible  solution  since  V(PP)  must  be 
infinite.  Thus,  we  solved  (Pp).  Suppose  that  V(LPP)  is  finite.  Then,  (Pp)  may  have 
feasible  solutions,  and  we  test  whether  the  optimal  solution  of  (Pp)  can  be  derived.  The 
sufficient  conditions  for  derivation  of  the  optimal  solution  of  (Vp)  from  the  solution  of 
(LPP),  or  equivalently  from  the  solution  of  (RPP),  are  as  follows: 

Theorem  4.19  Sufficient  Conditions  for  the  Derivation  of  the  Optimal  Solution 

Suppose  that  z*  and  Fj*  are  the  optimal  solution  values  for  z,  and  Fb  i =0,  ..., 
1+1 , of  (LPP)  assuming  that  V(LPP)  is  finite.  If  the  following  conditions  are  satisfied, 
then  V(PP)  = V(LPP),  and  an  optimal  solution  of  (Pp)  can  be  found. 

(i)  Zj*  and  Fj*,  i = 1,  ...,  I + 1 , are  all  integers. 

(ii)  For  each  i E (1,  ..,  I},  there  exists  a mode  m*  E M(i)p  such  that  CC^fo*)  = 

Cio.-.nz*). 
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Penalty  Cost 


Penalty  * 
(B_DUR  - 
DUE_DATE ) 


DUEDATE 


B_DUR 


Figure  4.10  Penalty  Cost  Function  in  LP  Formulation 


(iii)  For  all  k and  t,  V r.  £ b.  holds,  where  m*  is  a mode  chosen  from  (ii). 

1(111  *)K  K 


Proof:  Suppose  that  the  assumptions  and  conditions  (i)  to  (iii)  of  Theorem  4. 19  hold. 
Let  us  consider  the  following  solution  of  (Pp):  (zq*,  z,*,  ...,  zI+1*,  F0*,  Fj*,  ...,  FI+1*, 


of  the  optimal  solution  of  (LPP)  and  that  Xi(m*)*  = 1 for  i = 1,  I if  m*  is  a mode 
chosen  from  condition  (ii)  of  Theorem  4. 19  and  xi(m)*  = 0 if  m ^ m*  for  m E M(i)  and 
i = 1,  ..,  I.  Assume  that  xI+1*  is  0 if  FI+1*  < DUE  DATE  and  1 otherwise.  Since 
constraints  (2),  (4),  (9),  (10)  and  (12)  of  (Pp)  are  explicitly  considered  in  (LPP),  the 
solution  satisfies  those  constraints.  Since  condition  (ii)  of  Theorem  4.19  considers  m* 
E M(i)p,  i = 1,  I,  and  since  we  set  xi(mt)*  = 1,  while  setting  xi(m)*  = 0,  m ^ m* 


s, 


*).  Assume  that  z*  and  Fj*,  i = 0,  ...,  1 + 1,  are  components 


for  m E M(i)  and  i = 1,  ..,1,  constraints  (5),  (6)  and  (13)  of  (Pp)  are  satisfied  by  the 
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solution.  From  conditions  (i)  and  (iii),  constraints  (3)  and  (8)  are  satisfied.  Since  x1+1* 
is  set  to  0 if  FI+1*  < DUE_DATE  and  1 otherwise,  the  constraint  (11)  is  satisfied.  Let 
us  consider  any  mode  variable  from  the  solution.  Suppose  that  xi(m)*  = 0.  Then, 
constraint  (7)  is  satisfied  since  xi(m)*  yi(m)c  < xi(m)*z*  < xj(m)*yi(m)n  are  reduced  to  0 < 
0 < 0.  Suppose  that  xi(m.)*  = 1.  Since  V(LPP)  is  assumed  to  be  finite,  CC^Xz*)  must 
be  finite.  Condition  (ii)  implies  that  CCi(p)(z,*)  = CKm.)(p)(Zi*).  Thus,  Ci(m«)(p)(zi*)  is 
finite.  From  the  definition  of  as  shown  in  Definition  4.12,  yi(m*)c  < z*  < yi(m*)n 
must  hold  for  Cj^^Zj*)  to  be  finite.  Thus,  constraint  (7)  is  satisfied.  Since  all  the 
constraints  of  (Pp)  are  satisfied  by  the  solution,  the  solution  is  a feasible  solution  of  (Pp). 


From  condition  (ii)  of  Theorem  4. 19  and  from  the  values  of  xi(m)*  and  x^*,  the 
following  equalities  hold: 


(a) 


E E 


i=l  meM(i) 


+ (xI+1*)  * Penalty  * (FI+1*  - DUEDATE) 


£ Ci(m*)(zi*)  + (xI+1‘)  * Penalty  * (FI+1*  - DUE  DATE) 


i=l 


52  Ci(m^\z*)  + (x1+1*)  * Penalty  * (FI+1*  - DUE  DATE) 


i=l 


since  LB(Zi)p  < z * < UB(Zi)p  and  m*  6E  M(i)p  for  i E {1,..,!} 


Y CCi(p)(zi*)  + (xI+1*)  * Penalty  * (FI+1*  - DUE  DATE) 


i=l 


from  condition  (ii)  of  Theorem  4.19 


= V(RPP)  = V(LPP)  since  (LPP)  is  an  LP  formulation  of  (RPP) 
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Thus,  the  solution  value  of  (Pp)  with  the  solution  is  V(RPP)  = V(LPP).  This  and 
(a)  imply  that  the  solution  is  an  optimal  feasible  solution  of  (Pp)  and  that  V(PP)  is 
V(LPP).  Q.E.D. 

4,2.2  Optimality  Condition  Checking  Procedures  and  Branching  Rules 

In  this  section,  we  present  procedures  that  check  the  conditions  listed  in  Theorem 
4.19  and  create  children  nodes  from  the  current  node  if  any  condition  is  not  satisfied. 
Suppose  that  we  are  at  node  p and  some  of  conditions  in  Theorem  4.19  are  violated. 

Then,  we  will  create  children  nodes  such  that  F(PP)  = U { F(Pq)  } • (Pp)  and 

qeC-Node(p) 

(Pq),  for  some  q E C-Node(p),  are  identical  except  for  one  of  M(i),  H,  LB(Zi),  UB(Zi), 
LB(Fi)  and  UB(F;).  If  (Pp)  and  (Pq)  differ  in  one  of  M(i),  H,  LB(Zi),  UB(Zj),  LB(F;)  and 
UB(Fi),  then  we  will  say  that  the  component  is  updated  in  (Pq).  For  example,  if  M(i)q 
differs  from  M(i)p,  then  we  will  say  that  M(i)q  is  updated.  Let  us  consider  each  condition 
in  Theorem  4.19. 

The  first  condition  in  Theorem  4.19  is  whether  z*  and  F;*,  i = 0,  ...,  1 + 1, 
from  the  optimal  solution  of  (LPP)  are  all  integers  or  not.  Inspection  for  this  is  trivial. 
Notice  that  Zq  and  zI+1  are  forced  to  be  zero  by  constraint  (8)  of  (LPP)  and  that  F0  is 
forced  to  be  zero  by  constraint  (3)  of  (LPP).  Suppose  that  z*  for  some  i E {1,  ..,  1}  is 
not  an  integer.  Since  Zj*  is  a component  of  an  optimal  solution  of  (LPP),  z * satisfies 
constraint  (6)  of  (LPP),  which  is  LB(Zi)p  < z{  < UB(Zi)p.  Since  constraint  (8)  of  (Pp) 
requires  an  integer  value  of  zb  there  is  no  feasible  solution  of  (Pp)  with  z(  such  that  [z  *] 
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< Zj  < [Zj*]  + 1,  where  [x]  denotes  the  largest  integer  less  than  or  equal  to  x.  This 
implies  that  the  feasible  interval  of  zs  in  (Pp)  can  be  partitioned  as  follows: 

(i)  LB(Zj)p  < ^ < [z^]  and 

(ii)  [Zj*]  + 1 < ^ < UB(zi)p. 

In  a similar  way,  if  F+  for  some  i G {1,  1 + 1}  is  not  an  integer,  we  can 

partition  the  feasible  interval  of  Fj  as  follows: 

(i)  LBtFJP  < F(  < [F  *]  and 

(ii)  [F*]  + 1 < Fs  < UB(Fi)p. 

If  some  z*  or  F+  is  not  an  integer,  we  choose  one  among  non  integer  components 
and  create  two  children  nodes,  ql  and  q2.  If  z * is  chosen,  then  we  update  UB(z;)ql  = 
[z,*]  and  LB(Zi)q2  = [z+]  + 1.  If  F+  is  chosen,  then  we  update  UB(Fj)ql  = [F+]  and 
LB(Fi)q2  = [F+]  + 1.  Since  F(LPq)  is  bounded  by  LB(zi)q  and  UB(Zj)q  and  by  LB(Fj)q 
and  UBIFi)*1,  the  optimal  solution  of  (LPP)  cannot  reappear  as  an  optimal  solution  of 
(LPql)  or  (LPq2). 

The  second  condition  in  Theorem  4. 19  is  whether  there  exists  an  m*  G M(i)p  that 
satisfies  CC-^\z  *)  = Ci(mt)(p)(z  *).  Since  CC^  for  all  i and  Ci(m)(p)  for  all  i and  m are 
defined  for  all  nonnegative  durations,  CC^Zi*)  and  Cj^^Zi*)  are  readily  available.  If 
multiple  modes  satisfy  condition  (ii)  in  Theorem  4. 19,  then  choose  one  arbitrarily  among 
them.  Suppose  that  CC^Zj*)  ^ Ci(m)(p)(zi5t!)  for  any  m G M(i)p.  Since  CC^  < Ci(m)(p) 
for  all  m G M(i)p  and  for  all  zx  > 0,  CCi(p)(zi*)  ^ Ci(m)<p)(zi5ic)  for  any  m G M(i)p  implies 
that  the  solution  value  of  (Pp)  with  the  solution  derived  from  the  solution  of  (LPP)  as 
shown  in  Theorem  4.19  must  be  greater  than  V(LPP).  Then,  we  cannot  verify  that  the 
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derived  solution  is  an  optimal  solution  of  (Pp)  since  the  solution  value  of  (Pp)  with  the 
derived  solution  is  strictly  bigger  than  V(LPP).  One  example  for  this  case  is  shown  in 
Figure  4.11. 


CCjtP) 


Figure  4.11  CC^z*)  < Cj^z,*)  for  all  m E M(i)p 

We  will  say  that  there  exists  an  underestimation  gap  for  activity  i at  duration  z-* 
if  CC^Zj*)  < Cjdn^Zj*)  for  all  m E M(i)p.  If  an  underestimation  gap  exists,  we  will 
create  two  children  nodes,  ql  and  q2,  and  update  UB(Zj)ql  = z*  and  LB(z,)q2  = z *. 

The  following  two  Figures,  Figure  4.12  and  Figure  4.13,  illustrate  the  effect  of 
this  branching  operation.  If  we  compare  Figure  4.11  to  Figure  4.12  and  4.13,  at  any 
duration  between  the  activity  duration  bounds,  Q00  and  Ci(m)(p)  have  not  changed,  but 
CCi(p)  has  changed  to  be  a tighter  underestimation  function  by  partition. 
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Cost 


c i(l)  (yi(l)c) 
ci(2)  (yi(2)c) 

ci(l)  (yi(l)n) 
ci(2)  & i(2)n ) 


Function 

C id) 

Ci(2) 

> 

yi(l)c  yi(2)c  yi(l)n 

y i(2)n 

LB(Zi)  = zj' 

Figure  4.12  Activity  Cost  Functions  with  Z;  > z ■* 


Cost 


Figure  4.13  Activity  Cost  Functions  with  zs  < z ■* 

The  third  condition  in  Theorem  4.19  checks  feasibility  with  respect  to  the 
resource  constraints.  This  checking  procedure  is  performed  if  conditions  (i)  and  (ii)  of 
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Theorem  4.19  are  met.  Note  that  the  solution  from  (LPP)  provides  F * and  z*  for  each 
activity  and  that  m*  for  each  activity  is  found  during  the  checking  process  for  (ii)  of 
Theorem  4.19.  Feasibility  of  the  resource  constraints  can  be  done  by  investigating  all 
the  resource  requirements  at  the  start  time  of  each  activity  in  the  project  since  the 
resource  requirements  increase  when  an  activity  starts  to  be  processed.  If  the  resource 
requirements  do  not  exceed  the  capacity  of  the  resources  at  the  start  time  of  each  activity, 
then  we  have  an  optimal  solution  for  (Pp).  If  the  resource  requirements  exceed  the 
capacity  of  any  resource  at  the  start  time  of  some  activity,  then  we  will  say  that  there 
exists  a resource  conflict  at  the  start  time  of  the  activity.  A resource  conflict  set  for  the 
period,  RC,  is  defined  as  a set  of  activities  in  progress  during  the  period.  If  RC  exists, 
then  we  have  failed  to  derive  a feasible  solution  with  Fs*,  z*  and  m*.  In  order  to 
develop  a branching  rule  for  this  case,  we  need  the  following  Characterization. 

Characterization  4.20  The  Characterization  of  a Feasible  Solution  with  RC(A,t) 

Suppose  that  we  have  a schedule,  A,  of  the  project,  where  each  activity’s  mode 
selection  and  duration  are  restricted  by  M(i)p  and  by  LB(Zj)p  and  UB(Zj)p,  respectively, 
for  i = 1,  ...,  I.  Assume  that  at  some  period,  t,  the  resource  requirements  are  greater 
than  the  capacity  of  some  resource.  Let  RC(A,t)  be  a set  of  activities  in  progress  during 
the  period  t in  the  schedule  A.  Then,  any  feasible  solution  of  (Pp)  is  characterized  in 
terms  of  the  RC(A,t)  in  two  mutually  exclusive  cases: 

(a)  All  activities  in  the  RC(A,t)  are  processed  simultaneously  in  at  least  one  period. 

(b)  Not  all  the  activities  in  the  RC(A,t)  are  processed  simultaneously  in  any  period. 
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Proof:  Obvious. 


Figure  4.14  shows  a partial  schedule  for  (Pp) 

Period 

|-  1 -|-  2 -]-  3 -j-  4 -j-  5 -|-  6 -j-  7 -|-  8 -|-  9 -|-  10- | 

| <====!====> | | <===========4=========> | 

| <====2 ================> I 1 <====5================> | 

j | <======3========>| <====6====>| 

j |< — 7 >|  |< — 10— — >| 

| | <======8========> | | <=====9=========>| 

Activity  1:  Start  time  = 0,  Duration  = 2,  Finish  Time  = 2 

Activity  3:  Start  time  = 2,  Duration  = 3,  Finish  Time  = 5 

Figure  4.14  A Partial  Schedule  of  a Feasible  Solution  in  Gannt  Chart 

In  Figure  4.14,  activities  4,  5,  6 and  7 are  simultaneously  processed  in  at  least 
one  period,  and  activities  3,  6,  and  7 are  never  processed  simultaneously  in  any  period. 
The  following  two  lemmas  show  the  condition  for  each  case  of  Characterization  4.20. 
Lemma  4.21  summarizes  the  necessary  condition  for  activities  in  the  RC(A,t)  to  be 

processed  simultaneously  in  at  least  one  period,  and  Lemma  4.22  summarizes  the 

necessary  and  sufficient  condition  for  all  the  activities  in  the  RC(A,t)  not  to  be  processed 
simultaneously  in  any  period. 
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Lemma  4.21  The  Necessary  Condition  for  Simultaneous  Processing 

Assume  that  RC(A,t)  is  defined  in  Characterization  4.20.  Suppose  that  in  a 
feasible  solution  of  (Pp)  all  the  activities  in  the  RC(A,t)  are  processed  simultaneously  for 
at  least  one  period.  Then,  there  must  exist  some  mode  m**  G M(i)p  for  each  activity 
in  the  RC(A,t)  such  that  V r „ s b,  must  hold  for  all  k. 

ieRC(A,t) 

Proof:  Assume  that  RC(A,t)  is  defined  in  Characterization  4.20.  Suppose  that  there  is 
some  period  t*  in  some  feasible  solution  of  (Pp)  such  that  all  the  activities  in  the  RC(A,t) 
are  in  progress.  This  implies  that  the  RC(A,t)  is  a subset  of  St»,  which  is  a set  of 
activities  in  progress  at  period  t*  in  the  feasible  solution  of  (Pp).  By  contradiction, 
assume  that  there  is  no  mode  m**  G M(i)p  for  some  activity  in  the  RC(A,t)  such  that 

Er.  . ^ b.  holds  for  all  k.  Since  the  RC(A,t)  is  a subset  of  St„,  constraint 

i(m**)k  k 

ieRC(A,t) 

(13)  of  Definition  4. 15  implies  that  constraint  (3)  of  Definition  4. 15  cannot  be  satisfied. 
Then,  the  schedule  is  not  feasible.  Contradiction.  This  proves  Lemma  4.21.  Q.E.D. 

Lemma  4.22  The  Necessary  Condition  and  Sufficient  Condition  for  Non  Simultaneous 
Processing 

Assume  that  RC(A,t)  is  defined  in  Characterization  4.20.  Suppose  that  in  a 
feasible  solution  of  (Pp),  activities  in  the  RC(A,t)  are  not  processed  simultaneously  at  any 
period.  Then,  there  must  exist  at  least  one  pair  of  activities,  (i,j),  in  RC(A,t)  such  that 
the  start  time  of  activity  j is  no  earlier  than  the  finish  time  of  activity  i. 
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Proof:  Obvious. 

Characterization  4.20  shows  that  any  feasible  solution  of  (Pp)  can  be  characterized 
with  the  RC(A,t)  in  two  mutually  exclusive  cases.  Thus,  two  branching  procedures  are 
presented.  Each  of  the  procedures  deals  with  one  of  the  two  cases  of  Characterization 
4.20.  For  the  convenience  of  notation,  we  will  denote  RC(A,t)  as  RC. 


Procedure  4.23  Branching  Procedure  for  the  Simultaneous  Processing  Case 
Step  1.  Given  an  RC  and  M(i)p  for  each  i E RC,  find  all  the  permutations  of  the 

modes  of  the  activities  in  the  RC.  The  number  of  permutations  is 

n i Msr  i ■ 

ieRC 

Step  2.  For  each  permutation, 

2.1  Find  the  necessary  resource  requirements  when  all  the  activities  in  the  RC 
are  processed  simultaneously  with  the  modes  specified  by  the  permutation. 
Let  NEC(k)  be  the  amount  of  resource  type  k required  for  simultaneous 

processing  for  the  activities  in  the  RC  with  modes  specified  by  the 
permutation.  Then,  NEC(k)  = ri(m*)k  ’ w^ere  m*  *s  determined  by 

ieRC 

permutation. 

2.2  If  NEC(k)  < bk  for  all  k E {1,..,K},  then  create  a child  node  under  the 
current  node  p.  Assume  that  the  newly  created  node  is  called  node  q. 
Update  M(i)q  as  {m*},  where  m*  for  each  i E RC  is  determined  by  the 


permutation. 

If  NEC(k)  > bk  for  any  k E {1,..,K},  then  do  nothing. 


72 


The  existence  of  a mode  combination  that  satisfies  the  condition  of  Lemma  4.21 
is  tested  in  two  steps.  In  Step  1,  all  the  possible  mode  permutations  are  found.  In  Step 
2,  each  permutation  found  in  Step  1 is  tested  whether  the  amount  of  resource  required 
by  the  permutation  does  not  exceed  the  amount  of  available  resource  for  each  resource 
type.  If  NEC(k)  < bk  for  all  k,  then  the  necessary  condition  of  Lemma  4.21  is  satisfied. 
Thus,  we  create  a child  node  under  the  current  node  p with  a restriction  that  each  activity 
in  the  RC  must  be  processed  with  the  mode  specified  from  the  mode  permutation.  An 
example  of  this  procedure  is  presented. 

Example  4.24  An  Example  for  Procedure  4.23 

Suppose  that  RC  = {2,3,8},  and  all  the  relevant  information  is  as  follows: 
Activity  M(i)p 

2 1,2,3 

3 1,2 

8 3 

Resource  Requirement  for  each  Activity(Mode) 


ri(m)k  Resource  Type 

1 

2 

3 

4 

r2(l)k 

1 

2 

1 

1 

r2(2)k 

0 

1 

2 

2 

r2(3)k 

1 

1 

1 

1 

f3(l)k 

2 

1 

1 

1 

r3(2)k 

3 

1 

2 

2 

r8(3)k 

1 

2 

0 

0 

Available  Resource  Level: 

4 

4 

3 

3 
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Step  1.  Permutation  Activity(Mode) 


1 2(1), 3(1), 8(3) 

2 2(1), 3(2), 8(3) 

3 2(2), 3(1), 8(3) 

4 2(2), 3(2), 8(3) 

5 2(3), 3(1), 8(3) 

6 2(3), 3(2), 8(3) 


The  number  of  permutations  is  3*2*1  = 6. 

Step  2.  Permutation  NEC(k)  1 2 3 4 NEC(k)  < bk 

for  all  k 


1 

2 

3 

4 

5 

6 


4 5 2 2 

5 5 3 3 

3 4 3 3 

4 4 4 4 

4 4 2 2 

5 4 3 3 


No 

No 

Yes 

No 

Yes 

No 


In  this  example,  two  children  nodes  are  created  from  the  current  node.  Suppose 
that  the  first  created  node  is  called  node  q.  Then,  M(2)q  = {2},  M(3)q  = {1}  and  M(8)q 
= {3}. 


Procedure  4.25  Branching  Procedure  for  the  Non-Simultaneous  Processing  Case 
Step  1.  Find  all  the  possible  ordered-pairs  of  (different)  activities  in  the  RC.  The 


number  of  pairs  is  | RC  | P2- 
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Step  2.  For  each  ordered  pair  (i,j),  create  a child  node.  Assume  that  the  current 

node  and  the  newly  created  node  are  called  node  p and  node  q, 
respectively.  Update  Hq  as  Hp  U {(i,j)}. 

This  Procedure  is  quite  obvious.  An  example  for  this  procedure  is  presented. 

Example  4.26  An  Example  for  Procedure  4.25 

Suppose  that  RC  = {2,3,8}.  Then,  Step  1 will  produce  the  following  result. 

Ordered  Pairs 

1 (2,3) 

2 (2,8) 

3 (3,2) 

4 (3,8) 

5 (8,2) 

6 (8,3) 

In  Step  2,  six  new  children  nodes  are  branched  from  the  current  node.  Assume 
that  the  first  child  node  branched  by  this  procedure  is  node  q and  that  the  current  node 
is  p.  Then,  Hq  = Hp  U {(2,3)}. 

4.2,3  A Formal  Statement  of  the  Exact  Solution  Procedure 

In  this  section,  we  present  a formal  statement  of  the  exact  solution  procedure. 


All  the  notations  which  are  not  explicitly  defined  in  this  section  follow  the  definitions 
presented  in  the  earlier  sections. 
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Procedure  4.27  Branch  and  Bound  Procedure 
Step  0.  Initialization 

Let  UB  = + oo , where  UB  denotes  the  upper  bound  for  the  problem’s 
optimal  solution  value.  Let  SOLUTION  = {},  where  SOLUTION 
denotes  an  optimal  solution  for  (P°).  Set  p = 0,  where  p denotes  the 
current  node.  Let  P-Node(O)  = -1,  C-Node(O)  = {}  and  k = p,  where 
k denotes  the  last  node  created.  Let  M(i)°  be  M(i)  for  all  i.  Let  H°  be  H. 
Let  LB(Zi)°  be  0 and  UB(zi)°  be  + oo  for  each  activity  i 6 {1,  ..,  I}.  Go 
to  Step  1. 

Step  1.  The  Next  Node  Problem  to  Solve 

We  use  a depth  first  search. 

If  the  current  node  is  node  -1,  then  the  problem  is  solved. 

If  p 5^  -1  and  if  node  p is  not  a dangling  node  and 

if  node  p does  not  have  a dangling  child  node  but  has  children 
nodes,  then  V(PP)  = min  V(PC)  for  all  c G C-Node(p). 
if  node  p does  not  have  any  child  node,  then  V(PP)  = V(LPP). 

if  V(PP)  = V(LPl),  where  t = P-Node(p),  then  fathom  all 
nodes  q G C-Node(t)  and  q ^ p. 

Let  p — P-Node(p). 

Repeat  Step  1. 
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If  p ^ -1  and  if  node  p is  not  a dangling  node  and  if  there  are  some 

dangling  nodes  in  C-Node(p),  then  let  p be  the  node  with  the 
smallest  index  node  among  the  dangling  nodes  in  C-Node(p). 
Repeat  Step  1. 

If  node  p is  a dangling  node,  then  the  next  node  problem  to  solve  is  (Pp). 
Go  to  Step  2. 

If  p = -1,  then  we  have  solved  the  problem.  This  occurs  when  we  know 
V(P°)  and  the  optimal  solution  for  (P°).  This  optimal  solution 
value  for  the  problem  is  UB  and  the  solution  is  SOLUTION. 

Step  2.  Solving  the  Node  Problem  at  Node  p. 

Construct  (Pp)  and  (LPP)  as  described  in  Definition  4.15  and  Definition 
4.18.  Solve  (LPp). 

If  V(LPP)  > UB,  then  node  p is  fathomed.  Go  to  Step  1. 

If  V(LPP)  < UB,  then  go  to  Step  3. 

Step  3.  Testing  whether  we  can  derive  an  optimal  solution  of  (Pp)  from  the 

solution  of  (LPP).  If  we  can,  then  we  have  solved  (Pp).  If  we  cannot, 
branch  from  node  p. 

3.1  Testing  for  the  Noninteger  Solution 

From  i = 1,  test  whether  z*  or  Fj*  is  an  integer,  where  z*  and  F,*  are 
the  solution  value  for  activity  i from  the  optimal  solution  of  (LPP).  If  z * 
is  not  an  integer,  then  create  children  nodes. 
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i.  MAKECHILDNODE.  MAKE  CHILD  NODE  is  implemented 
as  follows:  Let  k = k +1,  C-Node(p)  = C-Node(p)  U {k},  P- 
Node(k)  = p. 

Let  UB(Zi)k  be  [z*],  where  [x]  denotes  the  largest  integer  greater 
than  or  equal  to  x. 

ii.  MAKE_CHILD_NODE.  Let  LB(Zi)k  be  [7*]  + 1 . 

Go  to  Step  1. 

If  F*  is  not  an  integer,  then  create  children  nodes. 

i.  MAKE  CHILD  NODE.  Let  UB(Fi)k  be  [F  *]. 

ii.  MAKE  CHILD  NODE.  Let  LB(Fj)k  be  [Fj*]  + 1 . 

Go  to  Step  1. 

If  all  the  activities  are  checked  and  if  there  is  no  noninteger  solution,  then 
go  to  Step  3.2. 

3.2  Testing  for  the  existence  of  an  Underestimation  Gap 

From  i = 1 , find  CC/^Zj*)  and  Ct(m)(p>(z*)  for  each  m £ M(i)p.  If  there 
is  no  m*  E M(i)p  such  that  CCi(p)(zi*)  = then  create  children 

nodes. 

i.  MAKE  CHILD  NODE.  Let  UB(Zi)k  be  z*. 

ii.  MAKE  CHILD  NODE.  Let  LB(Zi)k  be  z*.  Go  to  Step  1. 

If  all  the  activities  are  checked  and  if  there  is  no  underestimation  gap,  go 
to  Step  3.3. 
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3.3  Resource  Constraint  Checking 

If  there  are  multiple  modes  for  some  activity  i such  that  CC/^Zj*)  = 
Ci(m)(p)(Zi*),  then  choose  one  among  them.  With  the  start  time  and  the 
duration  found  from  the  solution  of  (LPP)  and  with  the  mode  found  in  step 
3.2,  check  whether  the  current  schedule  satisfies  all  the  resource 
constraints.  If  there  is  no  resource  conflict  at  any  period,  then  we  have 
solved  (Pp).  V(PP)  is  V(LPP).  Update  UB  as  V(PP)  and  SOLUTION  as  the 
solution  of  the  (Pp).  Go  to  Step  1.  If  there  is  a resource  conflict  during 
some  period,  find  the  earliest  period  when  a resource  conflict  occurs. 
Create  an  RC,  the  set  of  all  the  activities  in  progress  during  this  period, 
and  create  children  nodes. 

For  the  Simultaneous  Processing  Case: 

Find  all  the  mode  permutations  as  shown  in  Step  1 in  Procedure 
4.23.  For  each  permutation,  test  whether  NEC(k)  < bk  for  all  k 
€ {1,..,K}  as  shown  in  Step  2 in  Procedure  4.23.  If  so,  then 
MAKECHILDNODE.  Update  M(i)k  = m*  for  all  i 6 RC, 
where  m*  is  a mode  chosen  in  the  permutation. 

For  the  Nonsimultaneous  Processing  Case: 

Find  all  the  possible  ordered-pairs  of  activities  in  the  RC  as  shown 
in  Step  1 of  Procedure  4.25.  For  each  ordered-pair,  (i,j), 
MAKE  CHILD  NODE.  Update  Hk  = Hp  U (i,j). 


Go  to  Step  1. 
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4.3  Refinement  of  the  Procedure 

We  presented  our  exact  solution  procedure,  which  is  a branch  and  bound 
algorithm.  Generally,  the  computational  effort  for  solving  the  problem  with  a branch  and 
bound  technique  depends  on  the  number  of  nodes  to  be  solved  and  the  computational 
effort  for  solving  each  node  problem.  Thus,  the  computational  effort  is  determined  by 
the  branching  rule  (how  to  create  nodes),  the  pruning  rule  (how  to  fathom  nodes)  and  the 
solution  procedure  for  the  node  problem.  The  number  of  possible  feasible  schedules  is 
enormous  and  increases  exponentially  as  the  number  of  activities  in  the  project  and  the 
number  of  modes  for  each  activity  increase.  Thus,  a branching  rule,  which  partitions  the 
feasible  region,  is  quite  important.  In  this  section,  we  will  present  some  pruning  rules, 
which  identify  a feasible  region  that  cannot  contain  an  optimal  solution  or  an  improved 
solution,  and  we  will  introduce  a refined  branching  rule. 

4,3.1  Minimal  Resource  Conflict  Set 

Suppose  that  there  exists  an  RC  in  some  node  p.  Then,  the  number  of  children 
nodes  to  be  branched  from  node  p is  at  least  | RC  | P2  since  we  create  a child  node  per 
ordered  pair  of  activities.  If  there  exists  some  mode  permutation  that  satisfies  the 
necessary  condition  for  the  simultaneous  processing  case,  the  total  number  will  be  more 
than  | RC  | P2.  For  example,  if  the  size  of  an  RC  is  2,  then  the  number  of  children  nodes 
is  at  least  2.  If  the  size  of  the  RC  is  4,  then  the  number  of  children  nodes  is  at  least  12. 
If  the  size  of  the  RC  is  6,  then  the  number  of  children  nodes  is  at  least  30.  This  rapid 
growth  of  the  number  of  children  nodes  as  the  size  of  the  RC  grows  leads  us  to  consider 
a more  efficient  branching  rule. 
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Suppose  that  an  RC  in  some  node  p is  {i,j,  ...}  and  that  activities  i and  j cannot 
be  processed  simultaneously  with  the  modes  found  in  Step  3.3  of  Procedure  4.27.  Then, 
any  feasible  solution  for  (Pp)  can  be  characterized  as  Fj  + z,  < Fj  or  Fj  + zi  < Fj  for 
the  non-simultaneous  processing  case,  or  as  xi(mt)  = 1 and  xj(m+)  = 1 such  that  ri(m,)k  + 
rj(m.)k  < bk,  for  all  k 6 K,  for  the  simultaneous  processing  case.  Thus,  branching  with 
only  two  activities  can  result  in  a large  reduction  in  the  number  of  nodes  created  under 
node  p.  This  branching  concept,  which  uses  a subset  of  the  RC,  is  generalized  in  the 
following  definition. 

Definition  4.28.  Minimal  Resource  Conflict  Set 

Suppose  that  there  exists  some  nonempty  RC  in  some  node  p.  A minimal 
resource  conflict  set  is  defined  as  the  smallest  subset  of  the  RC  that  still  causes  resource 
constraint  violation  with  the  modes  assigned  from  Step  3.3  in  Procedure  4.27. 

It  is  easy  to  see  that  branching  with  a minimal  resource  conflict  set  is  a valid 
branching  rule.  The  rule  for  finding  a minimal  resource  conflict  set  is  summarized  in 
the  following  procedure. 

Procedure  4.29  Procedure  for  Finding  a Minimal  RC 

Step  1.  Let  p = | RC  | , q = 1.  Suppose  that  m*  is  a mode  assigned  to 

activity  i in  Step  3.3  in  Procedure  4.27. 

Let  q = q + 1.  If  q = p,  then  stop.  The  RC  itself  is  a minimal  RC. 
If  q < p,  then  go  to  Step  3. 


Step  2. 
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Step  3.  Find  all  subsets  of  the  RC  that  have  a set  size  that  is  equal  to  q.  The 

number  of  subsets  is  nCa.  For  each  set,  check  whether  V r „ > 

ieSUBSET 

bk  for  any  k £ If  there  exists  a constraint  violation,  then  we 

have  found  a minimal  resource  conflict  set.  If  there  is  no  such  set,  go  to 
Step  2. 

4,3.2  Node  Fathoming  Rules 

Fathoming  a node  occurs  when  the  feasible  region  of  a node  is  found  not  to 
contain  any  improved  feasible  solution.  The  Procedure  4.27  contains  two  fathoming 
rules  as  described  in  Steps  1 and  2.  The  efficiency  of  a branch  and  bound  procedure 
partially  depends  on  the  efficiency  of  fathoming  rules  since  fathoming  nodes  reduces  the 
number  of  node  problems  to  be  solved. 

Suppose  that  V(Ppl)  is  known,  or  suppose  that  node  pi  is  fathomed.  Further 
suppose  that  F(Pp2)  is  a subset  of  F(Ppl).  This  implies  that  if  the  objective  functions  of 
(Ppl)  and  (Pp2)  are  identical,  then  V(Pp2)  > V(Ppl)  holds.  Notice  that  the  objective 
functions  of  all  the  node  problems  in  our  branch  and  bound  tree  are  identical.  Let  us 
consider  the  following  two  examples  in  Figure  4.15. 

The  notations  in  the  Figure  4.15  are  as  follows: 

Np:  node  p 

{i,j}:  RC  = {i,j} 

i(m),  j(n):  M(i)q  = {m},  M(j)q  = {n},  where  q is  the  child  node  index. 

i-»j:  Fj  + Zj  < Fj  is  imposed  to  node  q,  where  q is  the  child  node  index. 
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Figure  4.15  A Sample  Branch  and  Bound  Tree 

Since  our  exact  solution  procedure  is  based  on  the  depth-first  search,  the  relaxed 
problems  are  solved  in  the  following  sequence:  node  0,  1,  2,  3,  5,  8,  9,  10,  6,  7,  4,  11, 
12  and  13.  Suppose  the  current  node  is  node  8.  This  means  that  (P1)  and  (P2)  have 
already  been  solved  or  at  least  verified  that  they  cannot  yield  a better  solution  than  some 
known  feasible  solution.  In  other  words,  the  best  known  feasible  solution  must  be  better 
than  V(P1)  when  the  algorithm  reaches  node  8.  All  the  additionally  imposed  constraints 
on  node  1 are  M(l)  = {1}  and  M(2)  = {2},  while  all  the  additionally  imposed 
constraints  on  node  8 are  Fj  + zx  < F2,  M(l)  = {1},  M(3}  = {2},  M(2)  = {2}  and 
M(4)  = {1}.  Thus,  F(P8)  is  a subset  of  FCP1),  and  this  implies  that  node  8 cannot  have 
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a better  solution  than  node  1.  Therefore,  (P8)  does  not  have  to  be  considered.  Let  us 
define  (q)-Constraint-Set  as  the  constraint  set  that  is  imposed  on  node  q when  node  q is 
branched  from  its  parent  node  p such  that  F(Pq)  = F(PP)  n (q)-Constraint-Set.  Now, 
we  will  consider  the  following  case. 


Figure  4.16  A Sample  Branch  and  Bound  Tree 

Suppose  that  node  problems  1,  2 and  3 have  been  and/or  have  been  fathomed. 
It  is  easy  to  see  that  F(P7)  £ F(P2).  Further,  F(P9)  F(P3).  Let  us  consider  any 

feasible  solution  of  (P9).  Then,  the  solution  must  satisfy  all  the  constraints  (P°)  and  also 
F2  + z2  < F3  ((3)-Constraint-Set)  and  F3  + z3  < F,  ((9)-Constraint-Set).  The  two 
additional  constraint  sets  imply  that  F2  + z2  < F3  < F3  + z3  < F,  since  z3  > 0.  Thus, 
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any  feasible  solution  of  (P9)  is  also  a feasible  solution  of  (P3).  Therefore,  node  7 and 
node  9 can  be  fathomed  before  solving  them. 

The  following  Remark  summarizes  the  conditions  for  the  subset  fathoming  rule. 

Remark  4.30  The  Conditions  for  the  Subset  Fathoming  Rule 

Suppose  that  the  solution  of  node  A is  found,  or  suppose  that  node  A is  fathomed. 
If  node  B has  the  same  parent  node  as  node  A and  if  some  node  branched  from  node  B 
includes  the  (A)-Constraint-Set,  then  node  B can  be  fathomed. 

Proof:  Obvious. 

It  is  worthwhile  to  note  that  the  subset  fathoming  rule  for  some  node  applies  only 
when  the  node  is  branched  due  to  a resource  conflict.  If  two  children  nodes  are  created 
due  to  an  underestimation  gap  of  activity  i at  z*  at  node  p,  then  the  feasible  regions  of 
the  children  nodes  differ  from  each  other  because  z{  < z * for  the  first  child  node  and 
z,  > z*  for  the  second  child  node.  It  is  easy  to  see  that  there  cannot  be  an 
underestimation  gap  at  zx  = UB(Zj)p  or  zf  = LB(Zi)p.  Thus,  the  feasible  region  of  the  first 
child  node  cannot  be  a subset  of  that  of  the  second  child  node,  and  vice  versa.  A simple 
procedure  for  the  subset  fathoming  rule  is  presented. 

Procedure  4.31  Procedure  for  the  Subset  Fathoming  Rule 

When  the  solution  for  node  p is  found  or  node  p is  fathomed: 

If  a (p)-Constraint-Set  consists  of  a mode  permutation  or  a precedence  constraint, 
then  include  the  (p)-Constraint-Set  in  a Fathomed-Set  of  each  non  fathomed  node  that  has 
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the  same  parent  node  as  node  p.  A Fathomed-Set  of  some  node  k consists  of  an 
additional  constraint  set  of  an  already  fathomed  node  that  has  the  same  parent  node  as 
node  k. 

When  we  test  whether  node  q is  a subset  of  an  already  solved  node  or  a fathomed 

node: 


Step  1. 

If  the  (q)-Constraint-Set  consists  of  a mode  permutation,  then  go  to  Step 
2.  Let  p be  P-Node(q).  If  the  (q)-Constraint-Set  consists  of  a precedence 
constraint,  then  go  to  Step  3.  Otherwise,  stop.  Node  q is  not  a subset  of 
any  solved  or  fathomed  node. 

Step  2. 

If  p =0,  then  stop.  Node  q is  not  a subset  of  any  solved  or  fathomed 
node.  If  p ^ 0,  then  check  the  Fathomed-Set  of  node  p.  If  the 
Fathomed-Set  in  node  p contains  mode  permutations,  then  check  each 
mode  permutation.  If  M(i)  in  the  permutation  equals  M(i)q  for  all  i in  the 
permutation,  fathom  node  q and  stop.  If  there  is  no  such  permutation  in 
the  Fathomed-Set  of  node  p,  then  let  p = P-Node(p),  and  repeat  Step  2. 

Step  3. 

Let  p be  P-Node(q). 

3.1 

Update  AS(i)  for  all  i,  where  AS(i)q  is  defined  as  a set  of  activities  that 
must  start  after  the  completion  of  activity  i at  node  q. 

3.2. 

If  p = 0,  then  stop.  Node  q is  not  a subset  of  any  solved  or  fathomed 
node.  If  p ^ 0,  then  check  the  Fathomed-Set  of  node  p.  If  the 
Fathomed-Set  in  node  p contains  additional  precedence  constraints,  check 
each  additional  precedence  constraint.  Suppose  that  the  additional 
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precedence  constraint  in  consideration  is  Fk  + zk  < Fm.  Then,  check 
whether  m £E  AS(k)q.  If  true,  fathom  node  q and  stop.  If  there  is  no 
such  additional  precedence  constraint,  then  let  p = P-Node(p)  and  repeat 
Step  3.2. 


CHAPTER  5 

A HEURISTIC  PROCEDURE 


Although  the  exact  solution  procedure  presented  in  Chapter  4 can  find  an  optimal 
solution  to  the  problem  defined  in  Chapter  3,  the  required  computational  effort  grows 
exponentially  with  the  size  of  the  problem.  This  led  us  to  devise  a heuristic  procedure 
that  finds  good  quality  solutions  with  modest  computational  effort.  This  chapter  is 
devoted  to  the  presentation  of  our  heuristic  procedure.  We  present  four  types  of 
improvements  of  a feasible  schedule  in  section  5.1.  In  section  5.2,  we  explain  how  to 
generate  a feasible  schedule.  A formal  statement  for  our  heuristic  procedure  is  presented 
in  section  5.3. 


5.1  Four  Types  of  Improvements 

Once  a feasible  schedule  is  constructed,  four  types  of  improvements  are  used  in 
the  heuristic.  For  ease  of  presentation,  the  following  definitions  and  notations  are 
needed.  All  the  definitions  and  the  notations  that  are  not  explicitly  defined  in  this  chapter 
follow  the  definitions  and  the  notations  defined  in  the  previous  chapters. 

Definition  5.1  Schedule,  Element  of  Schedule,  Feasible  Schedule  and  Solution 

A schedule  is  an  (1+2)-  tuple  of  quadruples  {(0,  1,  ...,  1+1),  (F0,  ...,  F1+1),  (zq, 
..,  zI+1),  (MODE[0],  ...,  MODE[I+l])},  where  MODE[i]  is  an  assigned  mode  for 


87 


88 


activity  i G {0,..,I  + 1}  and  MODE[0]  and  M0DE[I  + 1]  are  defined  as  0.  F„  z,  and 

MODE[i]  are  called  the  elements  of  a schedule.  Suppose  that  the  following  conditions 

are  satisfied  for  a given  schedule. 

(i)  All  the  precedence  and  the  resource  constraints  are  satisfied. 

(ii)  y i(MODE[i])c  — zi  — yi(MODE[i])n  an£^  MODE[i]  G M(i)  for  all  i G { 1 , - . , 1} 

(iii)  Fj  and  z,  are  integers  for  all  i G {0,..,I+1}. 

(iv)  F0  = Zq  = zI+1  = 0. 

(v)  Solution  is  finite,  where  Solution  is  defined  as 

i 

E Ci(MODE[i])(zi)  + Penalty  *max{0,  FI+1  - DUE_DATE}  • 

i=l 

Then,  the  schedule  is  called  a feasible  schedule. 

It  is  easy  to  see  that  a feasible  solution  to  (P)  can  be  derived  from  a feasible 
schedule. 


Notation  5.2 
FINISH^ 
COST*: 
SLACK[p][k]: 


Additional  Notations  in  the  Heuristic 

completion  time  of  activity  i in  the  schedule 

activity  cost  of  activity  i in  the  schedule.  COSTj  = C^ode^^. 

the  amount  of  the  unused  renewable  resource  type  k during  period 


p.  SLACK[p][k]  = bk  - £ ,mm  . 

ieSn 

P(i): 

S(i): 


a set  of  immediate  predecessors  of  activity  i 
a set  of  immediate  successors  of  activity  i 
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ES,: 

in  a given  feasible  schedule,  the  earliest  time  at  which  activity  i 
can  be  rescheduled  to  start  without  violating  any  precedence 
relations  while  assuming  that  none  of  the  all  the  elements  in  the 
schedule,  except  F;  and  zh  change.  ES(  = max  {Fj  + Zj  | j E 

P(i)}. 

ES^: 

in  a given  feasible  schedule,  ES(  based  on  P(i)\{k}  not  on  P(i). 
ESiN{k)  = max  {Fj  + z,  | j E P(i)\{k}}. 

LF;: 

in  a given  feasible  schedule,  the  latest  time  by  which  activity  i can 
be  rescheduled  to  finish  without  violating  any  precedence  relations 
while  assuming  that  none  of  the  elements  in  the  schedule,  except 
Fj  and  zh  change.  LFj  = min  {Fj  | j E S(i)}. 

LFiV{k}: 

in  a given  feasible  schedule,  LFj  based  on  S(i)\{k}  not  on  S(i). 
LFiV{k)  = min  {Fj  | j E S(i)\{k}}. 

CR  COSTji 

the  net  increase  of  COST;  caused  by  reducing  z{  one  period  in  a 
given  feasible  schedule. 

UNSAVINGj: 

CR  COST  - i C‘(MODEW)  **  y i(MODE[i])c  K Zi  ^ y i(MODE[i])n 
- 1 \ °°  if  Zj  = y i(MODE[i])c 

the  net  reduction  of  COST;  caused  by  increasing  z{  one  period  in 
a given  feasible  schedule. 

SUM_CR[t]: 

UN  SAVING  = | C*<M°DE[i])  if  yi(MODE[i])c  ^ Zi  < yi(MODE[i])n 

1 1 - 00  if  Zj  = y i(MODE[i])n 

sum  of  CR  COST;  over  all  i E St 
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SUM_UN[t]:  sum  of  UN  SAVINGj  over  all  i G St  and  UNSAVINGj  > 0 such 

that  SUM  _ UN[t]  = £ UN_SAVINGi 

ieSp  UN_S  AVlNGj  > 0 

ELIG:  a set  of  unscheduled  activities  whose  predecessors  have  been 

completed. 

5.1.1  The  First  Type  of  Improvement 

Suppose  that  we  have  a feasible  schedule  and  some  modifications  are  done  in  the 
feasible  schedule,  where  the  modified  schedule  is  also  feasible.  For  ease  of  presentation, 
we  will  refer  to  the  initial  feasible  schedule  as  the  initial  schedule  and  refer  to  the 
modified  feasible  schedule  as  the  modified  schedule.  In  order  to  distinguish  an  element 
in  the  initial  schedule  and  the  element  in  the  modified  schedule,  we  put  ’ in  the  element 
in  the  modified  schedule.  For  example,  z{  denotes  a duration  of  activity  i in  the  initial 
schedule,  and  zs’  denotes  a duration  of  activity  i in  the  modified  schedule. 

In  the  initial  schedule,  an  attempt  to  reschedule  one  activity  at  a time  is  made  by 
the  first  type  of  improvement.  If  we  consider  rescheduling  activity  i with  the  first  type 
of  improvement,  none  of  the  elements  of  the  initial  schedule,  except  Fh  zK  and  MODEfi], 
are  to  change.  If  activity  i can  be  rescheduled  with  mode  m*  and  with  duration  z*  while 
satisfying  all  the  precedence  and  the  resource  constraints  and  if  COST;  > Ci(m+)(z  *),  then 
the  reschedule  of  activity  i with  mode  m*  and  with  duration  z * will  improve  the  Solution 
by  COSTj  - Cj^/z*)  and  maintain  feasibility.  The  rule  is  summarized  in  the  following 


lemma. 
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Lemma  5.3  The  First  Type  of  Improvement 

Suppose  that  we  a feasible  schedule  and  i E Further  suppose  that  there 

exist  some  periods.  Pi  and  p2,  where  p!  < p2,  and  mode  m*  E M(i)  that  satisfy 
conditions  (i)  to  (iii). 

(i)  (precedence  relation  feasibility)  p!  > (ESi  + 1)  and  p2  < LFj 

(ii)  (resource  feasibility)  For  each  period  pj  < p < p2,  SLACK[p][k]  + ri(MODE[i])k  - 

ri(m*)k  > 0 for  all  k holds  if  Fj  -F 1 < p < Fs  + zi5  and  SLACK[p][k]  - ri(m.)k  > 

0 for  all  k holds,  if  not. 

(iii)  (improved  solution)  COSTj  > Ci(m.)(min  {p2  - Pi  +1,  yi(m*)n}). 

Then,  rescheduling  as  described  in  (iv)  improves  the  Solution  while  maintaining 
feasibility. 

(iv)  (rescheduling  and  updating  cost)  Fj’  = (p,  -1),  z^  = min  |p2  - pj  +1,  yi(m.)n}, 

MODE[i]’  = m*,  COST;’  = Ci(m.)(zi’).  Solution’  = Solution  - COST;  + 

COSTi’.  All  the  other  elements  in  the  modified  schedule  keep  the  values  they  had 
in  the  given  feasible  schedule. 

Proof:  Assume  that  the  assumptions  and  conditions  (i)  to  (iii)  of  Lemma  5.3  hold  and 
that  activity  i has  been  rescheduled  as  shown  in  (iv)  of  Lemma  5.3.  Since  only  activity 
i was  rescheduled,  all  the  precedence  relations  that  do  not  involve  activity  i must  be 
maintained.  Since  px  > (ESj  +1)  and  Fj’  = (p,  -1),  from  the  definition  of  ESi5  F;’  > 
Fj  + Zj  = Fj’  + Zj’  for  all  j E P(i).  Since  p2  < LF|  and  z^  < p2  - Pi  +1,  Fj’  + Zi’  < 
LFj.  Thus,  from  the  definition  of  LFi;  Fj’  + z^  < Fj’  for  all  j E S(i).  Therefore,  the 
precedence  relations  are  maintained  in  the  modified  schedule.  It  is  easy  to  see  that  (ii) 
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of  Lemma  5.3  guarantees  that  the  resource  constraint  feasibility  is  also  maintained.  The 
completion  time  of  the  modified  schedule  is  identical  to  the  completion  time  of  the  initial 
feasible  schedule,  and  all  activities’  start  times  and  durations  in  the  new  schedule  are 
integers.  Since  COSTj,  j ^ i,  is  not  affected  by  the  rescheduling  of  activity  i and  COST; 
> Ci(m*)(min  {p2  - Pi  +1,  yi(m«.)n})  is  assumed,  Solution’  < Solution  follows.  It  is  easy 
to  see  that  other  conditions  of  Definition  5.1  are  met.  This  proves  the  Lemma. 

Q.E.D. 

Since  a mode  cost  at  duration  that  is  less  than  its  crash  duration  or  longer  than 
its  normal  duration  is  defined  as  infinite,  the  duration,  which  is  given  by  the  min  {p2  - 
Pj  +1,  yi(m*)n},  is  used  for  the  comparison  between  COST;  and  C^^min  {p2  - Px  +1, 
yi(m,)n}),  and  z'x  is  updated  to  the  duration  if  conditions  (i)  to  (iii)  in  Lemma  5.3  are 
satisfied. 

5.1.2  The  Second  Type  of  Improvement 

The  second  type  of  improvement  is  related  to  the  duration  adjustments  of  two 
activities  i and  j such  that  FINISH;  = Fj  in  the  current  schedule.  The  next  rule  tests  the 
possibility  of  improving  the  Solution  by  reducing  zx  and  increasing  zy  while  maintaining 
the  current  values  of  F;  and  FINISH^  The  elements  whose  values  can  be  changed  in  this 
case  are  zh  F3  and  zy  To  illustrate  this  rule,  two  partial  schedules  are  provided  in  Figure 


5.1. 
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A Given  Feasible  Schedule,  FINISH2  = F4 
period 


123456789  10  11  12 

! ! ! ! _ _ 


<=  1 => 
<====== 


==  3 =========> 


A Schedule  after  Reducing  z2  by  One  period  and  Increasing  z4 
by  One  Period  while  Maintaining  F2  and  FINISH4  values. 

period 


123456789  10  11  12 

_ _ ! ! ! ! ! ! ! _ _ 


<=  1 => 
<====== 


===>  <=  5 => 


Figure  5.1  The  Second  Type  of  Improvement 

Lemma  5.4  The  Second  Type  of  Improvement 

Suppose  that  we  have  a feasible  schedule.  Assume  that  FINISH*  equals  Fj  in  the 
schedule,  where  i,  j E {1,..,I}.  Suppose  that  CRCOST*  < UN  SAVINGj  and  that 
there  exists  a period  p,  1 < p < FINISH*,  which  satisfies  conditions  (i)  to  (iii)  below. 

(i)  (precedence  feasibility)  p > (ESj.{i}  + 1) 

(ii)  (resource  feasibility)  for  all  p < t < FINISH*  and  for  all  k,  SLACK[t][k]  + 

ri(MODE[i])k  “ rj(MODE[j])k  — 0. 

(iii)  (duration  feasibility)  p -1  - F,  > y^DE^c-  FINISHj  - p + 1 < yj(MODEU])n- 
Then,  the  following  rescheduling  results  in  an  improved  feasible  schedule. 

(iv)  (rescheduling)  z*’  = p -1  - F*,  Fj’  = p -1,  Zj’  = Fj’  - FINISHj. 
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(updating)  COST;  — Ci(MoDE[i])(zi  )>  COSTj  — Cj^ode^^Zj  ).  Solution’  = 
Solution  - COST;  - COSTj  + COST;’  + COSTj’. 

All  the  other  elements  in  the  modified  schedule  keep  the  values  they  had  in  the 
initial  feasible  schedule. 

Proof:  Suppose  that  all  the  assumptions  and  conditions  (i)  to  (iii)  in  Lemma  5.4  hold. 
Since  p < FINISH;  and  the  value  of  F;  does  not  change,  shortening  z{  does  not  cause  any 
violation  of  any  precedence  relations  that  involve  activity  i.  Note  that  (p  -1)  > ESjx{i) 
and  the  value  of  FINISH,  does  not  change.  Also  FINISH;  is  still  Fj  in  the  modified 
schedule.  Therefore,  starting  activity  j at  period  p while  maintaining  the  same  finish 
time  does  not  cause  any  violation  of  any  of  the  precedence  relations  that  involve  activity 
j.  Thus,  the  new  schedule  is  feasible  with  respect  to  the  precedence  relations.  For  all 
periods  from  p to  FINISH;,  activity  i,  which  was  in  progress  in  the  previous  schedule, 
is  no  longer  active  in  the  modified  schedule.  The  reverse  is  true  for  activity  j.  Since 
condition  (ii)  guarantees  that  the  change  of  resource  requirements  in  those  periods  does 
not  cause  any  resource  constraint  violation,  the  new  schedule  is  feasible  with  respect  to 
the  resource  constraints.  Condition  (iii)  of  the  lemma  and  rescheduling  Z;’  = F;  - p -1 
and  Zj’  = Fj’  - FINISHj  ensure  that  Z;  > yi(MoDE[i])c  and  z,  < yj[MODEU])n-  Thus,  the  new 
schedule  is  feasible  with  respect  to  activity  duration.  Notice  that  CR  COST;  < 
UN  SAVINGj  in  the  initial  feasible  schedule.  This  implies  that  Cj(MODE[i])  < cj(MODEUJ). 
Notice  that  the  project  completion  time  has  not  changed.  Therefore,  crashing  activity  i 
and  uncrashing  activity  j by  the  same  number  of  periods  while  maintaining  precedence. 


95 


resource  and  duration  feasibilities  results  in  improvement  in  the  Solution’.  It  is  easy  to 
see  that  other  conditions  of  Definition  5.1  are  met.  This  proves  the  Lemma.  Q.E.D. 

5.1.3  The  Third  Type  of  Improvement 

The  third  type  of  improvement  is  a mirror  image  of  the  second  type  of 
improvement.  Like  the  second  type  of  improvement,  this  type  of  improvement  is  related 
to  the  duration  adjustments  of  two  activities  i and  j,  but  the  conditions  are  reversed.  In 
this  type  of  improvement,  FINISHj  equals  Fj,  activity  i will  be  uncrashed  and  activity  j 
will  be  crashed  by  the  same  number  of  periods  while  the  values  of  F;  and  FINISHj  are 
maintained.  The  rule  for  this  type  of  improvement  is  summarized  in  Lemma  5.5.  To 
illustrate  this  type  of  improvement,  two  partial  schedules  are  provided  in  Figure  5.2. 


A Given  Feasible  Schedule,  FINISH2  = F4 
period 
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A Schedule  after  Increasing  FINISH2  and  F4  by  One  Period 

period 
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Figure  5.2  The  Third  Type  of  Improvement 
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Lemma  5.5  The  Third  Type  of  Improvement 

Suppose  that  we  have  a feasible  schedule.  Assume  that  FINISH;  equals  Fj  in  the 
schedule,  where  i,  j E {1,..,I}.  Suppose  that  UNSAVINGj  > CRCOSTj  and  that 
there  exists  period,  p > FINISHj  +1,  that  satisfies  conditions  (i)  to  (iii)  given  below. 

(i)  (precedence  feasibility)  p < LFix{j} 

(ii)  (resource  feasibility)  for  all  (FINISHj  + 1)  < t < p,  SLACK[t][k]  + ri(MODE[i])k - 
rj(MODE[j])k  — 0- 

(iii)  (duration  feasibility)  p - F,  < yi(MoDE[i])n-  FINISHj  - p > yj(MODEQ])c- 

Then,  the  following  rescheduling  improves  the  schedule  while  maintaining 
feasibility. 

(iv)  (rescheduling)  z,’  = p - F„  Fj’  = p,  Zj’  = F/  - FINISHj. 

(updating)  COST;  = Ci(MODE[i])(zi),  COSTj’  = Cj(M0DE[j])(Zj).  Solution’  = Solution 
- COSTj  - COSTj  + COSTj’  + COSTj’. 

Proof:  Similar  to  the  Proof  of  Lemma  5.4. 

5.1.4  The  Fourth  Type  of  Improvement 

The  fourth  type  of  improvement  attempts  to  crash  or  uncrash  all  the  activities  in 
progress  during  some  period  t.  For  the  fourth  type  of  improvement,  we  need  the 
following  research.  The  next  proposition  applies  to  the  case  when  all  the  activities  in  St, 
t > 1 , are  crashed  by  one  period  and  each  activity  whose  start  time  is  equal  to  or  later 
than  t is  reassigned  to  start  one  period  earlier. 
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Proposition  5.6  Duration  Reduction  of  All  Activities  in  St  for  some  t G {1,..,FI+1} 
Suppose  that  we  have  a feasible  schedule.  Assume  that  t*  is  a period  such  that 
1 < t*  < FI+1  and  S*  is  not  empty.  If  we  reduce  the  duration  of  each  activity  in  St*, 
i 0,  i ^ (I  +1),  by  one  period  while  maintaining  its  start  time,  and  if  we  reassign  F;’ 
= Fj  -1  for  all  i such  that  F;  > t*,  then  the  schedule  is  still  precedence  feasible. 

Proof:  Assume  that  all  the  assumptions  of  Proposition  5.6  are  met.  Suppose  that  t*  is 
an  arbitrary  period  between  1 and  FI+1.  Let  us  categorize  all  the  nondummy  activities 
in  the  project  into  three  cases  with  respect  to  t*:  (a)  activity  whose  finish  time  is  earlier 
than  t*,  (b)  activity  in  St*,  and  (c)  activity  whose  start  time  is  equal  to  or  later  than  t*. 
It  is  easy  to  see  that  all  the  nondummy  activities  must  be  categorized  only  by  one  case. 
Assume  that  a new  schedule  has  been  created  by  using  the  modification  shown  in  the 
Proposition.  Then,  the  followings  hold: 

(i)  Suppose  that  activity  i and  j are  categorized  by  (a)  and  FINISH;  < F,.  Then, 
FINISH;’  < F/  since  FINISH;’  = FINISH;  and  F/  = Fj. 

(ii)  Suppose  that  activity  i is  categorized  by  (a)  and  activity  j is  categorized  by  (b). 
This  implies  that  FINISH;  < Fj.  Then,  FINISH;’  < Fj’  since  FINISH;’  = 
FINISH;  and  Fj’  - Fj. 

(iii)  Suppose  that  activity  i is  categorized  by  (a)  and  activity  j is  categorized  by  (c). 
This  implies  that  FINISH;  < Fj  -1.  Then,  FINISH;’  < Fj’  since  FINISH;’  = 
FINISH;  -1  and  Fj’  = Fj’  -1. 
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(iv)  Suppose  that  activity  i is  categorized  by  (b)  and  activity  j is  categorized  by  (c). 
This  implies  that  FINISH;  < Fj.  Then,  FINISH;’  < Fj’  since  FINISH;’  = 
FINISH;  -1  and  F/  = F/  -1. 

(v)  Suppose  that  activity  i and  j are  categorized  by  (c)  and  F;  < Fj.  Then,  FINISH;’ 
< Fj’  since  FINISH;’  = FINISH;  -1  and  F/  = Fj’  -1. 

Notice  that  activities  that  are  categorized  by  (b)  do  not  have  precedence  relations 
with  other  activities  that  are  also  categorized  by  (b).  Thus,  the  above  cases  (i)  to  (v) 
represent  all  the  precedence  relations  of  nondummy  activities.  Since  all  the  nondummy 
activity  durations  are  assumed  to  be  positive  and  an  integer,  Z;  > 1 for  all  i 6 {1, 

I}.  Thus,  it  is  easy  to  see  that  min  {F;’+Z;’  | i £ {1,  1+1}  > 0}.  This  implies 

that  all  the  precedence  relations  involving  activity  0 hold  in  the  new  schedule.  Also  it 
is  easy  to  see  that  F1+1  ’ = FI+1  -1  and  that  FINISH;’  < F1+1’  for  all  i £ {0,  ..,  I}. 
Thus,  the  new  schedule  is  precedence  feasible.  Q.E.D. 

The  next  proposition  applies  to  the  case  when  all  the  activities  in  St,  t > 1,  are 
uncrashed  by  one  period  and  each  activity  whose  start  time  is  equal  to  or  later  than  t is 
reassigned  to  start  one  period  later. 

Proposition  5.7  Duration  Increase  of  all  Activities  in  St 

Suppose  that  we  have  a feasible  schedule.  Assume  that  t*  is  a period  such  that 
1 < t*  < FI+1  and  St*  is  not  empty.  If  we  increase  the  duration  of  each  activity  in  St», 
i ^ 0,  i ^ (I  +1),  by  one  period  while  maintaining  its  start  time,  and  if  we  reassign  F;’ 
= F;  +1  for  all  i such  that  F;  > t*,  then  the  schedule  is  still  precedence  feasible. 
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Proof:  Similar  to  Proof  of  Proposition  5.6. 

The  first  case  of  the  fourth  type  of  improvement  is  used  when  the  project  due  date 
is  violated.  If  there  exists  some  period  p,  1 < p < FI+1,  such  that  SUM_CR[p]  is  less 
than  Penalty,  it  is  logical  to  reduce  the  duration  of  each  activity  in  Sp  in  order  to  reduce 
the  total  penalty  cost  by  shortening  the  project  completion  time.  The  rule  for  this  case 
is  presented  in  the  following  lemma.  To  illustrate  this  case,  Figure  5.3  is  provided. 


A Given  Feasible  Schedule,  where  DUE_DATE  = 9 
period 
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A Schedule  after  Crashing  all  the  Activities  E S3  by  Two  Periods 
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Figure  5.3  The  First  Case  of  the  Fourth  Type  of  Improvement 


Lemma  5.8  The  First  Case  of  the  Fourth  Type  of  Improvement 

Suppose  that  we  have  a feasible  schedule.  Also  suppose  that  FI+1  > DUE_DATE 
and  there  exists  a period,  p*,  such  that  SUM_CR[p*]  < Penalty  and  1 < p*  < FI+1. 
Then,  the  following  rescheduling  results  in  an  improved  feasible  schedule. 
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(i)  Starting  at  period  (p*  +1),  find  the  first  period  in  which  the  set  of  activities  in 
progress  differs  from  Sp*.  Let  p**  be  the  period  -1.  Then,  the  sets  of  activities 
in  progress  in  all  the  periods  between  p*  and  p**  are  identical. 

(ii)  Find  maxreduction  such  that  maxreduction  = min 

{{min  {zj  - yi(MoDE[i])c  i i € Sp*},  p**  - p*  +1,  FI+1  - DUEDATE}. 

(iii)  Reschedule  z/  = - max  reduction  for  all  i G Sp»,  F/  = F(  - max  reduction  for 

all  i such  that  Fj  > p*.  Update  COST/  = Ci(MODE[i])(Zi’)  for  all  i G Sp*  and 
Solution’  = Solution  - (SUM_CR[p*J  - Penalty)  * max  reduction. 

Proof:  Assume  all  the  assumptions  of  Lemma  5.8  hold.  Assume  that  rescheduling  is 
done  as  shown  in  the  lemma.  Notice  that  resource  feasibility  is  determined  by  the 
activities  in  St,  for  all  t > 1,  and  the  assigned  mode  for  each  activity.  It  is  easy  to  see 
that  St’  = St  for  each  1 < t < p*  -1  and  St’  = St  + max_reduction  for  all  p*  < t < FI+1’. 
Since  there  was  no  mode  change,  this  implies  that  the  new  schedule  is  feasible  with 
respect  to  the  resource  constraints.  Precedence  feasibility  can  be  proven  by  applying 
Proposition  5.6.  Although  one  period  reduction  is  assumed  in  Proposition  5.6,  a 
reduction  of  multiple  periods  can  be  induced  easily  from  the  definitions  of  p*  and  p**. 
Notice  that  max_reduction  < zf  - y((MoDE[i])c  f°r  all  i C Sp*.  Since  z/  = zi  - 
max_reduction,  z/  > yi([MoDE[i])c  for  all  1 ^ Sp*  ln  lhe  initial  schedule.  The  other 
conditions  of  Definition  5.1  are  easily  verified  since  the  initial  schedule  is  a feasible 
schedule  and  only  partial  rescheduling  is  done.  It  is  easy  to  see  that  COSTj  = COST/ 
for  all  i $ Sp*  and  COST/  = COSTj  + ci(MODE[j])  * max_duration  for  all  i G Sp*  in  the 
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initial  schedule.  Since  maxduration  < FI+1  - DUEDATE,  the  reduction  in  total 
penalty  cost  by  shortening  the  project  completion  by  maxjreduction  is  Penalty  * 
maxreduction.  Since  SUMCRfp*]  < Penalty  in  the  initial  schedule,  the  Solution  is 
improved  by  max  reduction  * (Penalty  - SUM_CR[p*]).  This  proves  the  lemma. 

Q.E.D. 


The  second  case  of  the  fourth  type  of  improvement  is  used  when  the  project  due 
date  is  violated.  Suppose  that  there  exists  some  period  p,  1 < p < FI+1,  such  that 
SUM_UN[p]  is  greater  than  Penalty.  Notice  that  SUM_UN[p]  is  the  sum  of  positive 
UN  SAVINGj,  where  i E Sp.  In  this  case,  it  is  logical  to  uncrash  each  activity  that  is 
in  Sp  and  whose  UN  SAVING  is  positive  in  order  to  reduce  the  activity  costs  at  the 
expense  of  a higher  total  penalty  cost  by  delaying  the  project  completion  time.  The  rule 
for  this  case  is  presented  in  the  following  lemma.  To  illustrate  this  case,  Figure  5.4  is 
provided. 


Lemma  5.9  The  Second  Case  of  the  Fourth  Type  of  Improvement 

Suppose  that  we  have  a feasible  schedule.  Also  suppose  that  FI+1  > DUE  DATE 
and  that  there  exists  a period  p such  that  SUM_UN[p]  > Penalty  and  1 < p < FI+1. 
Then,  the  following  rescheduling  results  in  an  improved  feasible  schedule. 

(i)  Find  max_increase  such  that  max_increase  = min  {yi(MoDE[i])n  - A \ i E Sp,  z(  < 


y i(MODE[i])n}  • 


(ii)  Reschedule  z/  = zx  + max_increase  for  all  i E Sp  and  zx  < yi(MODE[i])n-  F/  = Fi 
+ max_increase  for  all  i such  that  Fi  > p.  Update  COST/  = C^oDE^fe’)  for 


102 


A Given  Feasible  Schedule,  where  DUEDATE  = 9 and  z2  = y2(MODE[2])n- 

period 
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A Schedule  after  Uncrashing  all  the  Activities  G S3  by  Two  Periods 

period 
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Figure  5.4  The  Second  Case  of  the  Fourth  Type  of  Improvement 

each  i such  that  i G Sp  and  ^ yi(MODE[i])n-  Solution’  = Solution  - (Penalty  - 
SUM_UN[p])  * max_increase. 

Proof:  The  proof  is  similar  to  Proof  of  Lemma  5.8  except  applying  Proposition  5.7 
instead  of  Proposition  5.6. 

In  this  Lemma,  the  max_increase  is  based  on  only  yi(MODE[i])n  and  Zj  for  all  i G Sp. 
Since  the  due  date  is  violated  and  the  Lemma  delays  the  project  completion  time,  FI+1 
does  not  affect  the  maxincrease.  Notice  that  only  the  activities  in  Sp  whose  durations 
are  not  the  normal  durations  are  uncrashed  and  that  UN  SAVINGi  > 0 implies  z{  < 


yi(MODE[i])n- 
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The  third  case  of  the  fourth  type  of  improvement  is  used  when  the  project 
completion  time  is  earlier  than  the  due  date.  In  this  case,  it  is  logical  to  uncrash 
activities  until  the  project  completion  time  is  delayed  to  the  due  date.  The  rule  for  this 
case  is  presented  in  the  following  lemma. 

Lemma  5.10  The  Third  Case  of  the  Fourth  Type  of  Improvement 

Suppose  that  we  have  a feasible  schedule.  Also  suppose  that  FI+1  < DUE  DATE 
and  that  there  exists  a period  p,  1 < 1 < FI+1,  such  that  SUM_UN[p]  > 0.  Then,  the 
following  rescheduling  results  in  an  improved  feasible  schedule. 

(i)  Find  max_increase  such  that  max_increase  = min  {min  {yi(MODE[i])n  ~ z\  I ' G Sp, 
Zi  < y i(MODE[i])n}  > DUE_DATE  - F[  + 1}. 

(ii)  Reschedule  z{’  = z,  + max_increase  for  all  i such  that  i G Sp  and  Zj  < yi(MoDE[i])n- 
Fi’  = Fj  + max_increase  for  all  i such  that  Fj  > p.  Update  COST;’  = 
Ci(MODE[i])(zi’)  for  all  i such  that  G Sp  and  zx  < yi(MODE[i])n-  Solution’  = Solution  - 
(SUM_UN[pj)  * max_increase. 

Proof:  The  proof  is  similar  to  that  of  Lemma  5.9  except  that  the  total  penalty  cost  in  this 
case  is  0. 

In  this  Lemma,  the  max_increase  is  based  on  not  only  yi(MoDE[i])n  and  zi  for  all  i 
G Sp  but  also  (DUE_DATE  -FI+1).  Since  the  due  date  is  not  violated  and  the  lemma 
delays  the  project  completion  time  up  to  the  due  date,  FI+1  must  be  considered. 

The  last  case  in  the  fourth  type  of  improvement  is  used  when  there  exist  two 
periods,  0 < p,q  < FI+1,  such  that  SUM_UN[p]  > SUM_CR[q].  Then,  it  is  better  to 
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uncrash  the  duration  of  each  activity  that  is  in  Sp  and  whose  UNJSAVING  is  positive  and 
to  crash  the  duration  of  each  activity  in  Sq.  The  rule  for  this  case  is  summarized  in  the 
following  lemma.  In  order  to  illustrate  this  case,  Figure  5.5  is  provided. 


A Given  Feasible  Schedule,  where  z2  = y2(MODE[2]>n- 
period 


I 

I 


A Schedule  after  Uncrashing  each  activity  i G S3  such  that  zt  ^ zj(MODE[i])n  by  Two 
Periods  and  Crashing  all  the  Activities  G S7  by  Two  Periods 

period 


Figure  5.5  The  Fourth  Case  of  the  Fourth  Type  of  Improvement 

Lemma  5.11  The  Fourth  Case  of  the  Fourth  Type  of  Improvement 

Suppose  that  we  have  a feasible  schedule.  Also  assume  that  there  exist  two 
periods,  p and  q*,  such  that  SUM_UN[p]  > SUM_CR[q*],  where  1 < p,q*  < FI+1. 
Then,  the  following  rescheduling  results  in  an  improved  feasible  schedule. 

(i)  Starting  at  period  (q*  +1),  find  the  first  period  in  which  the  set  of  activities  in 
progress  differs  from  Sq*.  Let  q**  be  the  period  -1.  Then,  the  sets  of  activities 
in  progress  in  all  the  periods  between  q*  and  q**  are  identical. 
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(ii)  Find  max_change  such  that  max_change  = 
min  { min  {z;  - yi(MODE[i])o  I > ^ Sq.}, 

q**  - q*  +1, 

m'n  {yi(MODE[i])n  “ Zi  I ' ^ Sp,  Zj  < yi(MODE[l])n} } • 

(iii)  Reschedule  as  follows: 

Z;’  = z;  + max_change  if  i G Sp  and  zt  < yi(MoDE[i])n- 

Zj’  = Zj  - max_change  if  i G Sq* 

if  p < q*,  then  reschedule  as  follows: 

F,’  = Fj  + max  change  if  p < Fj  < q* 
if  q*  < p,  then  reschedule  as  follows: 

Fj’  = Fj  - max  change  if  q*  < Fj  < p 
COST;  = Ci(MODE[i](zi ) if  i G Sp  or  Sq* 

Solution’  = Solution  - (SUM_UN[p]  - SUM_CR[q*])  * max  change. 

Proof:  This  proof  can  be  induced  easily  from  Proof  of  Lemma  5.8  and  Proof  of  Lemma 
5.10. 


In  this  case,  max_change  is  based  on  max_reduction  in  Lemma  5.8  and  on 
max  increase  in  Lemma  5.10.  Neither  the  project  completion  time  nor  the  project  due 
date  affects  the  max_change  since  crashing  and  uncrashing  by  the  same  periods  are 


exercised. 
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5.2  Generation  of  a Feasible  Schedule 

The  four  types  of  improvement  routines  introduced  in  section  5.1  are  based  on 
the  assumption  that  a feasible  schedule  is  available.  In  this  section,  a procedure  for 
generating  a feasible  schedule  will  be  presented.  When  we  construct  an  initial  feasible 
schedule,  we  choose  the  mode  and  the  duration  of  each  nondummy  activity  with  the 
following  rules.  For  each  activity,  choose  the  mode  whose  mode  cost  for  its  normal 
duration  is  the  least  expensive  among  the  available  modes  to  the  activity,  and  set  the 
duration  to  the  normal  duration  of  the  chosen  mode  such  that  NCi(MODE[i])  = minimum 
{NCi(m)  | m E M(i)}  and  z,  = yi(MODE[i])n-  Thus,  we  choose  the  initial  mode  and  the 
duration  of  activity  i so  that  COST;  is  minimized.  Since  we  have  already  chosen  a mode 
and  duration,  COST;  for  all  i is  fixed.  Therefore,  the  only  remaining  factor  in  the  initial 
Solution  is  the  project  penalty  cost.  In  order  to  minimize  the  project  penalty  cost,  a 
smaller  project  completion  time  is  preferred  to  a longer  one.  We  now  have  a project 
completion  time  minimization  problem  with  a single  mode,  fixed  duration  and  multiple 
resources.  There  are  many  good  heuristics  for  this  problem.  The  minimum  latest  finish 
time  rule  is  a heuristic  that  gives  priority  to  the  activity  whose  latest  finish  time  is  the 
minimum  for  assignment  among  the  activities  in  the  eligible  set.  A brief  statement  of 
the  minimum  latest  finish  time  rule  is  presented  in  Procedure  5.12.  The  procedure 
assumes  that  initial  mode  and  initial  duration  are  assigned  to  each  activity.  Since  the 
minimum  latest  finish  time  rule  requires  a project  due  date  that  cannot  be  violated  and 
since  the  due  date  in  our  problem  can  be  violated  with  penalty  cost  per  period,  we  set 
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TEMPDUEDATE  as  the  sum  of  Zj,  i = 1,..,  I,  for  the  purpose  of  calculating  the  latest 
finish  time  of  each  activity. 

Procedure  5.12  Procedure  for  Generating  a Feasible  Schedule 
Step  1.  Set  TEMPDUEDATE  as  the  sum  of  z„  i = 1 I.  Based  on 

TEMPDUEDATE,  perform  the  critical  path  analysis,  and  find  the  latest 
finish  time  of  each  activity.  Set  F0  = 0.  t is  set  to  0.  St+1  = {}. 

Step  2.  Update  ELIG.  If  ELIG  is  empty  and  all  the  activities  are  scheduled,  go 

to  Step  5.  If  ELIG  is  empty,  but  if  some  of  activities  are  not  scheduled, 
go  to  Step  4.  If  ELIG  is  not  empty,  then  choose  the  activity  from  ELIG 
whose  latest  finish  time  is  the  smallest,  and  go  to  Step  3. 

Step  3.  If  the  sum  of  resources  required  by  the  activities  in  St+1  and  the  activity 
under  consideration  does  not  exceed  the  capacity  of  the  resource  for  all 
types,  then  the  activity  is  assigned  to  start  at  the  current  t.  The  activity 
joins  St+1.  If  not,  the  activity  is  not  considered  for  scheduling  at  this 
time. 

Step  4.  Update  t as  the  smallest  completion  time  (Fj+Zj)  of  the  activities  in  St+1. 

St+1  is  updated  so  that  all  the  activities  whose  completion  time  are  less 
than  or  equal  to  t are  removed  from  St+1.  Go  to  Step  2. 

Step  5.  Set  COSTj  = Ci(MODE[i])(zi)  for  i = 1,..,  I.  Set  COST0  = COSTI+1  = 0. 

Set  Solution  = £ Ci(MODE[i])(Zi)  + Penalty  * max  {0,  FI+1  - 

i=l 


DUEDATE}. 
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Although  we  employed  the  minimum  latest  finish  time  rule,  any  good  heuristic 
rule  for  solving  the  minimization  problem  of  the  project  completion  time  with  single 
mode,  single  duration  and  multiple  renewable  resources  can  be  used  for  generating  an 
initial  feasible  schedule. 

5.3  A Formal  Statement  of  the  Heuristic  Procedure 
Our  proposed  heuristic  is  presented  in  the  following  procedure. 

Procedure  5.13  Heuristic  Procedure 

Step  0.  (Initialization)  Set  Num  Of  Trial  = 1.  For  all  i E (1,.., I},  choose  each 

activity’s  mode  and  duration  such  that  NCi(MODE[i])  = minimum  {NCi(m)  | 
m E M(i)}  and  Z;  = yi(MoDE[i])n-  Set  MODE[0]  = MODE[I  + l]  = 0 and 
Zo  — zI+1  = 0.  BestSolution  = infinite.  Best_Schedule  = {}. 

Step  1.  If  Num  Of  Trial  < MaxTrial,  go  to  Step  2.  If  Num  Of  Trial  = 
MaxTrial,  then  stop.  Best  Solution  and  Best  Schedule  is  the  best 
solution  value  and  schedule  found  by  this  heuristic,  respectively. 

Step  2.  (Generation  of  a Feasible  Schedule)  Given  the  mode  and  the  duration 

assignment,  generate  a feasible  schedule  with  the  minimum  latest  finish 
time  rule  as  described  in  procedure  5.12.  Adjust  the  SLACK  arrays 
properly.  Pre  Solution  = - infinite.  Update  M_INDEX[i][m]  as  follows: 
If  m = MODE[i],  then  set  M_INDEX[i][m]  = Solution.  Otherwise,  set 
M_INDEX[i][m]  = infinite.  Go  to  Step  3. 
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Step  3.  (Applying  the  Four  Types  of  Improvements)  If  Pre_Solution  = Solution, 

then  go  to  Step  4.  Otherwise,  set  Pre  Solution  = Solution  and  i = 0. 
Go  to  Step  3.1. 

3.1  (Application  of  the  first  type  of  improvement)  If  i = I,  then  let  i = 0, 
and  go  to  Step  3.2.  Otherwise,  set  i = i + 1.  Find  ES;  and  LF;  in  the 
current  feasible  schedule. 

3.1.1  If  there  are  no  untested  modes  in  M(i),  then  go  to  Step  3.2.  Otherwise, 
choose  an  untested  mode  from  M(i),  and  let  m be  the  mode.  Go  to  Step 
3.1.2. 


3.1.2  Find  p,*  and  p2*  such  that  (p2*  - pt*)  is  the  maximum  among  (p;  - p2)  that 
satisfy  the  conditions  of  Lemma  5.3.  If  there  exist  P;*  and  p2*,  then 
reschedule  as  shown  in  (iv)  of  Lemma  5.3.  Adjust  SLACK  arrays.  Let 
M_INDEX[i][m]  = Solution.  Go  to  Step  3.1.1. 

3.2  (Application  of  the  second  type  of  improvement)  If  i = I,  then  let  i = 0, 
and  go  to  Step  3.3.  Otherwise,  set  i = i +1. 

3.2.1  If  there  is  no  untested  activity  j such  that  FINISH;  equals  Fj,  then  go  to 
Step  3.2.  Otherwise,  choose  an  untested  activity,  j,  such  that  FINISH; 
equals  Fj5  update  ES^,  and  go  to  Step  3.2.2. 

3.2.2  Find  p*  such  that  p*  is  the  minimum  among  p’s  that  satisfy  the  conditions 
of  Lemma  5.4.  If  p*  exists,  then  reschedule  as  shown  in  (iv)  of  Lemma 
5.4.  Adjust  SLACK  arrays,  and  go  to  Step  3.2.1. 
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3.3  (Application  of  the  third  type  of  improvement)  If  i = I,  then  let  i = 0, 
and  go  to  Step  3.4.  Otherwise,  set  i = i +1. 

3.3.1  If  there  is  no  untested  activity  j such  that  FINISH,  equals  Fi5  then  go  to 
Step  3.3.  Otherwise,  choose  an  untested  activity,  j,  such  that  FINISH, 
equals  Fb  update  LF^,  and  go  to  Step  3.3.2. 

3.3.2  Find  p*  such  that  p*  is  the  maximum  among  p’s  that  satisfy  the  conditions 
of  Lemma  5.5.  If  p*  exists,  then  reschedule  as  shown  in  (iv)  of  Lemma 
5.4.  Adjust  SLACK  arrays.  Go  to  Step  3.3.1. 

3.4  (Application  of  the  fourth  type  of  improvement) 

3.4.1  Construct  SUM  CRft]  and  SUM_UN[t]  for  all  t such  that  1 < t < FI+1. 

3.4.2  (Application  of  the  first  case  of  the  fourth  type  of  improvement) 

If  Fi+i  > DUEDATE,  then  find  p*  such  that  SUM_CR[p*]  is  the 
minimum  among  SUM_CR[t]  ,l<t<  FI+1,  recorded. 

If  SUM_CR[p*]  < Penalty,  find  p**  and  max  reduction,  and 

reschedule  as  shown  in  Lemma  5.8.  Adjust  SLACK, 
SUMJUN  and  SUM_CR  arrays.  Repeat  Step  3.4.2. 

If  not,  go  to  Step  3.4.3. 

If  not,  go  to  Step  3.4.3. 

3.4.3  (Application  of  the  second  case  of  the  fourth  type  of  improvement) 

If  FI+1  > DUE  DATE,  then  find  SUM_UN[p]  such  that  SUM_UN[p] 
is  the  maximum  among  SUM_UN[t],  1 < t < FI+1,  recorded. 
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If  SUM_UN[p]  > Penalty,  find  max_increase,  and  reschedule  as 
shown  in  Lemma  5.9.  Adjust  SLACK,  SUM_UN  and 
SUMCR  arrays.  Repeat  Step  3.4.3. 

If  not,  go  to  Step  3.4.4. 

If  not,  go  to  Step  3.4.4. 

3.4.4  (Application  of  the  third  case  of  the  fourth  type  of  improvement) 

If  FI+1  < DUEDATE,  then  find  SUM_UN[p]  such  that  SUM_UN[p] 
is  the  maximum  among  SUM_UN[t],  1 < t < FI+1,  recorded. 

If  SUM_UN[p]  > 0,  find  maxincrease,  and  reschedule  as  shown 
in  Lemma  5.10.  Adjust  SLACK,  SUM  UN  and  SUM  CR 
arrays.  Repeat  Step  3.4.4. 

If  not,  go  to  Step  3.4.5. 

If  not,  go  to  Step  3.4.5. 

3.4.5  (Application  of  the  fourth  case  of  the  fourth  type  of  improvement) 

Find  p such  that  SUM_UN[p]  is  the  maximum  among  SUM_UN[t],  1 < 
t < FI+1,  recorded.  Find  q such  that  SUM_CR[q]  is  the  minimum  among 
SUM_CR[t],  1 < t < FI+1,  recorded. 

If  SUM_UN[p]  > SUM_CR[q],  find  max  change,  and  reschedule  as 

described  in  Lemma  5.11.  Adjust  SLACK,  SUM  UN  and 
SUMCR  arrays.  Repeat  Step  3.4.5. 


If  not,  go  to  Step  3. 
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Step  4.  (Update  of  the  BestSolution  and  the  Best  Schedule) 

If  Solution  < Best  Solution,  update  the  Best_Solution  with  the  Solution 
and  the  Best_Schedule  with  the  current  feasible  schedule. 

Increase  NumOfTrial  by  one. 

(Reinitialization  of  MODE[i]  and  z) 

Choose  the  predetermined  number  of  activities  in  the  following  way:  If 
M_INDEX[i*][m*]  is  the  maximum  in  array  M INDEX  and  activity  i*  is 
not  chosen,  then  choose  i*,  and  set  MODE[i*]  = m*.  Set  zit  as  an 
integer  value  randomly  chosen  between  yi+(m„)n  and  yi*(m*)c.  If  ties  happen 
in  choosing  i*,  choose  one  arbitrarily.  Repeat  until  the  predetermined 
number  of  activities  are  chosen.  Go  to  Step  2. 

Once  a feasible  schedule  is  generated  in  Step  2,  the  proposed  algorithm  applies 
the  four  types  of  improvement  techniques  repeatedly  until  an  improved  schedule  cannot 
be  found  by  those  techniques.  In  order  to  check  whether  there  is  no  improvement  from 
the  application  of  the  four  type  of  techniques,  a control  variable,  Pre_Solution,  is 
employed.  When  a new  feasible  schedule  is  constructed  in  Step  2,  the  Pre_Solution  is 
set  to  -infinity.  Whenever  the  application  of  four  types  of  improvement  techniques  (Step 
3)  are  about  to  be  initiated,  Pre_Solution  value  is  compared  to  Solution.  If  Pre_Solution 
equals  Solution,  this  implies  that  no  further  improvement  was  obtained  in  Step  3 of  the 
previous  application.  Thus,  the  algorithm  returns  to  Step  4 in  order  to  update  the 
Best_Solution  and  the  Best_Schedule  if  necessary  and  to  change  the  initial  modes  and 


durations  of  some  activities.  The  reinitialization  of  some  selected  activities’  initial  mode 
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and  duration  is  necessary  to  avoid  the  same  initial  feasible  schedule.  For  the  selection 
of  activities  for  this  purpose,  a two  dimensional  array,  M_INDEX[i][m],  is  introduced. 
When  a new  initial  feasible  schedule  is  constructed,  this  array  is  reinitialized  by  setting 
M_INDEX[i][m]  = Solution  if  m = MODE[i]  and  M_INDEX[i][m]  = infinity  if  m ^ 
MODE[i].  This  array  is  updated  whenever  MODE[i]  is  updated  in  Step  3 such  that 
M_INDEX[i][MODE[i]]  is  set  to  the  Solution  of  the  current  feasible  schedule.  Suppose 
that  the  algorithm  reaches  Step  4 and  that  M_INDEX[i][m]  is  infinite.  This  implies  that 
mode  m has  not  been  selected  for  activity  i in  this  trial.  Thus,  the  maximum  value  of 
M_INDEX[i][m]  is  selected,  and  MODE[i]  is  set  to  m unless  activity  i is  already  chosen. 
In  order  to  ensure  a different  initial  feasible  schedule  in  the  next  trial,  the  initial  duration 
of  each  selected  activity  is  generated  randomly  between  yi(m)n  and  yi(m)c.  This  process  is 
repeated  until  the  predetermined  number  (Num_Max_Trial)  of  initial  feasible  schedules 
are  generated  and  have  been  tried  to  be  improved  upon  the  four  types  of  improvement 


techniques. 


CHAPTER  6 

COMPUTATIONAL  RESULTS 


In  this  chapter  we  report  some  computational  results  for  the  heuristic  and  the 
exact  solution  methods.  Along  with  our  heuristic  and  exact  solution  methods,  the 
backtracking  procedure  developed  by  Patterson  et  al.  [1989]  and  improved  by  Sprecher 
[1994]  is  tested  as  a heuristic  procedure  and  as  an  exact  solution  procedure.  We  compare 
our  procedures  with  the  backtracking  procedure  since  this  procedure  is  the  only  known 
procedure,  as  far  as  we  know,  that  can  solve  our  problem.  The  necessary  modification 
for  the  backtracking  procedure  to  solve  the  problem  is  described  in  Appendix  A. 

For  the  experiment,  a set  of  data  consisting  of  the  precedence  relations,  each 
mode’s  normal  and  crash  durations,  normal  and  crash  costs,  required  resources  per 
period,  each  resource’s  capacity,  the  project  due  date,  and  the  penalty  cost  were 
generated  by  procedures  that  will  be  presented  in  the  following  sections. 


6.1  Data  Generation 


For  the  generation  of  data,  the  following  notations  are  needed. 


Notation  6. 1 Notations  for  Generation  of  Data 


LB_NC[m](UB_NC[m]): 


the  lower  (upper)  bound  for  NCi(m) 


LB_ND[m](UB_ND[m]): 


the  lower  (upper)  bound  for  yi(m)n 
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LB_CR[m](UB_CR[m]): 

LB_CR_F[m](UB_CR_F[m]): 


CRASH  FACTOR: 


LBRES  [m]  [k]  (UB_RES[m]  [k]) : 
S Comp  (CComp): 


mul_due: 
mul_pen: 
mul  res: 


the  lower  (upper)  bound  for  ci(m) 

the  lower  (upper)  bound  for  the  CRASHF ACTOR 

for  mode  m 

a multiplier:  yi(m)c  is  a product  of  yi(m)n  by 

CRASHF  ACTOR.  If  yi(m)c  is  not  an  integer,  then 

the  fractional  portion  is  dropped. 

the  lower  (upper)  bound  for  ri(m)k 

the  project  completion  time  from  the  first  (second) 

critical  path  analysis 

a multiplier  used  for  determining  project  due  date 
a multiplier  used  for  determining  penalty  cost 
a multiplier  used  for  determining  resource  capacity 


6.1.1  Generation  of  Precedence  Relations 

The  precedence  relations,  which  describe  the  project  network,  were  generated  by 
the  algorithm  given  by  Arno  Sprecher  [ 1 994 j . Since  the  algorithm  is  very  well 
documented  by  Sprecher,  we  do  not  explain  it  here. 

6.1.2  Mode  Normal/Crash  Duration  and  Cost,  Crashing  Cost 

For  all  activities,  i = 1,..,I,  we  randomly  selected  an  integer  between  LB_NC[m] 
and  UB_NC[m]  as  NCi(m).  In  a similar  way,  we  randomly  selected  integers  between 
LB_ND[m]  and  UB_ND[m]  and  between  LB_CR[m]  and  UB_CR[m]  for  yi(m)n  and  ci(m), 
respectively.  In  order  to  determine  the  mode  crash  duration,  yi(m)c,  we  used  a 
CRASHFACTOR  such  that  yi(m)c  = yi(m)n  * CRASHF  ACTOR,  where 


116 


CRASHFACTOR  is  randomly  chosen  between  LB_CR_F[m]  and  UB_CR_F[m].  If 
yi(m)n  * CRASH  FACTOR  is  not  an  integer,  the  fractional  portion  is  dropped.  For  each 
resource  demand,  ri(m)k,  we  randomly  chose  an  integer  between  LB_RES[m][k]  and 
UB  RES[m][k]. 

6.1.3  Due  Date.  Penalty  and  Resource  Capacity 

Once  we  have  all  the  cost,  duration,  crashing  cost  and  resource  demand  data 
along  with  the  precedence  relations,  we  perform  the  critical  path  analysis  twice.  In  the 
first  analysis,  the  mode  and  the  duration  of  each  non-dummy  activity  are  fixed  to  m*  and 
yi(m„)c,  respectively,  such  that  yi(m*)c  = min  {yi(m)c  | m E M(i)}.  In  the  second  analysis, 
the  mode  and  the  duration  of  each  non-dummy  activity  are  fixed  to  m**  and  yi(m**)n, 
respectively,  such  that  NCi(m**)  = min  {NCi(m)  | m E M(i)}.  Thus,  the  first  analysis 
is  done  with  the  shortest  duration  of  each  activity,  and  the  second  analysis  is  done  with 
the  cheapest  activity  cost. 

The  project  completion  time  of  the  first  (second)  analysis  is  saved  as  S_Comp 
(C_Comp).  Due  date  is  determined  by  the  following  formula: 

DUE  DATE  = SComp  * mul  due  + CComp  * (1  - muldue). 

In  order  to  determine  the  capacity  of  each  resource  (bk,  k=  1,..,K)  and  the  penalty 
cost  (Penalty),  the  start  time  of  each  activity  in  both  schedules  is  examined.  In  each  start 
time,  the  sum  of  ci(m+)  and  the  sum  of  ri(m+)k  for  each  k are  recorded,  where  i(m+)  is  an 
assigned  mode  for  activity  i in  the  schedule.  The  maximum  (minimum)  of  the  sum  of 
ci(m+>  is  saved  as  MAX_CRASH  (MIN  CRASH),  and  the  maximum  (minimum)  of  the 
sum  of  ri(m+)k  is  recorded  as  MAX_RES[k]  (MIN_RES[kj).  Also,  for  each  activity  and 
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for  each  mode,  the  largest  ri(m)k,  for  k = 1,..,K,  is  examined  and  saved  as 
MAXIMUM[k],  k = 1,..,  K. 

The  penalty  cost  is  determined  by  the  following  formula: 

Penalty  = MAX  CRASH  * mul_pen  + MIN  CRASH  * (l-mul_pen). 

The  capacity  of  each  resource  is  determined  by  the  following  formula: 

bk  = MAX  RES[k]  * mul_res  + MIN_RES[k]  * (1  - mul  res)  for  k = 1,..,K. 

If  bk  is  less  than  MAXIMUM[k],  then  bk  is  modified  as  MAXIMUM[k], 

6.2  Experimental  Design 


We  followed  the  Ceteris  Paribus  design.  Table  6. 1 shows  the  constant  parameters 
used  in  our  experiment,  and  Tables  6.2  and  6.3  illustrate  the  parameters  used  to  generate 
the  Default  problem  set. 


Table  6. 1 Constant  Parameters  in  the  Experiment 


1 So  1 

i p,«  i 

K 

3 

3 

3 

Table  6.2  Parameter  Set  I used  for  Generation  of  the  Default  Data 


m = 1 

m — 2 

m = 3 

LB  NCfml,  UB  NC[ml 

25,65 

50,90 

25,90 

LB  ND|ml,  UB  ND[ml 

12,20 

9,15 

9,20 

LB  CR[ml,  UB  CR|ml 

3,10 

5,15 

3,15 

CRASH  FACTOR 

.30, .70 

.30, .70 

.30, .70 

LB  RESrmiril,  UB  RESfmlfll 

0,3 

3,6 

0,6 

LB  RESfmU21,  UB  RES[m][21 

3,6 

0,3 

0,6 

LB  RESrmir3],  UB  RES[mlf31 

1,5 

2,7 

0,7 
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Table  6.3  Parameter  Set  II  used  for  Generation  of  the  Default  Data 


mul  due 

mul  pen 

mul  res 

.2 

.1 

.4 

Each  problem  set  consists  of  four  problems.  A number  of  problem  sets  with 
different  parameters  were  generated  as  shown  in  Table  6.4.  In  the  table,  (x)  indicates 
that  x is  the  default  setting. 


Table  6.4  Variable  Parameters  in  the  Experiment 


I 

(10)  12  14  16  20  25 

MAXCRASHING 

2 3 Default  * .5  (Default) 

COMPLEXITY 

1.5  (1.8)  2.1 

1 M(i)  | 

1 (2)  3 

PENALTY 

Default  * .8  (Default)  Default  * 1.2 

DUEDATE 

Default  * .8  (Default)  Default  * 1.2 

CAPACITY 

Default  * .8  (Default)  Default  * 1.2  Default  * 1.4 

In  the  second  row,  MAXCRASHING  indicates  the  difference  between  yi(m)n  and 
yi(m)c.  The  i-th  problem  with  MAX_CRASHING  = 2 differs  from  the  i-th  default 
problem  such  that  yi(m)c  of  the  problem  with  MAX_CRASHING  = 2 is  set  to  yi(m)n  - 2, 
for  all  i(m).  All  the  problems  with  MAX_CRASHING  = 3 were  created  in  a similar 
way.  The  value  of  yi(m)c  in  the  i-th  problem  with  MAX  CRASHING  = Default  * .5  was 
created  by  setting  yj(m)c  as  the  average  between  yi(m)n  and  yi(m)c  of  the  i-th  default  problem. 
If  the  modified  value  of  yi(m)c  is  not  an  integer,  the  fractional  portion  was  dropped.  The 
penalty  cost,  the  due  date  and  the  resource  capacities  were  recalculated. 
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The  i-th  problem  with  | M(i)  | = 1 was  created  by  dropping  the  second  mode 
of  the  i-th  default  problem.  The  i-th  problem  with  | M(i)  | = 3 was  created  by 
generating  and  adding  the  third  mode  to  the  i-th  default  problem.  In  both  cases,  the 
penalty  cost,  the  due  date  and  the  resource  capacities  were  recalculated. 

The  last  three  rows  in  the  table  show  the  different  settings  for  the  penalty,  the  due 
date  and  the  resource  capacities.  These  problems  are  different  from  the  default  problems 
in  that  the  corresponding  parameter  in  the  problem  is  found  by  multiplying  the  parameter 
in  the  default  problem  by  .8,  1.2  or  1.4.  If  capacities  were  reduced  to  80  percent  of  the 
default  capacities  and  the  resulting  capacity  of  resource  type  k,  bk,  is  less  than  the  largest 
ri(m)k  in  the  problem,  then  bk  is  set  to  the  largest  ri(m)k  in  the  problem. 

6.3  Report  of  the  Computational  Results 

In  order  to  increase  readability,  our  exact  solution  and  heuristic  solution 
procedures  are  named  E.A.E.  (Erenguc  and  Ahn’s  Exact  Solution  Procedure)  and  E.A.H. 
(Erenguc  and  Ahn’s  Heuristic  Procedure),  respectively.  Since  the  backtracking 
procedure  is  used  as  an  exact  solution  and  as  a heuristic  procedure,  the  backtracking 
procedure  as  an  optimal  solution  procedure  and  the  backtracking  procedure  as  a heuristic 
procedure  are  named  B.E.  (the  backtracking  algorithm  as  an  exact  solution  procedure) 
and  B.H.  (the  backtracking  algorithm  as  a heuristic  procedure),  respectively.  When  we 
performed  some  preliminary  experiments  on  the  efficiency  of  B.H.,  we  found  out  that 
the  best  solution  found  by  B.H.  often  deviated  significantly  from  the  optimal  solution. 
Obviously,  B.H.  cannot  search  enough  of  the  feasible  region  in  a limited  time.  This  led 
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us  to  consider  a modification  of  B.H.  The  modification  is  that  only  the  normal  and  the 
crash  durations  are  searched  by  B.H.  This  modification  can  be  justified  as  follows:  if 
we  do  not  have  enough  time  to  investigate  all  the  possible  cases,  then  we  can  search  only 
the  most  representative  cases.  In  order  to  distinguish  B.H.  from  the  modified  B.H.,  we 
will  denote  the  modified  B.H.  as  B.H.M.  One  thing  to  be  mentioned  is  the  application 
of  the  bounding  rule,  left-shift  rule,  suggested  by  Arno  Sprecher  [1994J.  This  rule  was 
already  applied  successfully  by  Demeulemeester  et  al.  [1992]  for  the  singe  mode 
makespan  minimization  problem.  This  rule  prevents  further  searching  of  an  inferior 
branch  at  an  early  stage.  Although  this  rule  improves  the  performance  of  B.E.,  the 
effect  of  this  rule  for  B.H.  is  not  significant.  The  problem  with  the  application  of  this 
rule  is  that  the  search  can  be  terminated  by  the  imposed  time  limit  before  the  search 
reaches  the  dominant  branch.  As  a result,  in  a few  problems,  no  feasible  solution  was 
found  within  the  time  limit.  In  contrast  to  B.H.,  this  rule  showed  a considerable 
contribution  to  B.H.M.  In  order  to  maintain  consistency,  we  applied  this  rule  for  both 
B.H.  and  B.H.M,  and  if  no  solution  was  found  within  time  limit,  then  we  dropped  this 
rule  to  find  a feasible  solution.  Also  our  preliminary  experiments  on  E.A.E.  revealed 
that  checking  the  resource  constraints  before  checking  noninteger  solutions  or  the 
existence  of  an  underestimation  gap  is  the  most  efficient  in  terms  of  CPU  time.  The 
modified  Step  3 of  Procedure  4.27  (E.A.E.  Procedure)  is  described  in  Appendix  B. 

All  the  programs  are  coded  in  C.  An  IBM  mainframe  (IBM  ES/9000-831  with 
3 vector  facilities  running  MVS/ESA)  was  used  for  E.A.E.  and  B.E.  with  a maximum 
execution  time  of  3,600  seconds.  The  computer  used  for  E.A.H.,  B.H.  and  B.H.M.  was 
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an  IBM  compatible  PC  with  80386dx  processor  with  no  math  co-processor.  We  set  the 
maximum  execution  time  for  the  heuristic  procedures  to  (I  * I * 0.3),  where  I is  the 
number  of  the  nondummy  activities  in  the  project.  This  time  limit  setting  forces  each 
heuristic  procedure  to  solve  large  size  problems  in  a polynomially  bounded  time.  In 
E.A.H.,  the  number  of  activities  to  be  reinitialized  in  Step  4 of  Procedure  5.13  was  set 
to  (I  * 0.1).  Since  the  execution  time  limit  was  imposed  on  each  heuristic  procedure, 
the  maximum  number  of  trials  in  Procedure  5.13  was  set  to  infinity,  and  the  procedure 
was  terminated  when  the  maximum  execution  time  limit  was  reached. 

6.3.1  Comparison  of  the  Exact  Solution  Procedures 

In  order  to  investigate  the  efficiency  of  B.E.  and  E.A.E.,  two  sets  of  problems, 
with  MAXCRASHING  = 2 and  with  MAXCRASHING  = 3,  respectively,  were 
solved  by  the  procedures.  The  number  of  nondummy  activities  in  both  problem  sets  was 
10.  The  result  is  reported  in  Table  6.5. 


Table  6.5  Comparison  of  the  Execution  Times  of  B.E.  and  E.A.E. 


Prob  No 

B.E. 

E.A.E. 

B.E. /E.A.E. 

1 

1,042 

25 

41.68 

over  3,600 

34 

N/A 

2 

138 

9 

15.33 

561 

7 

80.14 

3 

554 

14 

39.57 

2,727 

11 

247.91 

4 

226 

2 

113 

1,960 

3 

653.33 
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Each  exact  solution  procedure’s  consumed  CPU  time  (in  seconds)  for  the  solving 
i-th  problem  is  reported.  In  the  second  and  the  third  columns,  the  top  number  indicates 
the  CPU  time  for  the  problem  with  MAX  CRASHING  = 2,  and  the  bottom  one 
indicates  the  CPU  time  for  the  problem  with  MAX_CRASHING  = 3.  The  last  column 
shows  the  comparison  factor,  which  is  (the  CPU  time  of  B.E.  / the  CPU  time  of 
E.A.E.).  Since  B.E.  failed  to  optimally  solve  the  first  problem  with  MAX_CRASHING 
= 3 in  3,600  seconds,  the  comparison  factor  for  that  case  is  listed  as  N/A  in  the  last 
column.  When  there  are  only  3 possible  durations  for  each  mode  (MAX_CRASHING 
= 2),  the  comparison  factor  ranges  from  15.33  to  113  as  shown  in  the  last  column. 
When  the  number  of  possible  durations  increases  by  one,  the  comparison  factor  ranges 
from  80.14  to  653.33  excluding  the  first  problem,  which  B.E.  failed  to  find  the  optimal 
solution  within  the  time  limit.  This  clearly  indicates  that  the  computational  effort  for 
solving  the  problem  by  B.E.  is  far  more  sensitive  to  the  number  of  possible  durations  for 
each  mode  than  E.A.E.  is. 

6.3.2  Exact  Solution  Procedure’s  Execution  Time  and  Comparison  of  the  Heuristics 
Since  our  preliminary  experiment  showed  the  inadequacy  of  B.E.  as  an  exact 
solution  procedure  as  shown  in  section  6.3.1,  no  further  effort  was  made  to  optimally 
solve  the  problems  with  B.E.  Thus,  only  the  computational  results  of  E.A.E.  are 
reported  here.  The  purpose  of  solving  problems  with  E.A.E.  is  twofold:  testing  the 
efficiency  of  E.A.E.  on  the  various  problems  and  providing  the  optimal  solutions  that 
will  be  used  for  the  efficiency  tests  of  B.H.,  B.H.M  and  E.A.H.  The  computational 
results  for  the  problems  with  I = 10  are  listed  in  Table  6.6. 
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Table  6.6  Computational  Results  on  Problems  with  I = 10 


I 

Problem 

Time 

Limit 

B.H. 

[min,  max  error] 
[avg,  std  error] 

B.H.M 

[min,  max  error] 
[avg,  std  error] 

E.A.H. 

[min,  max  error] 
[avg,  std  error] 

E.A.E 

[min,  max  sec] 
[avg,  std  sec] 

10 

Default 

30 

[63.17,  118.92] 
[83.94,  22.54] 

[1.54,  52.21] 
[27.48, 18.14] 

[0.00,  2.07] 
[0.88,0.91] 

[4,  67] 

[29.50,  25.85] 

10 

MAX  CRASHING 
= 2 

30 

[16.97,  67.70] 
[43.81,  19.80] 

[6.03,41.10] 
[19.05,  14.40] 

[0.00,  2.94] 
[1.23,  1.12] 

[2,  25] 
[12.50,8.38] 

10 

MAX  CRASHING 
= 3 

30 

[22.99,  85.29] 
[52.96,  22.69] 

[5.42,41.46] 
[19.76,  14.57] 

[0.00,  4.20] 
[1.51,  1.65] 

[3,  34] 

[13.75,  12.03] 

10 

MAX  CRASHING 
* .5 

30 

[41.44,94.11] 
[72.84,  20.03] 

[1.34,  47.82] 
[20.16,  18.43] 

[0.00,9.13] 
[2.70,  3.77] 

[2,  29] 

[15.25, 10.76] 

10 

COMPLEXITY 
= 1.5 

30 

[7.94,  65.13] 
[37.25,  20.34] 

[3.83,  29.03] 
[13.53,  9.45] 

[0.34,  2.92] 
[1.58,0.99] 

[2,  46] 

[24.25, 15.56] 

10 

COMPLEXITY 
= 2.1 

30 

[10.37,  102.96] 
[57.18,35.54] 

[2.92,  42.00] 
[23.25,  13.88] 

[0.00,  2.23] 
[0.68,  0.92] 

[2,24] 

[10.75,8.10] 

10 

1 M(i)  | = 1 

30 

[14.11,32.43] 
[27.02,  7.54] 

[0.63,  3.95] 
[2.45,  1.32] 

[0.00,  0.99] 
[0.38,0.41] 

[2,  41] 

[12.25, 16.60] 

10 

1 M(i)  | = 3 

30 

[37.55,  94.72] 
[63.65,  23.45] 

[3.97,  42.59] 
[28.58,  15.59] 

[1.26,2.06] 

[1.54,0.32] 

[4,  74] 

[30.25,26.93] 

10 

PENALTY  * .8 

30 

[44.15,95.51] 
[62.00,  20.62] 

[0.85,41.85] 
[20.75,  14.92] 

[0.00,  2.96] 
[1.60,  1.07] 

[4,  74] 

[32.00,28.70] 

10 

PENALTY  * 1.2 

30 

[74.96,  137.43] 
[99.39,24.80] 

[2.05,  58.84] 
[30.71,20.19] 

[0.00,  0.28] 
[0.10,0.11] 

[4,  70] 

[31.50,  26.22] 

10 

DUEDATE*  .8 

30 

[51.95,  101.51] 
[72.71,20.38] 

[7.80,45.41] 
[25.14,  13.57] 

[0.00,  3.60] 
[1.09,  1.46] 

[7,  80] 

[36.50,  30.76] 

10 

DUEDATE*  1.2 

30 

[59.48,  135.39] 
[90.84,28.89] 

[0.00,56.98] 
[26.63,  20.23] 

[0.00,  1.84] 
[0.46,  0.80] 

[0,  47] 

[23.50,  22.05] 

10 

CAPACITY  * .8 

30 

[53.74,  108.00] 
[80.81,24.32] 

[8.13,54.31] 
[36.22, 18.43] 

[0.00,  0.84] 
[0.21,0.36] 

[1,37] 

[23.50, 13.50] 

10 

CAPACITY  * 1.2 

30 

[27.50,  98.04] 
[67.46,  29.96] 

[3.93,  86.63] 
[32.26,  32.87] 

[0.00,4.67] 

[2.33,2.16] 

[4,  100] 
[33.50,  39.07] 

10 

CAPACITY  * 1.4 

30 

[21.57,47.12] 
[33.87,  10.40] 

[0.82,  29.51] 
[11.90,  10.88] 

[0.00,  0.67] 
[0.21,0.28] 

[0,  45] 

[15.00, 17.76] 
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The  first  and  the  second  column  show  the  number  of  nondummy  activities  and  the 
problem  set  characteristics,  respectively.  Each  problem  set  consists  of  four  problems. 
The  third  column  lists  the  imposed  time  limit  on  the  heuristic  procedures.  In  the  fourth 
column  to  the  sixth  column,  the  first  line  shows  the  minimum  error  and  the  maximum 
error  of  the  heuristic  values,  and  the  second  line  shows  the  average  error  and  the 
standard  deviation  of  the  errors.  All  the  values  associated  with  errors  are  percentages. 
Error  is  calculated  using  the  formula,  error  = 1 - (the  heuristic  value  -optimum  value) 
/ optimum  value,  where  the  optimum  value  is  found  by  E.A.E.  In  the  last  column,  the 
relevant  data  for  CPU  time  for  solving  each  problem  by  E.A.E.  are  listed.  In  the  first 
line  the  minimum  and  maximum  CPU  times  are  reported,  and  in  the  second  line  the 
average  CPU  time  and  the  standard  deviation  of  CPU  times  are  listed.  The  unit  of 
measure  for  the  CPU  times  is  second.  When  we  recorded  the  CPU  time  for  each 
problem,  we  used  an  integer  variable.  Thus,  each  CPU  time  is  an  integer  value  after 
truncating  the  fractional  portion.  Hence,  x seconds  of  CPU  time  in  the  table  indicates 
that  the  actual  CPU  time  is  equal  to  or  greater  than  x seconds  but  less  than  (x+1) 
seconds.  The  average  and  standard  deviations  of  the  CPU  times  are  based  on  the  integer 
values  of  CPU  times. 

15  sets  of  problems,  thus  a total  of  60  problems,  were  solved.  The  average  CPU 
time  for  E.A.E.  was  22.93  seconds  with  a maximum  of  100  seconds,  a minimum  0 
second  and  a standard  deviation,  23.67  seconds. 

Table  6.7  summarizes  the  efficiency  of  each  heuristic  procedure.  The  average 
errors  of  B.H.,  B.H.M  and  E.A.H.  were  63.06,  22.53  and  1.10,  respectively.  The  last 
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row  shows  the  number  of  times  the  corresponding  heuristic  procedure  found  the  optimum 
solution  value.  While  B.H.  and  B.H.M.  found  no  and  one  optimal  solution  respectively, 
E.A.H.  found  23  optimal  solutions  among  the  60  problems.  It  is  worthwhile  to  mention 
that  in  all  60  problems,  the  solution  value  found  by  E.A.H.  was  the  best  and  the  solution 
value  found  by  B.H.M.  was  the  second  best  except  in  only  one  problem,  where  B.H.M. 
and  E.A.H.  both  found  the  optimum  solution  value.  In  all  60  problems,  B.H.  performed 
the  worst.  Table  6.8  summarizes  the  distribution  of  errors  for  E.A.H.  for  the  problems 
with  1 = 10. 


Table  6.7  Efficiency  of  Each  Heuristic  Procedure  for  the  (1  = 10)  Problems 


Error  (%) 

B.H. 

B.H.M 

E.A.H. 

Minimum 

7.94 

0.00 

0.00 

Maximum 

137.43 

86.63 

9.13 

Average 

63.06 

22.53 

1.10 

STD 

30.958 

19.04 

1.60 

Optimum 

Found 

0 

1 

23 

Table  6.8  Distribution  of  Errors  of  E.A.H.  for  the  Problems  with  1 = 10 


Error  (%) 

0 

0 < x < 1 

1 < x < 2 

2 < x < 3 

x > 3 

Frequency 

23 

15 

10 

7 

5 

Table  6.6  shows  that  the  problem  sets  that  were  relatively  quickly  solved  by 
E.A.E.  were  solved  by  B.H.  and  B.H.M.  generally  with  small  errors.  The  positive 
correlation  among  the  performances  of  the  heuristic  procedures  based  on  the  backtracking 


126 


procedure  and  the  CPU  time  of  E.A.E.  can  be  understood  from  the  fact  that  both 
heuristic  procedures  are  based  on  B.E,  which  is  an  optimal  solution  procedure.  The  lack 
of  clear-cut  correlation  between  E.  A. H.  and  any  other  procedure  can  be  understood  from 
the  fact  that  E.A.H.  is  not  a truncated  optimal  solution  procedure.  Generally,  the 
problems  that  can  be  characterized  as  one  of  the  following  were  solved  by  E.A.E.  with 
small  CPU  times  and  solved  by  B.H.  and  B.H.M.  with  small  errors:  small  intervals 
between  the  mode  normal  and  crash  durations,  high  or  low  complexity,  small  number  of 
modes  available  to  each  activity,  high  or  low  capacity.  The  other  factors  such  as  penalty 
cost  and  due  date  affected  the  performances  of  B.H.  and  B.H.M.  and  the  performance 
of  E.A.E.  in  an  opposite  way  in  that  high  penalty  cost  and  long  due  date  worsened  the 
errors  of  B.H.  and  B.H.M.  while  high  or  low  penalty  costs  and  short  due  dates  worsened 
the  CPU  times  of  E.A.E.  The  relationship  between  the  performance  of  E.A.H.  and  the 
variation  of  the  parameters  was  not  systematic. 

The  computational  results  with  I = 12  and  I = 14  are  shown  in  Table  6.9. 
Among  the  four  problems  with  I = 14,  only  two  problems  were  optimally  solved  by 
E.A.E  within  the  maximum  time  limit  of  3,600  seconds.  The  first  problem  of  (1  = 14) 
was  solved  with  an  error  tolerance  of  10  percent,  which  implies  that  the  best  solution 
found  by  E.A.E.  is  deviated  from  the  optimal  solution  at  most  by  10  percent  of  the 
optimal  solution.  In  the  fourth  problem  of  (1  = 14),  the  error  tolerance  applied  was  2 
percent.  In  calculating  the  errors  of  each  heuristic  procedure,  the  lower  bound  was  used, 

where  the  lower  bound  was  found  as  follows: 

, , , , the  best  solution  found  by  E.A.E. 

the  lower  bound  = 

(1  + error  tolerance) 
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The  computational  results  on  the  problems  with  I > 16  are  reported  in  Table 

6.10. 


Table  6.9  Computational  Results  on  Problems  with  I = 12  and  I = 14 


I 

Problem 

Number 

Time 

Limit 

B.H. 

Error  (%) 

B.H.M 
Error  (%) 

E.A.H 
Error  (%) 

E.A.E 

CPU  (second) 

12 

1 

43 

47.22 

38.31 

0.79 

11 

12 

2 

43 

44.89 

65.79 

6.54 

1,331 

12 

3 

43 

41.41 

18.69 

2.40 

163 

12 

4 

43 

61.83 

35.56 

2.48 

37 

14 

1 

58 

47.58  * 

36.86  * 

11.11  * 

108  * 

14 

2 

58 

69.18 

40.97 

4.95 

375 

14 

3 

58 

48.38 

37.07 

5.04 

1,543 

14 

4 

58 

25.73  ** 

18.40  ** 

5.33  ** 

2,377  ** 

*(**)  indicates  that  the  error  tolerance  applied  on  E.A.E.  is  10  (2)  percent 


In  every  problem  with  I > 12,  E.A.H.  performed  the  best  and  B.H.  performed 
the  worst.  For  the  problems  with  the  number  of  nondummy  activities  greater  than  or 
equal  to  16,  only  heuristic  procedures  were  used  to  solve  the  problems.  In  order  to 
compare  the  heuristic  procedures,  the  relative  deviation  concept  is  employed  in  which  a 
relative  deviation  of  a solution  is  the  percentage  deviation  from  the  best  known  heuristic 
solution.  Table  6.10  shows  the  performances  of  the  three  heuristic  procedures.  For 
problems  with  I = 16  or  more,  the  average  deviations  of  B.H.  and  B.H.M.  from  the 
solution  values  of  E.A.H.  was  56.13  percent  and  31.21  percent,  respectively. 
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Table  6.10  Computational  Results  on  Problems  with  I > 16 


I 

Problem 

Time 

B.H 

B.H.M 

E.A.H. 

Limit 

(Solution) 

(Solution) 

(Solution) 

(Deviation  %) 

(Deviation  %) 

(Deviation  %) 

16 

1 

76 

1,760 

1,663 

1,130 

55.75 

47.17 

0.00 

16 

2 

76 

1,561 

1,378 

1,035 

50.82 

33.14 

0.00 

16 

3 

76 

1,555 

1,380 

1,086 

43.19 

27.07 

0.00 

16 

4 

76 

1,439 

1,282 

902 

59.53 

42.13 

0.00 

20 

1 

120 

1,555 

1,444 

1,175 

32.34 

22.89 

0.00 

20 

2 

120 

2,213 

1,457 

1,097 

101.73 

32.82 

0.00 

20 

3 

120 

1,590 

1,272 

1,198 

32.72 

6.18 

0.00 

20 

4 

120 

1,475 

1,418 

1,058 

39.41 

34.03 

0.00 

25 

1 

187 

2,287 

1,601 

1,213 

88.54 

31.99 

0.00 

25 

2 

187 

2,004 

1,659 

1,439 

39.26 

15.29 

0.00 

25 

3 

187 

2,415 

2,224 

1,393 

73.37 

59.66 

0.00 

25 

4 

187 

2,344 

1,825 

1,494 

56.89 

22.16 

0.00 

CHAPTER  7 

CONCLUSIONS  AND  FUTURE  RESEARCH 


In  this  research  we  investigated  the  resource  constrained  project  scheduling 
problem  with  time/cost  trade-offs,  where  multiple  modes  for  each  activity  are  assumed. 
While  research  dealing  with  resource  constrained  project  scheduling  problems  has 
concentrated  on  problems  where  activity  duration  is  solely  determined  by  the  mode  to  be 
used,  we  allowed  for  crashing  of  activity  durations  through  higher  direct  cost.  As  far 
as  we  know,  there  is  no  work  that  considered  the  resource  constrained  project  scheduling 
problem  with  crashable  mode  duration  for  the  nonpreemptive  case.  We  believe  our 
approach  is  the  first  realistic  one  that  bridges  the  gap  between  the  time/cost  trade-off 
problem  and  the  resource  constrained  project  scheduling  problem. 

If  we  consider  each  duration  of  a mode  as  an  individual  mode,  then  this  problem 
can  be  solved  by  the  backtracking  procedure  proposed  by  Patterson,  Slowinski,  Talbot 
and  Weglarz  and  accelerated  by  Sprecher.  But  its  computational  requirement,  both  as 
an  exact  solution  and  as  a heuristic  procedures,  is  enormous.  Thus,  we  presented  an 
efficient  exact  solution  procedure  and  a heuristic  procedure. 

Our  exact  solution  procedure  is  based  on  a branch  and  bound  technique.  Since 
the  activity  cost  function  is  nonlinear,  three  activity  cost  functions  are  defined  in  order 
to  construct  relaxed  (bounding)  problems.  The  relaxed  problem  in  the  branch  and  bound 
tree  is  constructed  after  dropping  the  resource  constraints  and  replacing  the  activity  cost 
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function  with  the  convex  activity  cost  function.  The  relaxed  problem  is  a linear  program. 
Once  the  relaxed  problem  is  solved,  then  a series  of  fathoming  tests  are  performed.  If 
fathoming  tests  fail,  branching  rules  are  employed  in  order  to  partition  the  feasible 
region.  As  a result,  we  will  have  tighter  relaxed  problems,  and  the  chance  of  finding  a 
solution  will  be  increased. 

We  proposed  a multi-pass  heuristic  procedure.  In  each  pass,  we  assign  a mode 
and  duration  to  each  activity  and  apply  the  minimum  latest  finish  time  rule  in  order  to 
construct  an  initial  feasible  schedule.  We  assign  the  mode  and  duration  to  each  activity 
such  that  the  initial  assignment  of  mode  and  duration  differs  from  the  initial  assignment 
of  the  previous  pass.  Although  the  same  initial  assignment  can  repeat  in  the  next  pass, 
the  vast  number  of  the  combinations  of  mode  and  durations  renders  a repeat  almost 
impossible.  Once  a feasible  schedule  is  made,  then  four  types  of  improvements  are 
applied  until  no  more  improvements  are  possible.  Then,  the  best  solution  found  so  far 
(incumbent)  and  the  solution  found  in  the  current  pass  are  compared,  and  the  incumbent 
is  updated  if  necessary.  At  the  end  of  each  pass,  reinitialization  of  the  initial  mode  and 
duration  is  done  by  changing  the  mode  and  the  duration  of  a predetermined  number  of 
activities.  First,  we  find  the  activity  whose  mode  has  not  changed  recently.  Then,  we 
assign  the  mode  to  that  activity,  and  the  duration  is  randomly  generated  between  the 
mode  normal  and  crash  durations.  The  activities  not  chosen  for  the  reinitialization  keep 
the  mode  and  duration  assignments  in  the  schedule  of  the  pass.  A new  feasible  schedule 
is  constructed  with  the  mode  and  duration  assignments  provided  from  the  previous  pass 
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using  the  minimum  latest  finish  time  rule,  and  the  whole  process  is  repeated  until  a 
predetermined  number  of  passes  is  tried  or  until  the  maximum  execution  time  is  reached. 

In  order  to  test  the  efficiency  of  our  exact  solution  procedure  and  heuristic 
procedure,  we  generated  a total  of  80  problems  and  solved  them  by  our  procedures  and 
by  the  backtracking  procedure  as  an  exact  solution  procedure  and  as  a heuristic 
procedure.  Our  initial  experiment  showed  that  the  backtracking  procedure  as  an  exact 
solution  method  requires  too  much  computational  effort.  When  the  number  of  problems 
to  be  solved,  the  number  of  nondummy  activities  in  the  problem,  the  number  of  modes, 
and  the  possible  duration  of  each  mode  are  4,  10,  2,  and  4,  respectively,  the  comparison 
factor  of  the  backtracking  procedure  over  our  exact  solution  procedure  for  3 problems 
ranges  from  80.14  to  653.33  with  respect  to  CPU  time  on  a mainframe.  The 
backtracking  procedure  failed  to  optimally  solve  one  of  the  4 problems  in  3,600  seconds. 
Our  exact  solution  procedure  solved  all  the  problems  where  the  number  of  nondummy 
activities  is  equal  to  or  less  than  12  within  the  maximum  time  limit  of  3,600  seconds. 

We  used  a truncated  backtracking  procedure  as  a heuristic  procedure.  Since  our 
initial  experiment  revealed  relatively  a wide  deviation  of  the  solutions  of  the  procedure 
from  the  optimal  solution,  we  devised  a modified  backtracking  procedure,  which 
considered  only  the  mode  normal  and  crash  durations.  We  tested  the  three  heuristic 
procedures  on  the  problems  with  a 386dx  PC.  The  maximum  execution  time  was  set  to 
(I  * I * 0.3)  seconds,  where  I is  the  number  of  nondummy  activities  in  the  problem.  On 
79  of  the  80  problems,  our  heuristic  procedure  was  the  best,  the  modified  backtracking 
procedure  was  the  second  and  the  backtracking  procedure  was  the  worst.  In  one 
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problem,  both  our  heuristic  procedure  and  the  modified  backtracking  procedure  found 
the  optimal  solution.  For  the  60  problems  with  10  nondummy  activities,  the  average 
deviation  of  each  heuristic  solution  from  the  optimal  solution  were  1 . 1 percent  for  our 
heuristic  procedure,  22.53  percent  for  the  modified  backtracking  procedure,  and  63.06 
percent  for  the  backtracking  procedure. 

For  further  research,  we  plan  to  work  on  the  NPV  maximization  problem  in  the 
resource  constrained  multi  mode  project  scheduling  with  crashable  mode  duration. 


APPENDIX  A 

NECESSARY  MODIFICATION  OF  THE  BACKTRACKING  PROCEDURE 


The  backtracking  algorithm  for  solving  a general  class  of  precedence  and 
resource-constrained  scheduling  problems  can  be  used  to  solve  our  problem  after  proper 
modification.  Since  the  algorithm  is  explained  by  Sprecher  [1994]  in  detail,  we  will 
present  only  the  necessary  modification. 

The  procedure  assumes  that  each  mode  of  an  activity  has  a unique  duration  and 
that  a set  of  renewable  resources  and  a fixed  amount  of  nonrenewable  resource  are 
required  to  process  the  activity  with  the  mode.  Further,  a fixed  project  due  date  that 
cannot  be  violated  is  assumed.  In  our  problem,  each  mode  of  an  activity  requires  a set 
of  resources  to  process  the  activity,  but  the  execution  time  of  an  activity  with  a mode  can 
be  shortened  from  the  mode  normal  duration  up  to  the  mode  crash  duration  through 
higher  direct  cost.  Thus,  we  modify  each  duration  of  an  activity  with  a mode  as  an 
individual  mode.  Each  newly  modified  mode  assumes  the  same  resource  requirements 
as  the  original  mode.  The  cost,  which  is  equivalent  to  the  nonrenewable  resource  in  the 
backtracking  procedure,  of  the  newly  modified  mode  is  based  on  the  original  mode 
normal  cost  plus  the  appropriate  additional  cost  for  the  crashing.  Let  us  denote  the  cost 
of  the  newly  defined  mode  as  COST(i,mode). 

Suppose  that  activity  i with  mode  m is  considered  for  assignment  in  the  original 
backtracking  procedure.  The  assignment  is  checked  through  two  tests.  In  the  first  test, 
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the  procedure  checks  whether  the  activity  i with  mode  m can  be  finished  no  later  than 
its  latest  finish  time  while  satisfying  the  precedence  and  resource  constraints  and  starting 
no  earlier  than  the  start  time  of  any  already  scheduled  activity.  (Bounding  on  the  start 
time  no  earlier  than  the  start  time  of  any  already  scheduled  activity  is  justified  by  the 
precedence  tree  concept.)  If  possible,  the  second  test  checks  whether  the  lower  bound 
of  the  current  partial  schedule  is  less  than  the  best  known  solution.  Notice  that  the  first 
test  is  based  on  the  latest  finish  time  that  requires  a nonviolatible  fixed  project  due  date. 
Since  the  project  due  date  can  be  violated  in  our  problem,  the  first  test  is  modified  by 
finding  the  earliest  finish  time  while  satisfying  the  precedence  and  resource  constraints 
and  starting  no  earlier  than  the  start  time  of  any  already  scheduled  activity.  In  the 
second  test,  the  lower  bound  of  the  objective  function  is  found  as  follows: 

£ COST(j,m*)+  £ LB( COST(j) ) +COST(i,m) 

jeAS  j«AS,  j*  i 

+ Penalty  * max  (LB(FI+1)  - DUEDATE,  0}, 
where  AS  is  a set  of  already  scheduled  activities  and  LB(.)  denotes  the  lower  bound  of 
function  (.).  Note  that  our  objective  function  is  the  sum  of  each  activity  cost  and  the 
penalty  cost.  Thus,  each  scheduled  activity  cost  follows  COST(j,m*),  where  m*  is  the 
mode  chosen  for  activity  j.  The  lower  bound  of  the  unscheduled  activity  cost, 
LB(COST(j)),  j ^ i,  is  based  on  the  lowest  COST(j,m)  among  the  available  modified 
modes  of  activity  j.  The  cost  of  activity  i is  COST(i,m)  since  activity  i with  mode  m is 
under  consideration  for  scheduling.  The  lower  bound  of  the  project  completion  time  is 
based  on  the  scheduled  activities’  finish  time,  the  earliest  finish  time  of  activity  i found 
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in  the  first  test  and  the  earliest  finish  time  of  unscheduled  activities.  When  we  consider 
the  earliest  finish  time  of  an  unscheduled  activity,  each  duration  of  the  activity  is  based 
on  the  shortest  mode  duration  of  the  activity.  Each  unscheduled  activity’s  start  time  must 
consider  all  the  precedence  constraints  and  must  not  be  earlier  than  the  start  time  of 
activity  i. 

Among  the  five  bounding  rules  proposed  by  Sprecher,  only  two  rules  are  used. 
The  first  bounding  rule  is  based  on  the  latest  finish  time,  which  requires  a nonviolatible 
project  due  date,  and  the  fifth  bounding  rule  is  applied  only  when  there  is  no 
nonrenewable  resources,  which  is  equivalent  to  the  activity  cost.  Thus,  two  bounding 
rules  are  dropped  in  the  modified  procedure.  The  fourth  bound  rule  is  based  on  the 
nonrenewable  resources,  which  we  already  incorporated  in  the  second  test  for  the 
assignment.  Therefore,  only  the  second  and  third  bound  rules  are  used  in  the  modified 
procedure.  The  second  bounding  rule  is  designed  to  avoid  the  construction  of  the  same 
scheduling  (all  the  same  start  time,  mode  but  different  assignment  order),  and  the  third 
bounding  rule  is  a left-shift  rule,  which  is  already  used  for  the  single  mode  case  by 


Demeulemeester  et  al.  [1992]. 


APPENDIX  B 

MODIFIED  STEP  3 OF  THE  BRANCH  AND  BOUND  PROCEDURE 

Step  3.  Testing  whether  we  can  derive  an  optimal  solution  of  (Pp)  from  the 

solution  of  (LPP).  If  we  can,  then  we  have  solved  (Pp).  If  we  cannot, 
branch  from  node  p. 

3.1  Let  Zj*  and  F**  be  the  solution  value  for  activity  i from  the  optimal 
solution  of  (LPP).  Let  Zj**  be  z*  for  i = 1,  I.  For  i = 1 to  I,  find 
m*  such  that  C^.^mm  {z**,  UB(z,)p,  yi(m.)n})  is  the  minimum  among 
m G M(i)p.  If  ties  happen,  choose  one  arbitrarily.  If  z^*  > min 
{UB(zi)p,  yi(m*)n},  then  let  z-**  = min  {UB^,  yi(m,)n}. 

3.2  Resource  Constraint  Checking 

With  z**,  Fj*  and  m*,  check  whether  the  current  schedule  satisfies  all  the 
resource  constraints.  If  there  is  no  resource  conflict  at  any  period,  then 
go  to  Step  3.3.  If  there  is  a resource  conflict  during  some  period,  find  the 
earliest  period  when  a resource  conflict  occurs.  Create  an  RC,  the  set  of 
all  the  activities  in  progress  during  this  period,  and  create  children  nodes. 
For  the  Simultaneous  Processing  Case: 

Find  all  the  mode  permutations  as  shown  in  Step  1 in  Procedure 
4.23.  For  each  permutation,  test  whether  NEC(k)  < bk  for  all  k 
G K}  as  shown  in  Step  2 in  Procedure  4.23.  If  so,  then 
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MAKECHILDNODE.  Update  M(i)k  = m**  for  all  i G RC, 
where  m**  is  a mode  chosen  in  the  permutation. 

For  the  Nonsimultaneous  Processing  Case: 

Find  all  the  possible  ordered-pairs  of  activities  in  the  RC  as  shown 
in  Step  1 of  Procedure  4.25.  For  each  ordered-pair,  (i,j), 
MAKE  CHILD  NODE.  Update  Hk  = Hp  U (i,j). 

Go  to  Step  1. 

3.3  Testing  for  the  Noninteger  Solution 

From  i = 1,  test  whether  z**  or  Fj*  is  an  integer.  If  z-**  is  not  an 
integer,  then  create  children  nodes. 

i.  MAKE  CHILD  NODE.  Let  UB(Zj)k  be  [z**]. 

ii.  MAKE  CHILD  NODE.  Let  LB(Zi)k  be  [z  **]  + 1. 

Go  to  Step  1. 

If  F*  is  not  an  integer,  then  create  children  nodes. 

i.  MAKE_CHILD_NODE.  Let  UB^  be  [F  *]. 

ii.  MAKE  CHILD  NODE.  Let  LB(Fi)k  be  [¥  *]  + 1 . 

Go  to  Step  1. 

If  all  the  activities  are  checked  and  if  there  is  no  noninteger  solution,  then 
go  to  Step  3.4. 

3.4  Testing  for  the  existence  of  an  Underestimation  Gap 

From  i = 1,  test  whether  CC^fe*)  = C^fe**).  If  CC>>(Zi*)  ^ 
Cl(mt)(p>(z**) , then  create  children  nodes. 
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MAKE  CHILD  NODE.  Let  UB(Zi)k  be  z,*. 
MAKE_CHILD_NODE.  Let  LB(Zj)k  be  z*.  Go  to  Step  1. 

If  all  the  activities  are  checked  and  if  there  is  no  underestimation  gap, 
then  we  have  solved  (Pp).  V(PP)  is  V(LPP).  Update  UB  as  V(PP)  and 
SOLUTION  as  the  solution  of  the  (Pp).  Go  to  Step  1. 

Suppose  that  we  found  m*  and  z**  as  shown  in  modified  Step  3 and  that  no 
constraint  violations  have  been  found  in  modified  Step  3.  Notice  that  z**  is  less  than 
or  equal  to  z*.  Since  the  solution  with  zx * and  F;*  is  precedence  feasible,  the  solution 
with  z**  and  F * is  also  precedence  feasible.  Resource  feasibility  and  integer  restriction 
are  checked  in  Steps  3.2  and  3.3.  Thus,  the  solution  (with  Fj*,  z***  and  m*)  is  feasible 
to  (Pp).  It  is  easy  to  see  that  the  solution  is  an  optimal  solution  of  (Pp)  if  CCi(p)(zi*)  = 
CKm*^\z**)  for  i = 1,  ...,  I.  Suppose  that  we  found  m*  and  z-**  as  shown  in  modified 
Step  3 and  that  a constraint  violation  has  been  found  in  modified  Step  3.  Since  modified 
Step  3 has  the  same  branching  rules  of  (original)  Step  3,  the  validity  of  each  branching 
rule  in  modified  Step  3 follows  that  in  (original)  Step  3. 
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